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

1、freqs模擬濾波器的頻率響應(yīng)語法：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)來計(jì)算頻率響應(yīng)hh, w = freqs(b, a, f) 挑選f個(gè)頻率點(diǎn)來計(jì)算頻率響應(yīng)h例子

2、：找到并畫出下面?zhèn)鬟f函數(shù)的頻率響應(yīng) Matlab代碼： a = 1 0.4 1;b = 0.2 0.3 1;w = logspace(-1, 1);logspace 功能：生成從10的a次方到10的b次方之間按對(duì)數(shù)等分的n個(gè)元素的行向量。n如果省略，則默認(rèn)值為50。freqs(b, a, w);You can also create the plot with: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

3、)To convert to hertz, decibels, and degrees, use:f = w/(2*pi);mag = 20*log10(mag);phase = phase*180/pi;算法：freqs evaluates the polynomials at each frequency point, then divides the numerator response by the denominator response:s = i*w;h = polyval(b,s)./polyval(a,s)一階低通濾波器頻響計(jì)算：%*%一階低通濾波器 % 取fH = 1/（2

4、*R*C）% s=jw=j*2*f % 1 1%RC電路的電壓增益 AVH = = （傳遞函數(shù)）% 1 + sRC 1 + j(f/fH)% 1 1 %電壓增益的幅值（模）|AVH| = = (幅頻響應(yīng))% / /% / 1 + （wCR）2 / 1 + （f/fH）2% V V%電壓增益的相位角 H = - arctg(wRC) = - arctg(f/fH) （相頻響應(yīng)）%*R=10000; %電阻的值 C=0.0000001; %電容的值 fZ=1; %H（ejw）表達(dá)式分子的系數(shù)向量 fM=R*C,1; %H（ejw）分母的系數(shù)向量 %logspace 功能：生成從10的a次方到10的

5、b次方之間按對(duì)數(shù)等分的n個(gè)元素的行向量。n如果省略，則默認(rèn)值為50。w=logspace(-1,5); %設(shè)定頻率的變化范圍 %freqs 根據(jù)系數(shù)向量計(jì)算返回一個(gè)模擬濾波器的H(ejw)的復(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) %畫相頻響應(yīng)的曲線 hw=line(0.1,1000,-45,-45);set(hw,'LineWidth',1);set(hw,'color','r');xlabel('角頻率/ (rad/s)'); ylabel('角度/ °'); grid on; subplot(2,1,2), %畫幅頻響應(yīng)的曲線 semilogx(w,20*log10(y)h=line(0.1,10000,-20,-20);set(h,'LineWidth',1);set(h,&#