matlab中想用filter函数实现IIR滤波器

来源：学生作业帮助网编辑：作业帮时间：2024/11/06 13:51:25

x��V_OY�*�f�)0C�q�Ӥ�1�� I�� t�Y�6b�� t�V�m��b��"�~��ޙy�+�s�aDw�]ý��=xS��M 51�i�ٯ��N*�"��Km��3�w��HZ[��&~�ss��!��Ŭ#+��x~2�i4;͹N��5�V��A/��v�]H}��u-��'��r��|aTK��}w�b��B��¼��^;�|d~��jh��s�\$s 6N^�� SAJQ�Y֘DsyԤ�W^ڪ&��91��R~�&s��h��z�u�z��(��M�y]�'[n�h\C#z�7��%2fӫ9cwc�s�1o�ɳ�|h��P��̍R��=y��k�� JG>;��d0�r��j_c09��R�<�C8��~[4��R ��F�tz�鱠��-PCh��Y�Kt(.f�� \X+-bw��T�ń)m\VT5��ܔ&L�BA :�N]-�k�/�I�34�]��w��&CyJcч��Np.ƢB:�� ~�� k�h9��F�<�aHz��0A2��tD��0��q�߉0Z+0��.�1R2�G��3B�G��G�$��ȲӶ*�5��8ce멤|H*��2��/F�5�d��O��y$'�%г��չ�'��o�Ԥ��&��Q.�2��S9�,��uc��-�%b��^x��V��c|��9ܾ�@�v�˩��c�\��ˉ �unv�BZF��5��\6�M�,��t!3��i~r��g�t"t_��N;Bx۸��M�~�F��ѹ�00pm0,��(�T��_&��"˸*�(U��e��c�^��9-��2�*��" 1[g�^n��Z��p��i��}��8 �,��o��-��!Uɥ��&F��=�"�aC!��tWۅ qJn��(�A)@Y��8g3�,�W��&Ǻ%!(7>�Ů�<�� є�0�蓢<��Dgs��=)Q��%PKu��#��D9��9)�e:��`�'HRA�R�ZJ��2ⵌ�[F\K� �B�_ĩ��/�8 ��w��I��C�^ɤ_�{��͸��_��߽�%x��d{�y2�$G��,�苅�y�ev/��Һ�'��cb|�{�@� !L�k�o:S��Ǘ�3m+�k��-��*G]B

matlab中想用filter函数实现IIR滤波器
matlab中想用filter函数实现IIR滤波器

matlab中想用filter函数实现IIR滤波器
说明：（1）为了使滤波器阶数尽可能低,每个滤波器的边界频率选择原则是尽量使滤波器过渡带宽尽可能宽.
（2）与信号产生函数mstg相同,采样频率Fs=10kHz.
（3）为了滤波器阶数最低,选用椭圆滤波器.（之后,再依次实现巴特沃斯、切比雪夫1、切比雪夫2数字滤波器）
按照图2 所示的程序框图编写的实验程序为exp1.m.
2、实验程序清单
%实验1程序exp1.m
% IIR数字滤波器设计及软件实现
clear all;close all
Fs=10000;T=1/Fs; %采样频率
%调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st
st=mstg;
%低通滤波器设计与实现
fp=890;fs=600;
wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指标（低通滤波器的通、阻带边界频）
[N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp
[B,A]=ellip(N,rp,rs,wp,'high'); %调用ellip计算椭圆带通DF系统函数系数向量B和A
y3t=filter(B,A,st); %滤波器软件实现
% 高低通滤波器设计与实现绘图部分
figure(4);
subplot(2,1,1);
myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线
yt='y_3(t)';
subplot(2,1,2);
tplot(y3t,T,yt); %调用绘图函数tplot绘制滤波器输出波形
function myplot(B,A)
%时域离散系统损耗函数绘图
%B为系统函数分子多项式系数向量
%A为系统函数分母多项式系数向量
[H,W]=freqz(B,A,1000);
m=abs(H);
plot(W/pi,20*log10(m/max(m)));grid on;
xlabel('\omega/\pi');ylabel('幅度(dB)')
axis([0,1,-80,5]);title('损耗函数曲线');
function tplot(xn,T,yn)
%时域序列连续曲线绘图函数
% xn:信号数据序列,yn:绘图信号的纵坐标名称（字符串）
% T为采样间隔
n=0:length(xn)-1;t=n*T;
plot(t,xn);
xlabel('t/s');ylabel(yn);
axis([0,t(end),min(xn),1.2*max(xn)])
function st=mstg
N=2000
Fs=10000;T=1/Fs;Tp=N*T;
t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;
fc1=Fs/10;
fm1=fc1/10;
fc2=Fs/20;
fm2=fc2/10;
fc3=Fs/40;
fm3=fc3/10;
xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t);
xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t);
xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t);
st=xt1+xt2+xt3;
fxt=fft(st,N);
subplot(3,1,1)
plot(t,st);grid;xlabel('t/s');ylabel('s(t)');
axis([0,Tp/8,min(st),max(st)]);title('(a) s(t)的波形')
subplot(3,1,2)
stem(f,abs(fxt)/max(abs(fxt)),'.');grid;title('(b) s(t)的频谱')
axis([0,Fs/5,0,1.2]);
xlabel('f/Hz');ylabel('幅度');

matlab中想用filter函数实现IIR滤波器 MATLAB中FILTER函数有关的问题（急）题目是这样的：考虑函数conv和filter之间的关系,试利用filter函数来实现离散时间信号的卷积. 在matlab里面如何实现对一串数据的移动平均数?用filtfilt,filter还是其他函数? 关于matlab filter这个函数^^y = filter(b,a,X)这个函数期中b和a代表什么呢?看过了注释,如果要实现fir 高通,如何设置,为什么呢? 关于matlab filter这个函数到底怎么用? matlab中filter和conv函数有什么区别? 用MATLAB编写递归函数实现计算n2时,f(i)=f(i-1)+f(i-2). matlab函数表示用matlab怎么实现啊谢谢 matlab中什么函数可以实现离散函数变为连续函数图像低通滤波作用请问图像为什么要低通滤波,低通滤波的作用是什么?matlab语句如下：% low-pass filter the imagefiltered = conv2fft(dvalue,filtro,'same');conv2fft是一个自编的函数.具体怎么实现不用去管,就 filter 求matlab中实现离散序列傅里叶变换的函数. matlab用一个数乘一个函数怎么实现 matlab 中hilbert函数具体实现的什么运算? 求gabor 小波变换实现函数,matlab 如何用matlab实现两个函数的卷积运算 matlab中如何实现Q函数的反函数 Matlab的函数中如何实现按引用传递参数?