信息工程概论实验.docVIP

《信息工程概论》实验 实验一 实验名称：傅里叶变换的实现 实验性质：验证性 实验目的：加深对短时傅里叶变换、离散傅里叶变换等的理解，并对其在系统信号中的应用有一定的了解。 实验要求：分组与独立相结合 实验内容： 傅里叶变换的实现。包括：短时傅里叶变换 离散傅里叶变换 连续傅里叶变换 信号分析中的傅里叶变换 建议语言：Matlab C++ Java 作业： 1、实现P81例3.1 例3.2 P99例3.8 例3.9 2、实验作业1中的逆变换 3、扩展WVD变换。P90例3.4 P91例3.6 P91例3.7验证两个信号和的WVD与各信号WVD和的关系。 4、Gabor展开 参考内容： 1利用matlab命令实现傅里叶变换 例：求 2利用matlab命令实现傅里叶逆变换 例：求的傅里叶逆变换f(t) 3.求DFT的幅值谱 例：求=5+6，以进行采样，求其DNF的幅值谱。 4. 求 f(t)=的短时傅立叶变换 解：取窗函数为 所以所求短时傅里叶变换为 5、WVD变换也可利用傅里叶变换来做 实验二 实验名称：小波变换的实现 实验目的：加深对小波变换的理解及其在信号处理中的应用 实验内容： 小波变换 信号分析中的小波变换 注：可以利用积分来实现 自由选定信号函数和母小波 有能力的同学练习多分辨分析 参考内容：已知Mexican帽小波 r= f（t）= 则f(t)的连续小波变换为： 实验三 实验名称：Koch（科克）曲线的生成 实验目的：理解科克曲线的生成原理 了解它的意义 实验内容： 科克曲线生成 有时间的同学练习sierpinski三角形 cantor三分集的生成 参考内容： 1、科克曲线： p=[8;0]; q=[6;0]; koch(p,q,4) 2、sierpinskin三角形 a=1;b=0.5;c=1; k=7; A=zeros(2,3^(k+1)); A(:,1:3)=[0 a b;0 0 c]; for n=1:k B=1/2*A; A(:,1:3^n)=B(:,1:3^n); A(:,3^n+1:2*3^n)=B(:,1:3^n)+1/2*[a;0]*ones(1,3^n); A(:,2*3^n+1:3^(n+1))=B(:,1:3^n)+1/2*[b;c]*ones(1,3^n); end for i=1:3^k patch(A(1,3*i-2:3*i),A(2,3*i-2:3*i),b); end 3、Cantor三分集 function cantor3(n) if nargin==0; n=6; end hold on [x,f]=ccp(0,0.5,n,0); axis off function [x,f]=ccp(x,f,n,y) plot([x-f,x+f],[y,y],r,linewidth,2); if n0 [x,f]=ccp(x-f/1.5,f/3,n-1,y-1); f=f*3; x=x+f/1.5; [x,f]=ccp(x+f/1.5,f/3,n-1,y-1); f=f*3; x=x-f/1.5; end 实验四 实验名称：信号的频谱图 实验目的：加深对不同种已调信号特点的理解 实验内容： 画出抑制载波双边带调幅信号在一个周期内的包络图 画出抑制载波双边带调幅信号的频谱 画出抑制载波单边带调幅各点的频谱 画出包络检波器输出的频谱 作业：P69 2.1 2.3 2.4 参考内容： 1、波形图： t=0:0.01:1; f1=20*cos(6.28*t)+50*cos(12.56*t); f2=cos(2*pi*10^6*t); f3=conj(f2); ； f4=f1.*f3; plot(t,f4) 频谱图： t=0:0.01:1; f1=20*cos(6.28*t)+50*cos(12.56*t); f2=cos(2*pi*10^6*t); f3=conj(f2); f4=f1.*f3; S=fft(f4) plot(t,S) 2、包络图：t=0:0.01:3; f1=20*cos(6.28*t)+50*cos(12.56*t); f2=cos(2*pi*10^6*t); f3=conj(f2); f4=f1.*f3; plot(t,f4) hold on ind=find(diff(sign(diff(f4)))==+2)+1; val=interp1(t(ind),f4(ind),t,lin