demofs.m

很多matlab的源代码

% ADSP Toolbox: Version 2.0 
% For use with "Analog and Digital Signal Processing", 2nd Ed.
% Published by PWS Publishing Co.
%
% Ashok Ambardar, EE Dept. MTU, Houghton, MI 49931, USA
% http://www.ee.mtu/faculty/akambard.html
% e-mail: [e-mail]
% Copyright (c) 1998


%FOURIER SERIES DEMO
echo off,clc
vx=matverch;

n=-1;while ((n ~= 0)), clc,help fsmenu
n = input('ENTER A NUMBER FROM 1 to 6 [or 0 TO QUIT] : ');
if n==0,
if vx < 4, eval('clg');else,eval('clf');end
%break
return
end
if n>0 & n<7
d=0;while (d<=0 | d>1),d=input('Duty ratio [Pulse width/Period](<=1): ');end
if n==1,if d<1,x='.5*(t==d)+(t>=0 & t<d)';else,x='1+0*t';end,end
if n==2,if d<1,x='0.5*(t==d)+t.*(t>=0 & t<d)/d';else,x='t';end,end
if n==3,x='2*t.*(t>0 & t<=d/2)/d+(2-2*t/d).*(t>d/2 & t<=d)';end
if n==4,c=input('Duration of flat portion: ');
while c>=d|c==0,disp(['Flat portion must be >0 & <' num2str(d)])
c=input('Duration of flat portion: ');end,t1=.5*(d-c);
t2=t1+c;x='t.*(t>0&t<=t1)/t1+(t>t1&t<=t2)+(t-d).*(t>t2&t<=d)/(t2-d)';end
if n==5,c=input('# of cycles over pulse width: ');
if d<1,x='sin(2*pi*t*c/d).*(t>=0 & t<d)+.5*sin(2*pi*c)*(t==d)';
else,x='sin(2*pi*t*c/d)';end,end
if n==6,c=input('Time constant: ');
if d<1,x='.5*exp(-d/c)*(t==d)+exp(-t/c).*(t>=0 & t<d)';
else,x='exp(-t/c)';end,end
td=input('Shift (+ or -) from origin (HIT ENTER IF 0): ');
if isempty(td),td=0;end,td=rem(td,1);
ty=input('Smoothing window IN SINGLE QUOTES e.g ''hamming'' else ENTER: ');
disp(' '),
vr=version;del=1024;
%if vr(1)=='S',del=512;end %%Check this???

t=0:1/del:1;yt=eval(x);
if isempty(ty),[a,b]=fseries(yt,[1 td]);
else,[a,b]=fseries(yt,[1 td],[1 3 5 32],ty);end
if isempty(b),p='Convergence rate = 1/';
if (n==1|n==2|n==6),disp([p 'k'])
elseif (n==3|n==4), disp([p '(k^2)'])
else
if c/0.5==round(c/0.5),disp([p '(k^2)']),else,disp([p 'k']),end
end
end
disp(' '),disp('NOTE: These results to k=32 are stored in the 32x7 matrix a')
vx=matverch;
if vx>=4,if exist('addfs4')==2,eval('addfs4'),end,end
disp(''),disp('STRIKE A KEY TO RETURN TO THE FS MENU'),pause
subplot,hold off,end,end,
clear b c d n td x t yt
clc