matlab仿真一階低通濾波器幅頻特性和相頻特性1

/freqs模擬濾波器的頻率響應語法：h=freqs(b,a,w)

[h,w]=freqs(b,a)

[h,w]=freqs(b,a,f)

freqs(b,a)描述：freqs返回一個模擬濾波器的H(jw)的復頻域響應（拉普拉斯格式）

請給出分子b和分母ah=freqs(b,a,w)根據(jù)系數(shù)向量計算返回模擬濾波器的復頻域響應。freqs計算在復平面虛軸上的頻率響應h，角頻率w確定了輸入的實向量，因此必須包含至少一個頻率點。[h,w]=freqs(b,a)自動挑選200個頻率點來計算頻率響應h[h,w]=freqs(b,a,f)挑選f個頻率點來計算頻率響應h例子：找到并畫出下面?zhèn)鬟f函數(shù)的頻率響應

Matlab代碼：

a=[10.41];b=[0.20.31];w=logspace(-1,1);logspace功能：生成從10的a次方到10的b次方之間按對數(shù)等分的n個元素的行向量。n如果省略，則默認值為50。freqs(b,a,w);Youcanalsocreatetheplotwith:h=freqs(b,a,w);

mag=abs(h);

phase=angle(h);

subplot(2,1,1),loglog(w,mag)

subplot(2,1,2),semilogx(w,phase)Toconverttohertz,decibels,anddegrees,use:f=w/(2*pi);

mag=20*log10(mag);

phase=phase*180/pi;算法：freqsevaluatesthepolynomialsateachfrequencypoint,thendividesthenumeratorresponsebythedenominatorresponse:s=i*w;

h=polyval(b,s)./polyval(a,s)一階低通濾波器頻響計算：%**************************************************************************%一階低通濾波器%取fH=1/（2Π*R*C）%s=jw=j*2Π*f%11%RC電路的電壓增益AVH=———————=—————————（傳遞函數(shù)）%1+sRC1+j(f/fH)%11%電壓增益的幅值（模）|AVH|=—————————=————————(幅頻響應)%/——————/————————%\/1+（wCR）^2\/1+（f/fH）^2%VV%%電壓增益的相位角ΦH=-arctg(wRC)=-arctg(f/fH)（相頻響應）%**************************************************************************R=10000;%電阻的值C=0.0000001;%電容的值fZ=1;%H（e^jw）表達式分子的系數(shù)向量fM=[R*C,1];%H（e^jw）分母的系數(shù)向量%logspace功能：生成從10的a次方到10的b次方之間按對數(shù)等分的n個元素的行向量。n如果省略，則默認值為50。w=logspace(-1,5);%設定頻率的變化范圍%freqs根據(jù)系數(shù)向量計算返回一個模擬濾波器的H(e^jw)的復頻域響應（拉普拉斯格式）g=freqs(fZ,fM,w);%求其頻率響應%angle()返回復數(shù)元素的相角，單位為弧度x=angle(g);%求相頻響應%abs()對復數(shù)返回其幅值y=abs(g);%求幅頻響應subplot(2,1,1),%loglog()繪制雙對數(shù)坐標，semilogx()/semilogy()繪制單對數(shù)坐標semilogx(w,x*180/pi)%畫相頻響應的曲線hw=line([0.1,1000],[-45,-45]);set(hw,'LineWidth',1);set(hw,'color','r');xlabel('角頻率/(rad/s)');ylabel('角度/°');gridon;subplot(2,1,2),%畫幅頻響應的曲線semilogx(w,20*log10(y))h=