untitled3.asv

来自「这个是我数字信号处理这门课的课程设计,文件包括用matlab编写的代码和设计报告」· ASV 代码 · 共 35 行

w1s=2;w1p=1;N=5;

w2s=0.5;w2p=1;
for i=1:3
    if m==1
for i=0:1:30
   w1=0:0.1:i*0.1;
   if i==0
       w2=3.2;
   elseif i==1
       w2(2)=3.1
   elseif i==2
       w2(3)=3.0
   elseif i==3
       w2(4)=2.8
   else
  w2(i+1)=1/w1(i+1);
   end
   
G=1./((j.*w1).^5+3.236.*(j.*w1).^4+5.2361.*(j.*w1).^3+5.2361.*(j.*w1).^2+3.2361.*(j.*w1)+1);
G1=abs(G);

subplot(2,1,1)
plot(w1,G1);
axis([0,3.5,0,1.2])
subplot(2,1,2)
plot(w2,G1)
axis([0,3.5,0,1.2])
pause(1)

end