function u = fw(t2,t1,w)
b=pi/t2;
a=b*cot(b*t1);
c=a/b;
y=exp(c*atan(1/c))*(sqrt(2*pi)*b)*(exp(-c*pi-1i*pi.*(w/b))+1)/(2*pi*((a+1i*w).^2)+b^2);
u=y;
end
t2=0.13;
ratio=0.1;
t1=ratio*t2;
w=0:1:400;
u=fw(t2,t1,w);
plot(w,u); hold on
你好,给你改了一个点除(矩阵向量之间对应位置相乘就是.*
,同理有./
和.^
),然后保存为main.m文档,就可以画图了(有帮助望采纳哟,谢谢啦)
t2=0.13;
ratio=0.1;
t1=ratio*t2;
w=0:1:400;
u=fw(t2,t1,w);
plot(w,u,'-'); hold on
function u = fw(t2,t1,w)
b=pi/t2;
a=b*cot(b*t1);
c=a/b;
y=exp(c*atan(1/c))*(sqrt(2*pi)*b)*(exp(-c*pi-1i*pi.*(w/b))+1)./(2*pi*((a+1i*w).^2)+b^2); % 这里改成了一个点除
u=y;
end
效果