1、t=-pi:pi;f=2*cos(3*t+pi/6);f(t)=2*cos(3*t+pi/6)axis(-pi,pi,-3,3)(3)f=2*sin(3*t+pi/6);f(t)=2*sin(3*t+pi/6)2、用MATLAB表示抽样信号(sinc(t)、矩形脉冲信号(rectpuls(t, width)及三角脉冲信号(tripuls(t, width, skew)。t=-3:3;f=sinc(t);f(t)=sinc(t)axis(-3,3,-0.5,1.1)f=rectpuls(t,4);f(t)=rectpuls(t,4)axis(-3,3,0,1.5)t=-2:2;f=tripuls
2、(t,1,0);f(t)=tripuls(t,1,0)axis(-2,2,0,1.1)3、编写如图3的函数并用MATLAB绘出满足下面要求的图形。图3源程序:clcclearcloset0=0:0.05:12;t1=-12:14;t2=0:t3=-6:1;t4=-3:23;f0=4*rectpuls(t0-6,12)+3*tripuls(t0-6,4,0);f1=4*rectpuls(-t1-6,12)+3*tripuls(-t1-6,4,0);f2=4*rectpuls(t2-8,12)+3*tripuls(t2-8,4,0);f3=4*rectpuls(1-2*t3-6,12)+3*tri
3、puls(1-2*t3-6,4,0);f4=4*rectpuls(0.5*t4+1-6,12)+3*tripuls(0.5*t4+1-6,4,0);subplot(3,2,1);plot(t0,f0);f(t)subplot(3,2,3);plot(t1,f1);f(-t)subplot(3,2,4);plot(t2,f2);f(t-2)subplot(3,2,5);plot(t3,f3);f(1-2t)subplot(3,2,6);plot(t4,f4);f(0.5t+1)五、实验总结实验二 离散时间信号时域分析1. 掌握常用离散信号的MATLAB表示方法。2、掌握用MATLAB计算卷积和的
4、方法。二、实验原理分析三、实验内容1、用MATLAB表示离散信号:,。k=-3:fk=2.(k);subplot(1,2,1);stem(k,fk,filledaxis(-4,4,0,9);h=-7:7;ak*uk)fh=2*sin(h);subplot(1,2,2);stem(h,fh,axis(-7.5,7.5,-2.5,2.5);A*sink2、已知离散序列波形如图4所示,试用MATLAB绘出满足下列要求的序列波形。(1) fk-2uk (2) f-k (3) f-k+2 (4) fk-2uk-2图4k=-6:5;f=0,1,3,6,10,15,14,12,9,5,0,0;stem(k,
5、f,fkh=0:f1=10,15,14,12,9,5,0,0;stem(h,f1,fk-2*uksubplot(3,2,4)stem(-k,f,f-ksubplot(3,2,5)stem(-k+2,f,f-k+2h=2:f1=14,12,9,5,0,0;fk-2*uk-23、若,计算x=1 1 1 1 0 0;y=conv(x,x);y执行结果:y =1 2 3 4 3 2 1 0 0 0 0实验三 连续时间系统时域分析1、掌握卷积计算方法。2、掌握函数lsim,impulse,step的用法,lsim为求取零状态响应,impulse为求取单位脉冲响应,step为求取单位阶跃响应。3、运用课堂
6、上学到的理论知识,从RC、RL一阶电路的响应中正确区分零输入响应、零状态响应、自由响应与受迫响应。安装有matlab6.5以上版本的PC机一台。1. 分别用函数lsim 和卷积积分两种方法求如图7 所示系统的零状态响应。其中L=1,R=2,e(t) = (t),i(0 ) = 2。方法(1):函数lsim法clear;ts=0;te=10;dt=0.01;t=ts:dt:te;f=exp(-1*t);a=1 2;b=1;y=lsim(b,a,f,t);plot(t,y);Time(sec)ylabel(i(t)方法(二):卷积积分法syms t xe=exp(-1.*(x);h=exp(-2.
7、*(t-x);i=int(e.*h,x,0,t);ezplot(i,0 8);exp(-1t)*(-2t)2. 求上述系统的冲激响应与阶跃响应。a:冲激响应求解impulse(b,a);axis(0,8,-0.2,1.1);b:阶跃响应求解step(b,a);axis(0,5,0,0.6);五、思考题1. 为什么连续时间系统的零状态响应为激励与冲击响应的卷积?答:一方面,零状态响应与系统的特性和外加激励有关;另一方面,根据卷积的定义,函数e(t)与函数h(t)相卷积,即对于某一t值时乘积e()h(t-)曲线下的面积,因此,连续时间系统的零状态响应为激励与冲击响应的卷积。2. 利用卷积积分法计算
8、系统响应应从几个方面进行?可以先求零输入响应,再求零状态响应,而零状态响应的求解需要将激励与冲击响应求卷积,最后,将零输入响应与零状态响应相加即得系统响应。六、实验总结实验四 离散时间系统时域分析掌握用MATLAB计算离散时间系统响应的方法。1、假设某系统的单位函数响应,系统激励信号,求系统的零状态响应,并画图表示。k=0:h=(0.8).k;e=1 1 1 1;y=conv(h,e);t=0:6;stem(t,y,零状态响应2、采用P61 recur.M 计算系统响应:画出时的系统响应。a=-0.5 -0.1;b=0 1 0;y0=1 0;x0=0 0;n=0:10;x=ones(1,len
9、gth(n);y=recur(a,b,n,x,x0,y0);stem(n,y,nyn系统响应recur.mfunctiony = recur(a,b,n,x,x0,y0)N=length(a);M=length(b)-1;y=y0 zeros(1,length(n);x=x0 x;a1=a(length(a):-1:1);b1=b(length(b):for i=N+1:N+length(n),y(i)=-a1*y(i-N:i-1)+b1*x(i-N:i-N+M);endy=y(N+1:N+length(n);四、实验总结实验五 连续信号的频域分析1. 掌握周期信号的频谱 Fourier 级数
10、的分析方法。2. 深入理解信号频谱的概念,掌握典型信号的频谱以及Fourier 变换的主要性质。3. 掌握调制与解调的基本原理及滤波器的使用。三、实验原理1. 求如图所示周期矩形脉冲信号的Fourier级数表达式,画出频谱图,并用前N次谐波合成的信号近似。syms t n;T=2*pi; %设T为2*piN=9; %设N为9f=heaviside(t)-2*heaviside(t-T/2)+heaviside(t-T); subplot(2,2,1);ezplot(f,0,2.1*pi);原函数h=exp(-j*n*2*pi/T*t); %用指数傅里叶级数表示,角频率2*pi/TA1=int(
11、f.*h,t,0,T);A=2/T*A1;for n1=-N:-1 C(n1+10)=subs(A,n,n1);for n1=1:NC(10)=0;%N=0时系数subplot(2,2,3);k=-N:N;stem(k,abs(C); %作出幅度谱Cn 的幅度Omegasubplot(2,2,4);stem(k,angle(C); %作出相位谱Cn的相位f1=0;for m=-N: f1=f1+1/2*C(m+10)*exp(j*m*t); %前N次谐波合成的信号近似subplot(2,2,2);ezplot(f1,0,2.01*pi);前N次谐波合成的信号近似2、试用fourier()函数求
12、下列信号的傅里叶变换F( j) ,并画出F( j)(1)f (t) = te3t (t)syms t x ;x=fourier(t*exp(-3*t)*heaviside(t);xz=abs(x);ezplot(z);%符号函数作图函数|F(j)|结果:x =1/(3+i*w)2(2)f (t) = sgn(t)x=fourier(2*heaviside(t)-1);%2*heaviside(t)-1即为sgn(t)x =-2*i/w3、调制信号为一取样信号,利用MATLAB 分析幅度调制(AM)产生的信号频谱,比较信号调制前后的频谱并解调已调信号。设载波信号的频率为100Hz。clc;clo
13、se;Fm=10;t1=0:0.00002:0.2;syms t v;x=sin(2.0*pi*Fm*t)/(2.0*pi*Fm*t);ezplot(x,0,0.2);Fx=fourier(x,v);subplot(3,2,2);ezplot(Fx,-50*pi,50*pi);axis(-50*pi,50*pi,-0.05,0.1);频谱y=x*cos(200*pi*t);b=subs(y,t,t1);plot(t1,b);调制后axis(0,0.2,-1,1);Fy=fourier(y,v);ezplot(Fy,-250*pi,250*pi);axis(-250*pi,250*pi,-0.0
14、5,0.1);z=y*cos(200*pi*t);Fz=fourier(z,v);G=-heaviside(v-20*pi)+heaviside(v+20*pi);%门函数Fx1=Fz*G;x1=2*ifourier(Fx1,v);%滤波过程中幅度减半且反向。ezplot(x1,0,0.2);解调后ezplot(2*Fx1,-50*pi,50*pi);1、根据试验1 的结果,解释Gibbs 现象。因为对于具有不连续点的函数,即使级数的项无限增大,在不连续处,级数之和不收敛于函数f(t);在跃变点附近的波形,总是不可避免的存在有起伏震荡,从而使跃变点的值超过一形成过冲,造成吉布斯现象。2、比较周
15、期信号与非周期信号的频谱。区别:1周期信号的频谱为离散频谱,非周期信号的频谱为连续频谱。2.周期信号的频谱为Fn的分布,表示每个谐波分量的复振幅;而周期信号的频谱为F(j)的分布,(F(j)/2表示合成谐波分量的复振幅,所以也将称为频谱密度函数。联系:1.都是反映将时域信号表示为正弦类信号时各谐波分量的分布特性。2.若周期信号是连续非周期信号的周期延拓,则两者的关系为F(j)= =3、调制与解调的基本原理是什么?为什么要进行调制?调制:调制就是使信号f(t)控制载波的某一个或某些参数(如振幅、频率、相位等),是这些参数按照信号f(t)的规律变化的过程。载波可以是正弦波或脉冲序列。以正弦型信号作
16、载波的调制叫做连续波调制。调制后的载波就载有调制信号所包含的信息,称为已调波。对于连续波调制,已调信号可以表示为:,它有振幅频率和相位三个参数构成。改变三个参数中的任何一个都可以携带同样的信息。因此连续波的调制可分为调幅、调相、和调频。解调:解调是调制的逆过程,它的作用是从已调波信号中取出原来的调制信号。调制过程是一个频谱搬移的过程,它将低频信号的频谱搬移到载频位置。如果要接收端回复信号,就要从已调信号的频谱中,将位于载频的信号频谱再搬回来。之所以进行解调,是因为无线电通信是通过空间辐射方式传输信号的,调制过程可以将信号的频谱搬移到容易以电磁波形势辐射的较高频率范围;此外,调制过程可以将不同的
17、信号通过频谱搬移托付至不同频率的载波上,实现多路复用,不至于互相干扰。实验六 连续系统的频域分析一、实验目的:掌握连续时间系统变换域分析的基本方法。四、实验内容:如图所示系统:(a)对不同的RC值,用freqs函数画出系统的幅频曲线。H(j)=1/(1+jRC),RC的取值依次为100、10、0.00001时的幅频曲线.b=0 1;for c=-5:2 RC=10c; a=RC 1; freqs(b,a); axis(10(-2),10(5),0.1,1); hold on得到一系列幅频曲线,从左到右依次为RC的取值依次为100、10、0.00001时的幅频曲线。(b) 信号包含了一个低频分量
18、和一个高频分量。确定适当的RC值,滤除信号中的高频分量并画出信号和在s范围内的波形。提示: |H( j)|为最大值的/2处对应的频率为通带截止频率c,首先求取|H( j)|并找到c和RC关系,然后根据题意选定c即可确定RC值。由(a)中的图可知,当RC=-2时符合题意。0.001:f=cos(100*t)+cos(2000*t);subplot(2,1,1);y1=cos(100*t)/(1+j*100*10(-2)+cos(2000*t)/(1+j*2000*10(-2);subplot(2,1,2);plot(t,y1)2、信号任选,分析以下几种情况下信号的频谱和波形变化:(1) 系统满足
19、线性不失真条件时;(2) 系统只满足恒定幅值条件时;(3) 系统只满足相位条件时;(4)系统两个条件均不满足时。利用fourier求取信号的傅立叶变换E(j),然后设计H( j) = H( j) e ( )使之满足不同条件,计算R( j)= E(j) H(j)并画频谱图。e=exp(-2*abs(t);subplot(2,3,1);ezplot(e,-3,3);axis(-3,3,-0.2,2);Fe=fourier(e,v);subplot(2,3,2);ezplot(Fe,-3,3);幅度谱axis(-3,3,0,2);i=1;for a=-3:0.02:3 R11=subs(Fe,v,a
20、); C(i)=angle(R11); i=i+1;b=-3:subplot(2,3,3);plot(b,C);相位谱axis(-3,3,-1,1);H1=2*exp(-j*v*1);R1=Fe*H1;r1=ifourier(R1,t);subplot(2,3,4);ezplot(r1,-3,3);满足线性不失真条件axis(-3,3,-0.2,2)subplot(2,3,5);ezplot(abs(8/(4+v2)*exp(-i*v)axis(-3,3,0,2.2); R11=subs(R1,v,a);subplot(2,3,6);axis(-3,3,-3,3);H1=(1-j*v)/(1+j*v);D=abs(R1
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1