matlab在信号处理中的应用课件.pptVIP

第4章MATLAB在信号处理中的应用;4.1信号及其表示;t=0:0.1:30; %对时间变量赋值

x=exp(-0.707*t).*sin(2/3.*t);%计算变量所对应的函数值

plot(t,x);%绘制函数曲线

grid;%图形窗口加网格

xlabel(time(sec));%标注x轴,y轴

ylabel(x(t));

axis([-0.05,30,-0.05,0.35]);%确定x轴，y轴显示范围;函数名;t=0:0.0001:1;

x1=sawtooth(2*pi*50*t,0.8)

plot(t,x1)

axis([0,0.2,-1,1])

;t=0:0.001:2

y=chirp(t,0,1,150)

figure(1)

plot(t,y)

axis([0,0.5,0,1]);t=0:0.01:1;

y1=tripuls(t);

y2=tripuls(t,0.6);

subplot(211);plot(t,y1);

subplot(212);plot(t,y2);

;4.1.3离散时间信号的表示;4.1.4几种常用离散时间信号的表示;3．实指数序列;functionfsxl(delta,omig,n1,n2)

e=delta+omig.*j;

n=[n1:n2];

x=exp(e*n);

x_real=real(x);%生成实部序列

x_imag=imag(x);%生成虚部序列

x_magnitude=abs(x);%生成幅度序列

x_phase=(180/pi)*angle(x);%生成相位序列;4.2信号的基本运算;n1=[-5:4];%序列x1(n)的起始及终止位置

n1s=-5;n1f=4;

x1=[231-13421-5-3];%序列x1(n)不同时间的幅度

n2=[0:9];%序列x2(n)的起始及终止位置

n2s=0;n2f=9;

x2=[1111111111];%序列x2(n)不同时间的幅度

ns=min(n1s,n2s);nf=max(n1f,n2f);%求新信号的时间起始终止位置

n=ns:nf;

y1=zeros(1,length(n));%延拓序列初始化

y2=zeros(1,length(n));

y1(find((n=n1s)(n=n1f)==1))=x1;%给延拓序列y1赋值x1

y2(find((n=n2s)(n=n2f)==1))=x2;%给延拓序列y2赋值x2

ya=y1+y2;%逐点相加

yp=y1.*y2;%逐点相乘

subplot(411);stem(n,y1,.);

subplot(412);stem(n,y2,.);

subplot(413);stem(n,ya,.);

subplot(414);stem(n,yp,.);;;;;4.2.4两序列的卷积运算;例如：求序列;4.3信号的能量和功率;非周期三角波信号能量的MATLAB计算：;4.4线性时不变系统;4.4线性时不变系统;注意：在MATLAB中，传递函数由分子、分母两个多项式的系数来表示，系数为降幂排列。例如：;4．极点留数模型;系统模型的MATLAB表示;4.4.2系统模型的转换函数;线性系统模型的变换函数;在???令窗口输入：

z=[3];p=[1,2];k=2;

[A,B,C,D]=zp2ss(z,p,k)

屏幕显示结果：

A=

3.0000-1.4142

1.41420

B=

1

0

C=

2.0000-4.2426

D=

0;[例4-19]求离散时间系统;4.4.3系统互联与系统结构;MATLAB实现函数parallel()

格式：[A,B,C,D]=parallel(A1,B1,C1,D1,A2,B2,C2,D2)或

[num,den]=parall