caiyang.asv

数字信号处理_matlab实现程序,上面有全套的实验代码

% 数字信号处理教程——MATLAB释义与实现%caiyang   演示采样定理，时域采样，频谱周期延拓
%同时演示采样频率小于2*fc时，产生的混叠效应
clf reset
set(gcf,'unit','normalized','position',[0.1,0.2,0.5,0.5]);
set(gcf,'defaultuicontrolunits','normalized');
haxes1=axes('position',[0.05,0.2,0.4,0.2],...
            'xlim',[-5,5],'ylim',[-5,5]);
haxes2=axes('position',[0.5,0.2,0.4,0.2],...
            'xlim',[-1,1],'ylim',[0,1]);
haxes3=axes('position',[0.05,0.45,0.4,0.2],...
            'xlim',[-5,5],'ylim',[-5,5]);
haxes4=axes('position',[0.5,0.45,0.4,0.2],...
            'xlim',[-1,1],'ylim',[0,1]);
haxes5=axes('position',[0.05,0.7,0.4,0.2],...
            'xlim',[-20,20],'ylim',[-5,5]);
        %绘出初始时的时域信号
        t=-0.05:0.001:0.05;
        fc=100;   %这里假设三角波的频带宽度为100Hz
        y=fc/2*sinc(fc*t/4).*sinc(fc*t/4);
        plot(t,y);
        
haxes6=axes('position',[0.5,0.7,0.4,0.2],...
            'xlim',[-1,1],'ylim',[0,1.1]);
        %绘出频谱
x1=1:-0.1:0;  x2=0:0.1:0.9;
x3=[x2,x1];
k=-10:10;
df=0.01;%这里假设fc=100,df代表10HZ 所以fc占用0.1
f=k*df;
plot(f,x3);axis([-1,1,0,1.1]);


htext1=uicontrol(gcf,'style','text',...
                 'position',[0.05,0.05,0.6,0.05],...
                 'backgroundcolor','white',...
                 'fontsize',10,...
                 'visible','on',...
                 'string','信号频带宽度为100Hz，输入采样频率（100<fs<900Hz）');
 htext2=uicontrol(gcf,'style','text',...
                 'position',[0.2,0.92,0.15,0.05],...
                 'fontsize',10,...
                 'string','时    域');
 htext3=uicontrol(gcf,'style','text',...
                  'position',[0.65,0.92,0.15,0.05],...
                  'fontsize',10,...
                  'string','频   域');
 hedit=uicontrol('style','edit',...
                 'position',[0.7,0.05,0.2,0.05],...
                 'backgroundcolor','white',...
                 'visible','on');
 set(hedit,'CallBack',['call_caiyang1(haxes4);',...
             'call_caiyang3(haxes3);',...
             'call_caiyang4(haxes1,fc);',...
            'call_caiyang2(haxes2);']);