西北工业大学Matlab实验报告

1、西北工业大学基于MATLAB的数字信号处理实验报告学院：学号：姓名：专业：实验时间：2016年6月指导教师：冯燕八、西北工业大学2016年06月实验一 MATLAB基本编程实验一、实验目的及要求1 .熟悉MATLA运行环境；2. 掌握MATLAB勺基本语法和函数；3. 掌握MATLAB勺基本绘图功能二、实验设备(环境)及要求1. 计算机2. MatIab软件编程实验平台三、实验内容与步骤1. 求下列线性方程组的解6x1 3x2 4x3 = 3-2x1 52 73 - -48x< - 4x? - 3x3 -72 .编写一个MATLA函数文件，该函数的功能是求出最大的 n值， 使得n! ：

2、real max0通过M文件或命令行调用该函数文件计算 realmax=50000 时的 n 值。3.某多项式的定义为：y = cos(m*cos'(x)其中X的值介于-1,1之间 当m的值由1变化到5,可以得到5条曲线。将这5条曲线画在同一图形 窗口上，并标明每一条曲线。四、设计思想1. 采用矩阵求逆的方法，将系数矩阵求逆，乘以单元矩阵，直接得到未知数。2. 在M文件中使用循环语句，判断是否符合条件，返回值，主程序进行调用3. 将5条曲线分别绘制，plot用不同颜色标志，采用hold On画在一个图中，五、程序代码及注释1. A=6,3,4;-2,5,7;8,-4,-3;%b=3;-

3、4;-7;x=i nv(A)*b%2. M文件：fun Cti On n = fun( max )%UNTITLED SUmmary of this function goes here% DetaiIed expla nati on goes herea=1;b=1;while b<maxa=a+1;b=b*a;endn=a-1;end主程序max=in PUtCPIeaSe in PUt the realmax:');n=fun( max)3.x=-1:0.01:1;y1=cos(acos(x);plot(x,y1, ''% 第一条hold Ony2=cos(

4、2*acos(x);plot(x,y2, ''% 第二条系数矩阵求逆hold Ony3=cos(3*acos(x);plot(x,y3,''%第三条hold Ony4=cos(4*acos(x);pot(x,y4,''% 第四条hold Ony5=cos(5*acos(x);pot(x,y5,''% 第五条六、实验结果1结果：» Ar 6, 3, 4:-2.5,7:8,-4*-3; b-C3;-4：-7;x" inv (A) *bX =0.60007.0000'5.40002结果：实验二MATLAB基本

5、编程实验一、实验目的及要求1. 回顾数字信号处理的主要内容；2. 掌握利用MATLA进行信号处理的方法；3. 了解信号处理工具箱中一些函数的功能；二、实验设备(环境)及要求1. 计算机2. MatIab软件编程实验平台三、实验内容1 . x(n)=8,5,3,11,7,0,-1,4,2,6,5,-5 乞 n 5 是原序列，设 y(n)是原序列受到噪声污染并移位了的序列 y(n) = x(n-3)+w(n),其中w(n)是均值为0,方 差为1的高斯随机序列。计算y(n)与x(n)的相关性。2. 一个线性时不变系统由下面差分方程描述y(n) = 0.8y(n -1) x(n)1) 求频率响应；2)

6、 画出其幅频和相频图；3) 如果系统输入为x(n) = cos(0.02二n)u(n),求出在0" <100内的响 应y(n),并画出输入序列和输出序列的曲线图说明。3. 设计IIR滤波器(带通，三种方法，fs = 2000HZ通带频率500700Hz, 阶数自选，画频率特性并分析比较)四、设计思想1. 由normrnd函数产生高斯序列，判断函数位进行赋值，产生移位序列，通过函数xcorr进行互相关处理2. 用freq函数求出频率相应，通过 abs和angle函数求出幅频和相频图，采用filter函数求出输出序列3.分别采用butterworth滤波器、椭圆滤波器和Cheby滤

7、波器，求出频率特性五、程序代码及注释1. x=8,5,3,11,7,0,-1,4,2,6,5;for n=1:1:11if *4yy( n)=0;else yy( n)=x( n-3);endendW=normrnd(0,1,1,11);y=yy+w;a=xcorr(x,y);2. fs=1000;b=1;a=1,-0.8;h,f=freqz(b,a,512,fs);mag=abs(h);ph=a ngle(h);SubPlOt(4,1,1);plot(f,mag);grid;Xlabel('freque ncy(hz)'); ylabel('mag nitude

8、9;);SubPlOt(4,1,2);plot(f,ph);grid;xlabel('freque ncy(hz)');ylabel('phase');n=1:100;x( n)=cos(0.02*pi* n); y=filter(b,a,x);subplot(4,1,3);Stem(x);subplot(4,1,4);stem(y)3.(1) n=10;Wn=500 7001000;b,a=butter( n,w n);h,f=freqz(b,a,256,1000); mag=abs(h); ph=a ngle(h);SubPlOt(2,1,1);plot(f

9、,mag);grid;XlabeI('freque ncy(hz)'); ylabel('mag nitude'); subplot(2,1,2);plot(f,ph);grid;xlabel('freque ncy(hz)'); ylabel('phase'); n=10;Wn=500 7001000;b,a=ellip( n,0.5,20,w n)h,f=freqz(b,a,256,1000); mag=abs(h); ph=a ngle(h); subplot(2,1,1);plot(f,mag);grid;xlabel(&

10、#39;freque ncy(hz)'); ylabel('mag nitude'); subplot(2,1,2);Pot(f,ph);grid;Xlabel('freque ncy(hz)'); ylabel('phase'); n=10;Wn=500 7001000;b,a=cheby1( n,0.5,w n);h,f=freqz(b,a,256,1000); mag=abs(h); ph=a ngle(h); subplot(2,1,1);plot(f,mag);grid;xlabel('freque ncy(hz)

11、9;); ylabel('mag nitude'); subplot(2,1,2);plot(f,ph);grid;xlabel('freque ncy(hz)'); ylabel('phase');六、实验结果1结果:I FiigLire 1FiIe Edit ieW IniSert IOolS esktop WindOW HIeIP口 匚丄 <3* I a I I 300112o2.结果:I File EdIt VieW InSert ToolS Desktop Vimdow 旦亡IPiP .j j T < ® ? Tl

12、a 3.(1)FigUre 11Ale EClit JJ iew Insert TOOlS DeSldtOP VVindlOW HelP曰UJll钦、紳国凰尿”包5dp=-raLUoo(2)25030035040045050100150200OIH QL A c ® «口 InI骂FigUre 1Flle Edit VieW InSert TOOlS DeSktOP IiyindOW HelP13FlgUrr 1XEile EdiC ylt LnSert IoaIS Qeskrap Aindo Kelp,n-DaafeLiL<f> ® a aAQOOQ

13、 250 frequency(hzIlII l>IIiliIlfcIll i>i>Iii<hiiiI*Il<I*fcI-hbFIFFl>l>PIiIlIiIl! Jl一 _ _ 一 -I _ _ I _ _ZyP H J VVWIIIIIlIII*9lFPiIiiJ i|IIiI|!FIii!l_. _IL_ _ _ _,一 -hiIIilIIIiiIIiIi>IiIi1>*Iil>IIi1>hIhIiIFIFIibIl>l>f0IlKi-ii iiPPIIiIP*Iii IiHhhhFFF I-P - hPi>

14、;P ,bF51实验三 基于MATLA的图像处理实验一、实验目的及要求1了解图像处理的基本概念和功能；2. 掌握利用MATLA进行图像处理的方法；3. 了解图像处理工具箱中一些函数的功能；二、实验设备(环境)及要求1. 计算机2. MatIab软件编程实验平台三、实验内容1. 对barbaran.tif图像进行 DCT变换，分别选取最大的 1/4、1/8、1/16 个变换系数(其余置为0),进行反DCT得到重构图像，显示原图像和各重构 图像并分别计算重构图像的峰值信噪比。2. 对Cameraman.tif图像至少用3种方法进行边缘增强，给出实验结果并 分析。四、设计思想1. 先用imread函

15、数读取，找出其中最大值，采用dct2函数进行dct变换，显示出来，判断小于最大值的变换系数的值赋零，再用idct反变换处理，改变变换系数，分别重构图像，最后再计算重构图像和原图像峰值信噪比2. 采用直方图处理、对比度扩展增强、对数变换这三种方式处理，分别采用histeq、imadjust 和 log 函数。五、程序代码及注释1. grb=imread('barbara n.tif);a=max(grb(:);i=rgb2gray(grb);SubPlOt(2,1,1);imshow(i);j=dct2(i)SubPlOt(2,1,2); imshow(j); b=a4;j(abs(j)

16、<b)=0;k=idct2(j);figure;imshow(k,0 255);A=double(i);B=double(k);m,n=size(i);D=SUm(SUm(A-B).2); MSE=D/(m* n);grb=imread('barbara n.tif); a=max(grb(:);i=rgb2gray(grb); subplot(2,1,1); imshow(i);j=dct2(i)subplot(2,1,2);imshow(j);b=a8;j(abs(j)<b)=0;k=idct2(j);figure;imshow(k,0 255);A=double(i)

17、;B=double(k);m,n=size(i);D=Sum(Sum(A-B).2);MSE=D/(m* n);grb=imread('barbara n.tif); a=max(grb(:);i=rgb2gray(grb);SUbPlOt(2,1,1);imshow(i);j=dct2(i)subplot(2,1,2);imshow(j);b=a16;j(abs(j)<b)=0;k=idct2(j);figure;imshow(k,0 255);A=double(i);B=double(k);m,n=size(i);D=SUm(SUm(A-B).2);MSE=D/(m* n);

18、2. i=imread('camerama n.tif);SUbPlOt(2,2,1);imshow(i);j,t=histeq(i);subplot(2,2,2);imshow(j);k=imadjust(i,0.15 0.9,0 1);subplot(2,2,3);imshow(k);i=double(i);l=log(i+1);SubPlOt(2,2,4);imshow(l,);六、实验结果1.结果：Q FigUfe T- Xfile Edrt jew IrKert IbOISNlndCW IHe-IP Ef3*Z- 口 Xrifk- E< yiew Lnfert I<

19、;ls Qestp JtIndiffW HeIPPdU k H沙歯跟卅尋口目和原图像的峰值噪声比为190.'IAgUrEi 1F Ie di ViSvw nPfir Tjc-Is j d -<和原图像的峰值噪声比为73W Figure 2口XE EIdit Jiew LnSert IQoI5 QCSktOP LViIndoWP Helps U idtV e © - ¾ o S 口2结果:和原图像的峰值噪声比为20JBI FigUre 1XFile £dit Yi亡W Insert Jols D÷sktop yrdov Hlpoaa teg-

20、a>¾uy;三种图像，直方图对比度明显，整体较为中和，对比扩展方式导致黑色更黑，黑 色的细节看不到，而对数增强方式，画面偏白，对比也明显。实验四 基于MATLAB申经网络编程实验、实验目的及要求1. 了解神经网络的基本概念和原理；2. 掌握用MATLA实现神经网络的思路和方法；3. 了解神经网络工具箱函数的功能。二、实验设备(环境)及要求1. 计算机2. MatIab软件编程实验平台三、实验内容1、利用2输入的MP模型是否可实现“与”、“或”、“非”和“异或”功能? 若能，编程实现并给出实现各种功能的权值和实验结果，若不能，请说明原 因。2、用人工神经网络拟合函数_0 I5xf

21、 (X)= Sin(x) + 0.46e .cos(1.2x), OC x c 4兀说明：1)网络结构为三层(输入层、1个隐层和输出层)2 )获取两组数据，一组作为训练集，一组作为测试集3 )用训练集训练网络4 )用测试集检验训练结果四、设计思想1. 使用MP模型解决此问题，选择初值，计算误差，更新全值，对所有模式 进行所有训练，直到输出满足。2. 使用三层BP网络来完成函数逼近。因为函数为一元函数，所以输入层有一个输入，输出层有一个输出，隐层神经元的个数为10。五、程序代码及注释1. p=0 0 1 1;0 1 0 1; % 与t=0 0 0 1;d=0 1 0 1;W仁 0.2;w2=0.

22、5;yita=-0.1;theta=0.1;while isequal(t,d)for i=1:1:4W仁 w1+yita*(d(j)-t(j)*p(1,j); w2=w2+yita*(d(j)-t(j)* p(2,j); theta=theta+yita*(d(j)-t(j); d(j)=w1*p(1,j)+w2*p(2,j)+theta; d(j)=rou nd(d(j);endendp=0 0 1 1;0 1 0 1;%或t=0 1 1 1;d=0 1 0 1;W仁 0.2;w2=0.5;yita=-0.1;theta=0.1;while isequal(t,d)for j=1:1:4W

23、仁 w1+yita*(d(j)-t(j)*p(1,j); w2=w2+yita*(d(j)-t(j)* p(2,j); theta=theta+yita*(d(j)-t(j); d(j)=w1*p(1,j)+w2*p(2,j)+theta; d(j)=rou nd(d(j);endendp=0 0 1 1;0 1 0 1;%异或t=0 1 1 0;d=0 1 0 1;W仁 0.2;w2=0.5;yita=-0.1;theta=0.1;while isequal(t,d)for j=1:1:4W仁 w1+yita*(d(j)-t(j)*p(1,j); w2=w2+yita*(d(j)-t(j)*

24、 p(2,j); theta=theta+yita*(d(j)-t(j); d(j)=W1*p(1,j)+W2*p(2,j)+theta; d(i)=ro Un d(d(i);end end2. x=0:0.01:4*pi;%画出模拟曲线y=si n(x)+0.46*exp(-0.15*x).*cos(1.2*x); figure(1);PlOt(X,y,'g');hold on;N=50;%初始化预测集样本数P仁0:0.2:4*pi;%训练集样本点P2=(4*pi)*ra nd(1,N);% 预测集样本点%计算训练集和预测集对应的目标向量T1,T2T1= Sin (P1)+0

25、.46*exp(-0.15*P1).*cos(1.2*P1); T2= Sin (P2)+0.46*exp(-0.15*P2).*cos(1.2*P2);%建立网络n et = n ewff(P1,T1,1 6);%设置训练参数n et.tra in Param.show = 50;n et.tra in Para m.lr = 0.05;n et.tra in Param.mc = 0.95;n et.tra in Param.epochs = 400;n et.tra in Param.goal = 1e-8;n et,tr=trai n(n et,P1,T1); %训练网络A=sim(

26、net,P2); %对预测集进行预测PlOt(P1,T1,'r+',P2,A,'b.');六、实验结果1. 与：W仁0.2 w2=0.4 theta=0或：w1=0.3 w2=0.5 theta=0.2无法实现非功能，由于非只需要一个输入异或：W仁0.1w2=0.1 theta=0.22.M Figure 1己Pd 4食Q甥樹晏N -良I因I 口OeOeffJfA1« Edit View InSert TOOIS DeSktGP Window HeIP0.40.20 2-0.4102468'1012140.C图中的蓝色曲线是所需模拟函数的曲线，

27、红色 +'为训练集的样本点， 而蓝色的.是预测集P2经过网络预测的结果A,即（P2, A）。我们可以 发现模拟的结果很好。实验五 MATLAB的综合实验、实验目的及要求培养学生利用MatIab解决专业问题的能力、实验设备（环境）及要求1.计算机2.MatIab软件编程实验平台三、实验内容（三题中选择一题）1、实现一个数字信号处理的仿真系统，通过该系统可以方便的验证信号处理 所学的基本内容，要求程序具有界面并实现以下功能：1）能产生或得到各种数字信号（Sin、方波、三角波、语音、噪声及其叠 加）；2）具有根据需求设计低通、高通、带通和带阻滤波器的功能；3）具有对数字信号进行DFT、DCT

28、和DWT变换和经典功率谱估计功能；4）产生或输入数字信号，显示时域、变换域及功率谱估计波形，在这个 信号上叠加噪声，通过滤波器对其降噪，得到输出信号的频域特性和 时间序列。2、 基于数字图像处理，实现一个汽车标志定位和分割的仿真系统。要求系统具有界面并实现以下功能：1）读入包含汽车标志的汽车图像；2）预处理及汽车标志定位；3）多种图像分割功能；4）形态学滤波等功能；3、结合自己的研究方向，实现一个神经网络应用的实例，要求用至少3种方 法实现（可以包括不同神经网络方法或经典方法），并对实验结果进行比 较分析。四、设计思想1.实现一个数字信号处理的仿真系统根据数字信号处理的仿真系统要求，设计界面框

29、图，安排界面布局图像二值化用GUl按照设计进行界面交互组件添加并设置组件属性对回调函数进行编程，实现题设功能调整界面并调试程序五、程序代码(界面除外)及注释function VarargOut = renliang(varargin)% Begin initialization code - DO NOT EDIT gui_S in glet on = 1;gui State = StrUCtegUi Name',mfilename, .'gui_S in glet on', gui_S in glet on,'gui_Ope nin gFc n', r

30、en Iia ng_OPe nin gFc n,'gui_OUtPUtFC n', ren Iia ng_OUtPUtFC n,'gui_LayOUtFC n',. 'gui_Callback', );if nargin && ischar(varargin1)gui_State.gui_Callback = str2fu nc(varargi n1);endif n argoutvarargout1: nargout = gui_mai nfcn( gui_State,Vararg in :);elsegui_mai nfcn(

31、 gui_State, Varargi n:);end% End initialization code - DO NOT EDlT% - EXeCUteS just before renIiang is made visible.function ren Iia ng _OPe nin gFc n( hObject, eve ntdata, han dles, Varargi n)% Choose defaUlt comma nd Iine OUtPUt for ren Iia nghan dles.oUtpUt = hObject;% UPdate han dles StrUCtUregU

32、idata(hObject,ha ndles);% UIWAIT makes ren Iia ng Wait for USer respo nse (See UIRESUME)% - OUtPUtS from this fUn ctio n are retUrned to the comma nd line.fUnction VarargOUt = renIiang _OUtPUtFCn(hObject, eventdata, han dles)% Get defaUlt comma nd Iine OUtPUt from han dles StrUCtUreVarargoUt1 = han

33、dles.oUtpUt;% - EXeCUteS on bUtt on PreSS in CheCkbox1.fUn ctio n CheCkbox1_CaIIback(hObject, eve ntdata, han dles)%正弦波输入%- EXeCUteS dUring ObjeCt Creation, after Settingall PrOPerties.fUnction CheCkbox1_CreateFC n( hObject, eve ntdata, han dles)fUnction edit1_Callback(hObject, eventdata, handles)%正

34、弦波幅值%- EXeCuteS during ObjeCt Creation,after Settingallfun ctio n edit1_CreateFC n(hbject, eve ntdata, han dies) if ispc && isequal(get(hObject,'BackgroUndColor'), get(0,'defaultUico ntrolBackgrou ndColor') set(hObject,'Backgro Un dColor','white');endfunction

35、edit6_Callback(hObject, eventdata, handles)%正弦波频率%- EXeCUteS during ObjeCt Creation,after Settingallfunction edit6_CreateFCn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroUndColor'), get(0,'defaultUico ntrolBackgrou ndColor') set(hObject,'Backgro Un

36、dColor','white');end% - EXeCUteS on butt on PreSS in CheCkbox1.fun ctio n CheCkbox2_Callback(hObject, eve ntdata, han dles)%方波输入%- EXeCUteS during ObjeCt Creation,after Settingallfunction CheCkbox2_CreateFC n( hObject, eve ntdata, han dles)function edit2_Callback(hObject, eventdata, hand

37、les)%方波幅值%- EXeCUteS during ObjeCt Creation,after Settingallfunction edit2_CreateFCn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroUndColor'), get(0,'defaultUico ntrolBackgrou ndColor') set(hObject,'Backgro Un dColor','white');endfun ctio

38、 n edit7_Callback(hObject, eve ntdata, han dles)%方波频率%- EXeCUteS during ObjeCt Creation,after Settingallfunction edit7_CreateFCn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroUndColor'), get(0,'defaultUico ntrolBackgrou ndColor') set(hObject,'Backgro

39、 Un dColor','white');PrOPerties.PrOPerties.PrOPerties.PrOPerties.PrOPerties.end% - EXeCuteS On butt On PreSS in CheCkboxi.fun ctio n CheCkbox3_Callback(hObject, eve ntdata, han dies)%角波输入%- EXeCUteS during ObjeCt Creation,after Settingallfunction CheCkbox3_CreateFC n( hObjeCt, eve ntdata

40、, han dies)fun ctio n edit3_CaiibaCk(hObject, eve ntdata, han dies)%角波幅值%- EXeCUteS during ObjeCt Creation,after Settingallfun ctio n edit3_CreateFC n(hObject, eve ntdata, han dies) if ispc && isequal(get(hObject,'BackgroUndColor'), get(0,'defaultUico ntrolBackgrou ndColor')

41、set(hObject,'Backgro Un dColor','white');endfun ctio n edit8_CaiibaCk(hObject, eve ntdata, han dies)%角波频率%- EXeCUteS during ObjeCt Creation,after Settingallfun ctio n edit8_CreateFC n( hObject, eve ntdata, han dies) if ispc && isequal(get(hObject,'BackgroUndColor'), g

42、et(0,'defaultUico ntrolBackgrou ndColor') set(hObject,'Backgro Un dColor','white');end% - EXeCUteS on butt on PreSS in CheCkboxi.fun ctio n CheCkbox4_CaiibaCk(hObject, eve ntdata, han dies)%语音信号输入%- EXeCUteS during ObjeCt Creation,after Settingallfunction CheCkbox4_CreateFC n

43、( hObject, eve ntdata, han dies)fun ctio n edit4_CaiibaCk(hObject, eve ntdata, han dies)%语音信号幅值%- EXeCUteS during ObjeCt Creation,after Settingallfun ctio n edit4_CreateFC n(hObject, eve ntdata, han dies) if ispc && isequal(get(hObject,'BackgroUndColor'), get(0,'defaultUico ntrol

44、Backgrou ndColor') set(hObject,'Backgro Un dColor','white');PrOPerties.PrOPerties.PrOPerties.PrOPerties.PrOPerties.end% - EXeCuteS On butt On PreSS in CheCkboxi.fun ctio n CheCkbox5_Callback(hObject, eve ntdata, han dies)%白噪声输入%- EXeCUteS during ObjeCt Creation, after Settingall

45、PrOPerties.function CheCkbox5_CreateFC n( hObjeCt, eve ntdata, han dies)fun ctio n edit5_CaiibaCk(hObject, eve ntdata, han dies)%白噪声幅值%- EXeCUteS during ObjeCt Creation, after Settingall PrOPerties.fun ctio n edit5_CreateFC n(hObject, eve ntdata, han dies)if ispc && isequal(get(hObject,'

46、BackgroUndColor'), get(0,'defaultUico ntrolBackgrou ndColor')set(hObject,'Backgro Un dColor','white');end% - EXeCUteS on butt on PreSS in PUShbUtt on1.function PUShbUtt on 1_CaiibaCk(hObject, eve ntdata, han dies)%信号输入控制global N;%信号取样数N=110033;global n;%信号取样点n=1:N;Sgn 1=g

47、et(ha ndles.checkbox1,'value');%5Sgn 2=get(ha ndles.checkbox2,'value');Sgn 3=get(ha ndles.checkbox3,'value'); sgn4=get(ha ndles.checkbox4,'value');Sgn 5=get(ha ndles.checkbox5,'value');rag1=str2 nu m(get(ha ndles.edit1,'stri ng');%5 rag2=str2 nu m(get(

48、ha ndles.edit2,'stri ng');rag3=str2 nu m(get(ha ndles.edit3,'stri ng'); rag4=str2 nu m(get(ha ndles.edit4,'stri ng');rag5=str2 nu m(get(ha ndles.edit5,'stri ng');frq1= str2 nu m(get(ha ndles.edit6,'stri ng');%5 frq2=str2 nu m(get(ha ndles.edit7,'stri ng

49、9;); frq3=str2 nu m(get(ha ndles.edit8,'stri ng');类信号选择与否类信号幅值输入类信号频率输入x1=si n( frq1*pi* n/N*2);%产生正弦波x2=SqUare(frq2*Pi* n/N*2);%产生方波x3=SaWtOOth(frq3*pi* nN*2,0.5)%产生三角波global x4;%读取语音信号global fs;x4,fs=wavread('speech_dft.wav'); x4=x4'x5=2*ra nd(1,N)-1;%产生白噪声global x;%选择信号叠加的输入信号

50、x=rag1*sg n1.*x1+rag2*sg n2.*x2+rag3*sg n3.*x3+rag4*sg n4.*x4+rag5*sg n5.* x5;axes(ha ndles.axes1);plot(x)title('输入信号时间序列)；% - EXeCUteS on butt on PreSS in radiobutt on1.fun ctio n radiobutt on 1_Callback(hObject, eve ntdata, han dles)%选择DFT变换% - EXeCUteS on butt on PreSS in radiobutt on2.fun ct

51、io n radiobutt on 2_Callback(hObject, eve ntdata, han dles) %选择DCT变换% - EXeCUteS on butt on PreSS in radiobutt on3.fun ctio n radiobutt on 3_Callback(hObject, eve ntdata, han dles) %选择DW变换% - EXeCUteS on butt on PreSS in PUShbUtt on2. function PUShbUtt on 2_Callback(hObject, eve ntdata, han dles) %对

52、输入信号进行 DFT DCT和DWT变换global N;%全局变量global n;global x;slect1=get(ha ndles.radiobutto n1,'value');%3种变换选择slect2=get(ha ndles.radiobutton 2,'value');slect3=get(ha ndles.radiobutton 3,'value');if slect1% 选择DFT变换 y=fft(x,N);%DFT 变换 axes(ha ndles.axes2); plot(abs(y(1:(N+1)/2) title(

53、'输入信号的DFT变换');else if SleCt2%选择 DCT变换y=dct(x);%DCT变换axes(ha ndles.axes2);plot(abs(y)title('输入信号的DCT变换');else if slect3%选择 DW变换cA,cD=dwt(x,'db1');%DWT变换%近似系数CA1和细节系数CD1axes(ha ndles.axes2);PlOt(1:(N+1)/2,abs(cA),1:(N+1)/2,abs(cD) lege nd('cA','cD');title('输

54、入信号的DW变换');endendend% - EXeCUteS On butt On PreSS in radiobutt on4.fun ctio n radiobutt on 4_Callback(hObject, eve ntdata, han dles) %选择低通滤波set(ha ndles.edit9,'E nable','off);%截止低频输入框灰暗显示set(ha ndles.edit10,'E nable',' on');%截止高频输入框使能% - EXeCUteS on butt on PreSS in r

55、adiobutt on5.fun ctio n radiobutt on 5_Callback(hObject, eve ntdata, han dles)%选择低通滤波一set(ha ndles.edit9,'E nable',' on');%截止低频输入框使能set(ha ndles.edit10,'E nable','off)%截止高频输入框灰暗显示% - EXeCUteS on butt on PreSS in radiobutt on6.fun ctio n radiobutt on 6_Callback(hObject, eve ntdata, han dles)%选择低通滤波一set(ha ndles.edit9,'E nable',' on');%截止低频输入框使能set(ha ndles.edit10,'E nable',' on');%截止高频输入框使