matlab仿真一阶低通滤波器幅频特性和相频特性

/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)根据系数向量计算返回模拟滤波器的复频域响应。freqs计算在复平面虚轴上的频率响应h，角频率w确定了输入的实向量，因此必须包含至少一个频率点。[h,w]=freqs(b,a)自动挑选200个频率点来计算频率响应h[h,w]=freqs(b,a,f)挑选f个频率点来计算频率响应h例子：找到并画出下面传递函数的频率响应

Matlab代码：

a=[10.41];b=[0.20.31];w=logspace(-1,1);logspace功能：生成从10的a次方到10的b次方之间按对数等分的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=———————=—————————（传递函数）%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）表达式分子的系数向量fM=[R*C,1];%H（e^jw）分母的系数向量%logspace功能：生成从10的a次方到10的b次方之间按对数等分的n个元素的行向量。n如果省略，则默认值为50。w=logspace(-1,5);%设定频率的变化范围%freqs根据系数向量计算返回一个模拟滤波器的H(e^jw)的复频域响应（拉普拉斯格式）g=freqs(fZ,fM,w);%求其频率响应%angle()返回复数元素的相角，单位为弧度x=angle(g);%求相频响应%abs()对复数返回其幅值y=abs(g);%求幅频响应subplot(2,1,1),%loglog()绘制双对数坐标，semilogx()/semilogy()绘制单对数坐标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=