How do i shorten this code?
显示 更早的评论
Hi, i know this is a noob-grade question, but i've tried making a code for a low pass butterworth filter magnitude response. Only problem is that i can't seem to compress the code so that i can have multiple lines of factor N into the graph without repetitions.
omc=1;
om=[0:0.1:4];
x=om/omc;
N=1;
y=x.^(2*N);
z=(1+y).^(0.5);
H=1./z;
figure;plot (om,H,'b');
xlabel('frequency in radians/sec');ylabel('Magnitude');
title('Low Pass Butterworth magnitude response');
N=2;
y=x.^(2*N);
z=(1+y).^(0.5);
H=1./z;
hold on;plot(om,H,'r');hold off;
N=3;
y=x.^(2*N);
z=(1+y).^(0.5);
H=1./z;
hold on;plot(om,H,'g');hold off;
N=4;
y=x.^(2*N);
z=(1+y).^(0.5);
H=1./z;
hold on;plot(om,H,'k');hold off;
legend('N=1','N=2','N=3','N=4');
采纳的回答
Cris LaPierre
2018-11-14
omc=1;
om=[0:0.1:4];
[N,om]=meshgrid(1:4,om);
x=om/omc;
y=x.^(2*N);
z=(1+y).^(0.5);
H=1./z;
plot(om,H)
legend('N=1','N=2','N=3','N=4');
xlabel('frequency in radians/sec');ylabel('Magnitude');
title('Low Pass Butterworth magnitude response');
2 个评论
Cris LaPierre
2018-11-14
If colors are important, change your plot command to the following:
h = plot(om,H);
set(h, {'color'}, {'b'; 'r';'g';'k'});
Deepak Aswani
2018-11-14
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Butterworth 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
