还剩11页未读，继续阅读
本资源只提供10页预览，全部文档请下载后查看！喜欢就下载吧，查找使用更方便
文本内容:
第章抽出的代码2MATLAB%ch2j
23.m%ch2j
23.mt=-10:
0.01:10;、“兀’中函数表示的是xt=sinc t;%MAT LABsine tntplot t,xtxlabel*t/s1ylabel1x t1抽样函数信号，title1%ch2j
24.m%ch2j
24.mt=-10:
0.01:10;为关系运算，成立时输出结果为不成立时输出结果为ut=t=0;%t=0t=0l0zplot t,ut;xlabel1t/s1ylabel1ut1单位阶跃信号，;title1%指定横、纵坐标的范围axis[-1010,-
0.
21.2]z z%ch2j
26.m%ch2j
26.mt=-10:
0.01:10;utl=-l.*t0;ut2=t=0;ut=utl+ut2;plot t,ut;xlabel1t/s1ylabel1ut1，符号函数信号，；title%cp21iti
4.m%cp21iti4,m t=-4:
0.001:4;xtl=
0.5*t+
2.*t0t=-4;xt2=t=0t
2.*2;xt3=-t+
4.*t=2t=4;%原信号xt=xtl+xt2+xt3;figure subplot221为绘制子图，表示在一张圆布上画行列共张%subplot subplot221224图，第位数表示第几张图plot t,xt grid3%绘制原信号on xlabel1t/s1;%绘制网格title1x t1;subplot222plot-t,xt%绘制原信号的反折信号%cp6jl
051.m改变抽样频率倍数；w0=2*pi*50;mi=[5,3,2,l];for%m=5,3,2,li=l:lengthmim=mi i;改变抽样频率；Ts=l./50*13*m;%m=5,3,2,1抽样点数N=2*pi/Ts/w0;%N=65,39,26,13;代码改变抽样频率%MATLAB%cp6j
1051.mn=0:1:N-1;・xnl=10*sinw0*n*Ts+
0.1*pi;xn2=5*sinw0*n*Ts*5+02*pi;xn3=3*sinw0*n*Ts*9+
0.3*pi;xn4=sinw0*n*Ts*13+
0.4*pi;%建立xn=xnl+xn2+xn3+xn4;x n1=i-1*2;时域波形；subplot4,3,I+i stem n,xn;%x nxlabel1n1;ylabel1x n1;title[!N=1num2strN];z%X=fftx,N;k=0:1:N-l;用按定义计算点Xk=xn*exp-j*2*pi/N*n**k;N DFT%取式f=k/Ts/N;6-123b%幅值频谱图形subplot4,3,I+i+1;stem f,abs Xk/N;xlabel1f/Hz1;ylabel1I X k|/N1;；1111title[fs=,num2strm,*650Hz]%幅值频谱局部放大图fsubplot4,3,1+i+2;stem f,abs Xk/N;title局部放大图|Xk|D;xlabel1f/Hz1;ylabel1I Xk|/N1;axis[0,800,0,6];end%%%%%%%%%%%%%%%%%%%%%%%%%%%%cp6jl
052.m代码改变抽样时长%MATLAB%cp6j
1052.m w0=2*pi*50;m=5;需改变抽样时长Ni=[65,3317];rfor i=l:lengthNi Ts=l./50*13*m;N=Nii;n=0:1:N-1;xnl=10*sinw0*n*Ts+
0.1*pi;xn2=5*sinw0*n*Ts*5+
0.2*pi;xn3=3*sinw0*n*Ts*9+
0.3*pi;xn4=sinw0*n*Ts*13+
0.4*pi;%建立xn=xnl+xn2+xn3+xn4;x n1=i-1*1;工+时域波形subplot3,2,i;stem n*Ts*1000,xn;%x nxlabel1nT/ms;ylabelT x n f;%title[1N=*,num2str N];axis[0,20,-20,20];%按£计算点%X=fft xN;1N DFTzk=0:1:N-l;%按定义计算点Xk=xn*exp-j*2*pi/N*n**k;N DFT招取式f=k/Ts/N;6-123b%幅值频谱图形subplot3,2,I+i+l;stem f,abs Xk/N;1111xlabelf/Hz;ylabelI Xk|/N;%r=2A i-1;%title[*T=T5o/,num2str r];oaxis[01/Ts,0,10];end%%%%%%%%%%%%%%%%%%%%%%%%%%%z%cp6jl
053.m代码抽样信号补%MATLAB0%cp6j
1053.m w0=2*pi*50;Ts=l./50*13*5;豺由样点数N=2*pi/Ts/w0;n=0:1:N;xnl=10*sinw0*n*Ts+
0.l*pi;xn2=5*sinw0*n*Ts*5+
0.2*pi;xn3=3*sinw0*n*Ts*9+
0.3*pi;xn4=sin w0*n*Ts*13-i-
0.4*pi;xn=xnl+xn2+xn3+xn4;%时域补x2n=[xn1:1:33zeros1,32];0务时域波形1subplot1,2,1;stem n*Ts*1000,x2n;xlabelnT/ms-20T;ylabelf xn f;axis[0,20,f20];%X=fftx2,N;k=0:1:N;1Xk=x2n*exp-j*2*pi/N*n*k;%DFTf=k/Ts/N;%频域波形subplot1,2,2;stem f,abs Xk/N;f1111xlabel/Hz;ylabel IXkl/N;axis[0,l/Ts0,5];z%cp6jl
054.m代码非整周期抽样%MATLAB w0=2*pi*50;Ts=l./50*13*5;抽样点数Ni=[65,80,95];%N=65,80,95;for i=l:lengthNi;N=Ni i;n=0:1:N-1;xnl=10*sinw0*n*Ts+
0.1*pi;x2=5*sinw0*n*Ts*5+
0.2*pi;xn3=3*sinw0*n*Ts*9+
0.3*pi;xn4=sinw0*n*Ts*13+
0.4*pi;xn=xnl+x2+xn3+xn4;1=i-1*1;超时域波形subplot3,2,I+i;stem n*Ts*1000,xn;1111xlabelnT/ms;ylabelx n;11title[N=,num2str N];axis[0,30,-20,20];%X=fft x,N;k=0:1:N-l;X=xn*exp-j*2*pi/N*n**k;%DFTf=k/Ts/N;%频域波形；subplot3,2,1+i+l stemf,abs X/N;1f11xlabelf/Hz;ylabel|Xk|/N;axis[0,1/Ts,0,6];end%绘制反折信grid onxlabel*t/s1;title1x-t,;subplot223Plot-t4-4xtz号的时移信号，波形沿时间轴向右平移个单位4%grid onxlabel*t/s1;title1x-t+41;subplot224plot
0.5*-t+4,xt信号沿时间轴压缩到原来的1/2grid onxlabel*t/s1;title*x-2t+4*;%cp21iti
215.m%生成信号%cp21iti
215.m dt=
0.01;kl=0:dt:1;xi t%生成信号xl=ones1,lengthkl;k2=0:dt:2;X2tx2=onesl lengthk2;z%计算卷积函数为计算卷积的函数，要注y=dt*conv xl,x2;xl t*x2t,conv意用求连续信号离散序列的卷积的不同convk0=kl1+k21;k3=lengthkl+lengthk2-2;故确定卷积的非零样值的时间向量k=k0:dt:k0+k3*dt;xi t*x tplotk,y xlabel
2、卷积结果，t/sz;title第章抽出的代码3MATLAB例3-7%cp31iti
7.mA=1;T=2;tao=l;t=-2:
0.001:2;N=input1Number ofharmonic=*X0=A*tao/T;wO=2*pi/T;X=X0*ones1,length t;for k=l:1:NX=X+2*X0*sinck*w0*tao/2/pi*cosk*w0*t;endplott,X例3-9%cp31iti
9.msyms x t衰减正弦信号x=exp-2*t*sin2*pi*t*heavisidet subplot311ezplotx title11X=fourier x;「幅度谱X=simplifyX subplot312ezplotabsX title D subplot313ezplotangleX「相位谱titleD例3-20%cp31iti
20.m%cp31iti
20.mw=0:
0.01:5;b=[l];a=[l221];h=freqsb,a,w;%subplot2,1,1;plot wabs h;zxlabel1\omegarad/s;ylabel1IHj\omega|*grid on;%subplot2,1,2;plotw,angle h;xlabel\omegarad/s1;ylabel\phi\omega|grid on;第章抽出的代码4MATLAB实现的代码见MATLAB cp4liti
5.m,%cp41iti
5.mxtl=exp-2*t;xt2=exp-5*t;xt=2*xtl+5*xt2;Xs=laplace xt;运行结果为Xs=2/s+2+5/s+5实现的代码见MATLAB cp4liti
6.m,%cp41iti
6.msyms tk Tpositive xt=diract-k*T;Xs=laplace xt运行结果为Xs=exp-T*k*s%cp41iti
20.m%分子系数，按降基排列num=
[1];denl=conv
[15],
[15];den2=conv
[15],
[13];%分母系数,按降幕排列den=conv deni,den2[r,p,k]=residuenum,den运行结果为r=-
0.1250-
0.2500-
0.
50000.1250P=-
5.0000-
5.0000-
5.0000-
3.0000k=[]%cp41iti20a.mnum=
[24];den=[l040];[r,p,k]residue num,den运行结果为-
0.5000-
0.5000i-
0.5000+
0.5000i
1.0000P=0+
2.0000i0-
2.0000i0k=[]%cp41iti22b.msyms s;Xs=2*s+
4./sA3+4*s国拉普拉斯反变换xt=ilaplace Xs运行结果为xt=sin2*t-cos2*t+1%cp41iti23a.msyms sX=s-4/s/s+1A3x=ilaplaceX%cp41iti23b.mnum=[1-4];a=conv
[10]
[11];fb=conv
[11]
[11];den=conva,b;z[r,p,k]=residuenum,den运行结果为r=
4.
00004.
00005.0000-
4.0000P=-
1.0000-
1.0000-
1.00000k=[]%cp41iti
26.msyms st;*拉普拉斯变换Xs=laplace exp-2*t%拉普拉斯变换Hs=laplace exp-5*tYs=Xs*Hs yt=ilaplaceYs*拉普拉斯反变换运行结果为Xs=1/s+2Hs=1/s+5Ys=1/s+2*s+5yt=exp-2*t/3-exp-5*t/3%cp41iti
29.m%分子多项式系num=
[26]den=[l数%分母多项式系32]z=rootsnum数%求零点%求极p=rootsden点运行结果为z=-3P=-2-1%cp41iti
34.mk=3;w=-k*pi:
0.01:k*pi;num=
[14];den=
[134];H=freqrespnum,den,sqrt-1*w;subplot2,1,1;plotw/pi,absH;xlabel1\omega/\pi1;ylabel1IHj\omega|1;一幅频特性曲线，;titlegrid on;subplot2,1,2;plotw/pi,angleH/pi;xlabel1\omega/\pi1;ylabel1\phi\omega/\pi*;「相频特性曲线，;titlegrid on;%cp41iti35b.mnum6=
[1];%极点p=0;阶极点den61=conv[1-p]/
[1]%1阶极点den62=conv[1-p],den61%2阶极点-p],den62den63=conv[1%3-p],den63阶极点den64=conv[1%4%画阶极点的〃/的波形impulsenum6,den63;3%cp41iti36a.mb=[l4];a=[l320][r pk]=residueb,az z运行结果为r=[1-32]p=[-2-10]k=[]%cp41iti36b.m%画单位冲激响应的时域波形第分子多项式系数num=
[14];%分母多项式系数den=[l320];sys=tfnum,den impulsesys;%cp41iti36c.m%画系统的零、极点图num=
[14];den=[l320];sys=tfnum,den;pzmapsys;第章抽出的代码5MATLAB%cp5j
21.m%cp5j21,mn=-5:1:5;x=0*n=-5n0+1*n==0+0*n0n=5;为绘制火柴棒图，这里主要用来绘制离散时间图形stem n,x Jfilled1;%stem axis[-6,6,-
0.2,
1.2]%cp5j
22.m%cp5j
22.mn=-5:1:5;x=0*n0+1*n=0;stemn,x,1filled1;axis[-66-
0.2,
1.2]z z第章抽出的代码6MATLAB%cp6j
23.m代码周期矩形序列的%MATLAB DFS%周期N=20;%半脉宽Nl=2;M=N1*2+1;n=-2*N+M+l/2:1:2*N-M+l/2;fO=zeros;fl=ones1,M;x=[fO fl,fO,fl,fO,fl,fO];%产生zXN nsubplot2,1,1;stemn,x;%绘制图XNFIxlabel1n1;title1xNn T;n=-2*N+M+l/2:1:2*N-M+l/2;k=-2*N+M+l/2:1:2*N-M+l/2;WN=exp-j*2*pi/N;nk=n1*k;Xk=x*WN.Ank/N;%依定义式计算系数DFS Xksubplot2,1,2;stemk,Xk,b1;N%绘制图X kxlabel1k1;title1XNk1;Nhold on;plotk,Xk,1r1;hold off;%绘制包络图XNQ%cp61itil
0.m代码的频率分辨率%MATLAB DFTNl=50;fs=200;Ts=l/fs;nl=0:1:N1-1;xl=cos49*2*pi*nl*Ts+cos51*2*pi*nl*Ts;%用计算Xl=fft xlNl/Nl;fftDFTf兀kl=0:l:Nl-l;wl=2*pi/Nl*kl;%Q=k*2/Nsubplot331;stem%iSi xl nxlabel1n1;title1xln1;axis[0,200,-
2.5,
2.5];subplot332;五为单位画stemwl/pi*100,abs Xl,
1.1;%f=Q/*100,Hz|X1k|xlabel1f/Hz^/titleC1|Xlk||;axis[0200,0,
0.55];zsubplot333;%画|局部放大stemwl/pi*100,abs Xl,
1.f;|X1k「|局部放大1xlabel ff/Hz，|X1k;axis[43,57,0,
0.55];%%%%%%%%%%%%%%%%%%%N2=200;n2=0:l:N2-l;x2=[xl1:1:50zeros1,150];X2=fftx2,N2/N2;k2=0:l:N2-l;兀w2=2*pi/N2*k2;%Q=k*2/N；subplot334;stemn2,x2,,xlabel1n1;title1x2n1;axis[0,200,-
2.5,
2.5];兀；subplot335;stemw2/pi*100abs X2,.%f=Q/*100HzAxlabel1f/Hz,;title,|X2k|f;axis[0,200,0,
0.55];；subplot336;stemw2/pi*100absX2,.z|局部放大，xlabel ff/Hz，|X2k axis[43,57,0,O.bbJ;%%%%%%%%%%%%%%%%%%%N3=200;n3=0:l:N3-l;x3=cos49*2*pi*n3/200+cos51*2*pi*n3/200;X3=fftx3N3/N3;rk3=0:l:N3-l;兀w3=2*pi/N3*k3;%Q=k*2/N・subplot337;stemn3,x3,;xlabel1n1;title1x3n1;axis[0,200,-
2.5,
2.5];subplot338;stemw3/pi*100,absX3,
1.1;%f=Q/n*100Hzxlabelff/Hz*/title1|X3k|1;axis[0,200,0,
0.55];subplot339;stemw3/pi*100absX3,.;z|局部放大xlabel1f/|X3k D；axis[43,57,0,
0.55];%%%%%%%%%%%%%%%%%%%%cp6jl
01.m代码信号的表示%MATLAB%基波频率w0=2*pi*50;w0t=0:
0.00001:
0.04;xtl=10*sinw0*t+
0.1*pi;xt2=5*sinw0*t*5+
0.2*pi;xt3=3*sinw0*t*9+
0.3*pi;xt4=sinw0*t*13+
0.4*pi;%设置xt=xtl+xt2+xt3+xt4;x t%画乂七波形plott,xt;xlabel1t/s1;title1xt1;grid on;%cp6jl
02.m代码信号的抽样%MATLAB%基波频率w0=2*pi*50;w0%抽样周期Ts=l./50*13*5;Ts%抽样点数N=2*pi/Ts/w0;Nn=0:1:N-l;xnl=10*sinw0*n*Ts+
0.l*pi;xn2=5*sinw0*n*Ts*5+
0.2*pi;xn3=3*sinw0*n*Ts*9+
0.3*pi;xn4=sinw0*n*Ts*13+
0.4*pi;xn=xnl4-xn2+xn3+xn4;%建立x n%画乂口波形stemn,xn;xlabel1n1;title1xn1;grid on;%cp6jl030m代码变换%MAT LABDFT%用函数求点变换%Xk=fft xn,N;fftN DFTk=0:1:N-1;%按定义求点变换Xk=xn*exp-j*2*pi/N*nf*k;N DFT%画幅频特性图stem k,abs Xk;xlabel1k1;ylabel1I Xk|1;%cp6jl
04.m代码检查频谱的幅值及相位%MATLAB DFTk=[l5,9,13];zfor i=l:lengthkki=kiAk=absXkki+1*2/N pk=angleXkki+1/pi end。
个人认证
优秀文档
获得点赞 0
