求助matlab公式计算,公式中已知s,n,r,L,i 表达式 和α的取值范围,能否画出e'(t)和

求助matlab公式计算,公式中已知s,n,r,L,i 表达式 和α的取值范围,能否画出e'(t)和

%题目说的有问题吧,应该是作e'(t)和α'之间的关系图吧,因为对α做定积分的结果肯定已经不含α了,在α'为某一常数的情况下结果也是一常数.所以应该是作e'(t)和α'之间的关系图把.L=0.8;r=0.2;n=2000;s=3.75*10^(-4);mu0=1;%假定为1t=2;%假定t=2didt=600*cos(pi*100*t)*pi*100;a=-mu0*s*n*r/(2*pi);% syms alpha alpha0% str=(L^2+r^2-2*L*r*cos(alpha))^(-0.5)*cos(alpha+atan(r*sin(alpha)/(L-r*cos(alpha))));% et=a*didt*int(str,alpha,0,alpha0)%积分没有解析解,只能寻求近似数值解alpha0=0.01:0.01:2*pi;%取值范围et=zeros(size(alpha0));for i=1:length(alpha0)alpha=0:0.001:alpha0(i);str=(L^2+r^2-2*L*r*cos(alpha)).^(-0.5).*cos(alpha+atan(r*sin(alpha)./(L-r*cos(alpha))));et(i)=a*didt*trapz(alpha,str);%数值积分的一种endplot(alpha0,et)gridxlabel('alpha0')ylabel('e(t)')

这下我知道了