资源简介
冲击电流发生器的matlab程序,需先在matlab命令窗口输入需要的R/L/C参数,或自行增加参数输入语句,仅作参考
代码片段和文件信息
function[imim1tmtm1tftt]=cjdlfsq_1(RLC) %P262冲击电流发生器FORTRAN改编程序
dtf=0.05;
dtt=0.1;
u=122.5;
nm=280;
i=zeros(999);
t=zeros(999);
for k=1:nm
t(k+nm)=k*dtt;
t(k)=k*dtf;
end
for k=1:nm
t(k+2*nm)=(nm+k)*dtt;
end
af=R/(2.0*L);
w0=1.0/sqrt(L*C);
afd=sqrt(af*af-w0*w0);
w=sqrt(w0*w0-af*af);
if ((af-w0)<0)
for j=1:3*nm
i(j)=u*exp(-af*t(j))*sin(w*t(j))/(w*L);
end
bt=atan(w/af);
tm1=bt/w;
disp([‘计算峰值时间Tm1为: ‘num2str(tm1)‘ μs‘]);
im1=u*exp(-af*bt/w)/sqrt(L/C);
disp([‘计算电流峰值Im1为: ‘num2str(im1)‘ kA‘]);
elseif ((af-w0)==0)
for j=1:2*nm;
i(j)=u*t(j)*exp(-af*t(j))/L;
end
tm1=sqrt(L*C);
im1=u*sqrt(C/L)*exp(-1.0);
else
p1=-af+afd;
p2=-af-afd;
for j=1:2*nm
i(j)=u*(exp(p1*t(j))-exp(p2*t(j)))/(L*(p1-p2));
end
tm1=log(p2/p1)/(p1-p2);
im1=u*(exp(p1*tm1)-exp(p2*tm1))/(L*(p1-p2));
end
im=i(1);
for j=1:nm
if(i(j)>im)
im=i(j);
end
end
disp([‘实际电流峰值Im为:
相关资源
- 高灵敏度GPS接收机MATLAB仿真,附捕获
- 基于MATLAB的质点弹道计算与外弹道优
- 阵列天线的matlab仿真
- MATLAB 经典程序源代码大全
- MATLAB小波软阈值去噪代码33473
- 天线阵的波束形成在MATLAB仿真程序及
- 非线性SVM算法-matlab实现
- 《MATLAB 智能算法超级学习手册》-程序
- 组合导航matlab程序
- 读取txt文件内容matlab代码实现
- Matlab实现基于相关的模板匹配程序
- matlab优化工具箱讲解
- 基于MATLAB的快速傅里叶变换
- 光纤传输中的分布傅立叶算法matlab实
- 基于matlab的图像处理源程序
- matlab 椭圆拟合程序
- 算术编码解码matlab源代码
- optical_flow 光流法 matlab 实现程序
- 引导图像滤波器 Matlab实现
- 分形几何中一些经典图形的Matlab画法
- OFDM系统MATLAB仿真代码
- SVM工具箱(matlab中运行)
- 图像小波变换MatLab源代码
- LU分解的MATLAB实现
- 冈萨雷斯数字图像处理matlab版(第三
- 替代数据法的matlab程序
- 用matlab实现的多站定位系统性能仿真
- 通过不同方法进行粗糙集属性约简m
- k近邻算法matlab实现
- matlab识别系统
评论
共有 条评论