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

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

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

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

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

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

Matlab代碼：

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