xian.m

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

% xian(haxes1,haxes2,haxes3,haxes4) 绘制图象
%function []=xian(haxes1,haxes2,haxes3,haxes4,x1,x2)

%x1=get(gcbo,'string');
%x1=str2num(x1);
%x2=get(gcbo,'string');
%x2=str2num(x2);
x1=input('x1=');
%x1=[0 1 2 3 5]
%x2=[1 2 3 5 6]
x2=input('x2=');
haxes1=axes('position',[0.10,0.05,0.63,0.18]);
haxes2=axes('position',[0.10,0.29,0.63,0.18]);
haxes3=axes('position',[0.10,0.53,0.63,0.18]);
haxes4=axes('position',[0.10,0.77,0.63,0.18]);

lx1=length(x1);
lx2=length(x2); 

lmax=max(lx1,lx2); 
if lx1>lx2 nx1=0; nx2=lx1-lx2;               
   elseif lx1<lx2 nx1=lx2-lx1; nx2=0;
else nx1=0; nx2=0; 
end 

%axes(haxes4);hold off
dt=1;              % dt 运算时间间隔 ; 
lt=lmax;           % 取长者为补零长度基准 
% 先将 x 补得与 h 同长，再两边补上同长度的零 
x=[zeros(1,lt),x1,zeros(1,nx1),zeros(1,lt)]; %设置x为3倍卷积长度
t1=(-lt+1 : 2*lt)*dt;                        % total length of t1 is 3*lmax 
% 将 h 补得与 u 同长，再前面补上 2*lt 长度的零 
stem(t1,x,'.r');
%xmin=min(t1),xmax=max(t1),ymin=min(x),ymax=max(x);
axis([min(t1),max(t1),min(x),max(x)]), hold on; 
%axis([xmin xmax 0 1]), hold on; 

axes(haxes3);hold off
%h=[zeros(1,2*lt),x2,zeros(1,nx2)]; 
%hf=fliplr(h); % 将 h 的左右翻转 ， 记为 hf 
%y=zeros(1,3*lt); % y 数组初始化 
for k=0: 2*lt 
axes(haxes3);hold off
h=[zeros(1,2*lt),x2,zeros(1,nx2)]; 
hf=fliplr(h); % 将 h 的左右翻转 ， 记为 hf 
y=zeros(1,3*lt); % y 数组初始化 
p=[zeros(1,k),hf(1: end-k)]; % 使 hf 向右循环移位 ,hf(1: end-k)=hf(1:3*lt-k) 
%y1=x.*p*dt; % 使输入和翻转移位的脉冲过渡函数逐项相乘 ， 再乘 dt 
%yk=sum(y1); % 相加，等效于积分 
%y(k+lt+1)=yk; % 将结果放入数组 y 
if k==0
   stem(t1,hf,'.r');
    axis([-lt*dt, 2*lt*dt, min(p), max(p)]); 
    pause
end
stem(t1,p,'.r') 
axis([-lt*dt, 2*lt*dt, min(p), max(p)])

axes(haxes2);hold off
p=[zeros(1,k),hf(1: end-k)]; % 使 hf 向右循环移位 ,hf(1: end-k)=hf(1:3*lt-k) 
y1=x.*p*dt; % 使输入和翻转移位的脉冲过渡函数逐项相乘 ， 再乘 dt 
%yk=sum(y1); % 相加，等效于积分 
%y(k+lt+1)=yk; % 将结果放入数组 y 
stem(t1,y1,'.r')
%axis([-lt,2*lt,min(y1),max(y1)]), 

axes(haxes1);
p=[zeros(1,k),hf(1: end-k)]; % 使 hf 向右循环移位 ,hf(1: end-k)=hf(1:3*lt-k) 
y1=x.*p*dt; % 使输入和翻转移位的脉冲过渡函数逐项相乘 ， 再乘 dt 
yk=sum(y1); % 相加，等效于积分 
stem(k*dt,yk,'.r')
%axis([-lt*dt,2*lt*dt,floor(min(y)+eps),ceil(max(y)+eps)]),
hold on
pause(1) 
end 
hold on