matlab中画波特图的某程序的一些疑问b=[0.4914];a=[1,0.9884,1.2386,0.4914];[h w]=freqs(b,a);phaiw=angle(h);subplot(211);plot(w,abs(h) ,'-k','linewidth',0.5,'markersize',3);hold on;xlabel('角频率(rad)');ylabel('幅度');title('幅频响

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/05 03:01:33
matlab中画波特图的某程序的一些疑问b=[0.4914];a=[1,0.9884,1.2386,0.4914];[h w]=freqs(b,a);phaiw=angle(h);subplot(211);plot(w,abs(h) ,'-k','linewidth',0.5,'markersize',3);hold on;xlabel('角频率(rad)');ylabel('幅度');title('幅频响

matlab中画波特图的某程序的一些疑问b=[0.4914];a=[1,0.9884,1.2386,0.4914];[h w]=freqs(b,a);phaiw=angle(h);subplot(211);plot(w,abs(h) ,'-k','linewidth',0.5,'markersize',3);hold on;xlabel('角频率(rad)');ylabel('幅度');title('幅频响
matlab中画波特图的某程序的一些疑问
b=[0.4914];
a=[1,0.9884,1.2386,0.4914];
[h w]=freqs(b,a);
phaiw=angle(h);
subplot(211);
plot(w,abs(h) ,'-k','linewidth',0.5,'markersize',3);
hold on;
xlabel('角频率(rad)');
ylabel('幅度');
title('幅频响应');
subplot(212)
plot(w,unwrap(phaiw) ,'-k','linewidth',0.5,'markersize',3);
hold on;
xlabel('角频率(rad)');
ylabel('相位');
title('相频响应');
这个程序中的w应该是角频率了,那是不是传递函数也必须要H(w)?如果是H(s)需要除以个2π?

matlab中画波特图的某程序的一些疑问b=[0.4914];a=[1,0.9884,1.2386,0.4914];[h w]=freqs(b,a);phaiw=angle(h);subplot(211);plot(w,abs(h) ,'-k','linewidth',0.5,'markersize',3);hold on;xlabel('角频率(rad)');ylabel('幅度');title('幅频响
freqs 函数返回的是 角频率 rad/s,输入的是H(s)的系数,
注意是系数,除什么都无关结果(分子分母消掉了)

b=0.4914;
a=[1,0.9884,1.2386,0.4914];
[h, w]=freqs(b,a);
freqs(b,a)
phaiw=angle(h);
figure
subplot(211);
plot(w,abs(h) ,'-k','linewidth',0.5,'markersize',3);
hold on;<...

全部展开

b=0.4914;
a=[1,0.9884,1.2386,0.4914];
[h, w]=freqs(b,a);
freqs(b,a)
phaiw=angle(h);
figure
subplot(211);
plot(w,abs(h) ,'-k','linewidth',0.5,'markersize',3);
hold on;
xlabel('角频率(rad)');
ylabel('幅度');
title('幅频响应');
subplot(212)
plot(w,unwrap(phaiw) , '-k','linewidth',0.5,'markersize',3);
hold on;
xlabel('角频率(rad)');
ylabel('相位');
title('相频响应');
按上面语句试试,图2可参照第一个改

收起