clc,clear,close all
%% define parameter
i1=21.9;
i2=10.1;
t1=1.3;
t2=1.7;
t3=6;
t4=58;
n=3;
%% count k1 k2
k1=exp(-(t1/t2)*(n*t2/t1).^(1/n));
k2=exp((-t3/t4)*(n*t4/t3).^(1/n));
%% count I
t=-10:.1:10;
I=i1/k1.*(((t./t1).^n)/(1+(t./t1).^n)).*(exp(-t./t2)) +i2/k2*(((t./t3).^n)/(1+(t./t3).^n))*(exp(-t./t4))
%% plot
figure,plot(t,I),grid on ,title('my plot')
matlab运行过程见 录像5.exe
内容全部来源于网络收集，如有侵权，请联系网站删除！
QQ：24596024