为啥无法画图啊,有人看看么


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 

结果画图出不来图,我看得到的y值只有一个值而不是随着自变量变化的

你好,给你改了一个点除(矩阵向量之间对应位置相乘就是.*,同理有./.^),然后保存为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

效果

img