caiyang.m

数字信号处理 各种滤波器的设计

%caiyang(haxes1,haxes2,haxes3,haxes4,haxes5,haxes6)
%function []=caiyang(haxes1,haxes2,haxes3,haxes4,haxes5,haxes6)
figure(3);
haxes1=axes('position',[0.05,0.70,0.40,0.21]);%position后面的矩阵四个数据分别表示组建的坐标位置的左边、底部、宽、高
                                              %整个图形窗口的左下为（0，0），右上为（1，1）
haxes2=axes('position',[0.53,0.70,0.40,0.21]);
haxes3=axes('position',[0.05,0.40,0.40,0.21]);
haxes4=axes('position',[0.53,0.40,0.40,0.21]);
haxes5=axes('position',[0.05,0.10,0.40,0.21]);
haxes6=axes('position',[0.53,0.10,0.40,0.21]);

global m
%fs=get(gcbo,'string');
fs=input('频率等于：');
%fs=str2num(fs);
m=1/fs;
axes(haxes1);
N=5000;
dt=0.2/N;
t=[-N:N]*dt;
f=100*(sinc(100*t/4)).*(sinc(100*t/4))/2;
plot(t,f,'linewidth',1.5),grid

axes(haxes2);
N=5000;
dt=0.2/N;
t=[-N:N]*dt;
f=100*(sinc(100*t/4)).*(sinc(100*t/4))/2;
wf=5000;
Nf=512;
w1=linspace(0,wf,Nf);
F1=f*exp(-j*t'*w1);
w=[-fliplr(w1),w1(2:Nf)];
F=[fliplr(F1),F1(2:Nf)];
plot(w,abs(F),'linewidth',1.5),grid

axes(haxes3);
t1=[-0.2:m:0.2];
y=ones(1,length(t1));
stem(t1,y,'.'),grid
axis([-0.2,0.2,0,1.3]);

axes(haxes4);
t1=[-0.2:m:0.2];
y=ones(1,length(t1));
wf=5000;
Nf=512;
w1=linspace(0,wf,Nf);
F1=y*exp(-j*t1'*w1);
w=[-fliplr(w1),w1(2:Nf)];
F=[fliplr(F1),F1(2:Nf)];
plot(w,abs(F)/max(abs(F)),'linewidth',1.5),grid

axes(haxes5);
t1=[-0.2:m:0.2];
f=100*(sinc(100*t1/4)).*(sinc(100*t1/4))/2;
stem(t1,f,'.'),grid

axes(haxes6);
t1=[-0.2:m:0.2];
f=100*(sinc(100*t1/4)).*(sinc(100*t1/4))/2;
wf=5000;
Nf=512;
w1=linspace(0,wf,Nf);
F1=f*exp(-j*t1'*w1);
w=[-fliplr(w1),w1(2:Nf)];
F=[fliplr(F1),F1(2:Nf)];
plot(w,abs(F)/max(abs(F)),'linewidth',1.5),grid