上面是u和ap,下面是w,
wn=2*500*pi;
ipson=0.038;
k=4500;
k_t=2400;
k_r=0.9;
N=2;
en=0;
ex=1.5708;
ayy=0.5*(-cos(2*ex)-2*k_r*ex-k_r*sin(2*ex))-0.5*(-cos(2*en)-2*k_r*en-k_r*sin(2*en));
fid=fopen('e:result2\qieshen3.txt','w');
fid1=fopen('e:result2\zhuansu3.txt','w');
for m=0:4
for w=2000:1:10000
s=i*w;
d=w/wn;
G_y=wn^2/(k*(wn^2+s^2+2*ipson*wn*s));
ap=1/((N/2*pi)*ayy*k_t*real(G_y));
fprintf(fid,'%f\r\n',ap);
n=60*w/(N*(2*m*pi+2*pi-2*atan((d^2-1)/2*ipson*d)));
fprintf(fid1,'%f\r\n',n);
end
end
fclose(fid1);
fclose(fid); |