求助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;%假定为1
t=2;%假定t=2
didt=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)')