多媒体应用基础_第四节_MATLAB.pptVIP

* * * * * * * * * * 信息管理与信息系统 延边大学 MATLAB基础 主讲：崔基哲 延边大学经济管理学院 办公室：人文楼306 多媒体应用基础 离散信号 周期序列 ex4 产生一个幅度为2，频率为4Hz,相位为 的正弦信号 A=2; phi=pi/6; omega=2*pi/12; n=-10:10; x=A*sin(omega*n+phi); stem(n,x,fill); grid on; 单位脉冲序列 和单位阶跃序列u(n) 可以用ones(1,n)和zeros(1,n)来生成单位脉冲序列和单位阶跃序列 ones(1,n)命令产生1行n列的1值 zeros(1,n)命令产生1行n列的0值 ex5 产生单位脉冲序列 和单位阶跃序列u(n-3) n=-2:30; x=[zeros(1,5),1,zeros(1,27)]; y=[zeros(1,5),ones(1,28)]; subplot(2,1,1); stem(n,x,fill);grid on; subplot(2,1,2) stem(n,y,fill);grid on; 差分与求和、微分与积分 例:已知三角波f(t)，画出其微分与积分的波形 %differentiation %integration h=0.001;t= -3:h:3; t= -3:0.1:3; y1=diff(f2_2(t))*1/h; for x=1:length(t) plot(t(1:length(t)-1),y1) y2(x)=quad(f2_2, -3,t(x)); end plot(t,y2) 差分 y=diff(f); 求和 y=sum(f(k1:k2)); 微分 y=diff(f)/h; h为数值计算所取时间间隔 定积分 quad(‘function_name’,a,b); %function_name为被积函数名，a和b指定积分区间。 其中函数 f2_2(t)为： function ft2=f2_2(t) ft2=tripuls(t,4,0.5); 已知一个系统的传递函数为 求系统单位冲激响应h[n]，以及h[n]的幅频相频响应图 N=64; x=[1 zeros(1,N-1)]; %产生单位冲激函数 b=[0.008 -0.033 0.05 -0.033 0.008]; a=[1 2.37 2.7 1.6 0.41]; y=filter(b,a,x); %计算单位冲激响应 figure(1); n=1:N; stem(n,y,fill);grid on; title(单位冲激响应); figure(2) Fs=1024; freqz(b,a,N,Fs);grid on; %做出幅频和相频响应曲线 离散时间系统单位脉冲响应的求解 h=impz(b,a,k) b, a 分别是差分方程左、右端的系数向量 ? k 表示输出序列的取值范围  h 就是单位脉冲响应 求系统y[k]+3y[k-1]+2y[k-2]=10f [k]的单位脉冲响应 k=0:10; a=[1 3 2]; b=10; h=impz(b,a,k); stem(k,h); FIR数字滤波器的设计 IIR数字滤波器最大缺点：不易做成线性相位 现代图像、语音、数据通信对线性相位的要求 线性相位的FIR得到大力发展和应用 线性相位FIR数字滤波器的特点 系统函数为： 其频响为： 线性相位的FIR滤波器是指其相位函数 满足线性方程： α式中表示系统群时延，β表示附加相移 当h(n)为实数且偶对称时，FIR滤波器为 恒相时延 相位曲线为一条过原点、以-(N-1)/2为斜率的直线 当h(n)为实数且奇对称时，FIR滤波器为 恒群时延 相位曲线为一条截距为π、以-(N-1)/2为斜率的直线 由h(n)的奇偶二种对称性以及N等于奇偶不同情况可以有四种频率特性与之对应。 阶数N的一般选择方法 ω π Δω ω π ωp ωc 0 1 ：通带纹波