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

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

ID:55953564

大小:36.00 KB

页数:4页

时间:2020-06-18

matlab仿真一阶低通滤波器幅频特性和相频特性.doc_第1页
matlab仿真一阶低通滤波器幅频特性和相频特性.doc_第2页
matlab仿真一阶低通滤波器幅频特性和相频特性.doc_第3页
matlab仿真一阶低通滤波器幅频特性和相频特性.doc_第4页
资源描述:

《matlab仿真一阶低通滤波器幅频特性和相频特性.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

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

2、reqs(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),l

3、oglog(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=po

4、lyval(b,s)./polyval(a,s)一阶低通滤波器频响计算:%**************************************************************************%一阶低通滤波器%取fH=1/(2Π*R*C)%s=jw=j*2Π*f%11%RC电路的电压增益AVH=———————=—————————(传递函数)%1+sRC1+j(f/fH)%11%电压增益的幅值(模)

5、AVH

6、=—————————=————————(幅频响应)%/————

7、——/————————%/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功能:生成从1

8、0的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()/

9、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=line([0.1,10000],[-20,-20]);set(h,'LineW

10、idth',1);set(h,'color','r');xlabel('角频率/(rad/s)');ylabel('幅值/dB');gridon;

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。