addfs4.m

很多matlab的源代码

%FOURIER Graphics demo of Fourier series expansion.
%NOTE: a is from a=fseries(...), ty is window type specified in fseries


% 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


vx=matverch;if vx<4,disp('This routine is NOT for MATLAB 3.x '),return,end

clf,[ma,na]=size(a);Nmax=ma-1;N=min(9,Nmax);colormap(cool);vw=[-345 30]; 
t = 0:.01:1;ak=a(:,2);bk=a(:,3);lt=length(t);x0 = ak(1)*ones(1,lt);
if exist('ty')==0,ty='rect';else,if isempty(ty),ty='rect';end,end
if vx<5, ytk='yticklabels';else,ytk='yticklabel';end

s1 = ['x=x0;y = zeros(N+1,lt);y(1,:)=x0;' ...
	'win=window(ty,2*N+1);win=win(N+1:2*N+1);ys=x0*win(1);'...
	'for k=2:N+1,m=k-1;sm=ak(k)*cos(2*pi*m*t)+bk(k)*sin(2*pi*m*t);' ...
	'ys=ys+sm*win(k);x=x+sm;y(k,:)=x;end,subplot(2,1,1),' ...
	'zl=min(min(y));zh=max(max(y));ax3=[0 1 0 N zl zh];'...
	'ax2=[0 1 zl zh];plot(t,y(N+1,:),t,ys),axis(ax2);' ...
	'title(''Final & smoothed Reconstruction''),' ...	
	'nn=int2str(N);ln=length(nn);n0=[blanks(ln-1) ''0''];'...
	'subplot(2,1,2),surf(t,(0:N),flipud(y)),axis(ax3);view(vw),'...
	'set(gca,''ytick'',[0 N],ytk,[nn;n0],''box'',''on''),'...
	'grid,title([int2str(N) '' harmonics''])'];

eval(s1),disp(' '),disp('USE POPUP MENU TO CHANGE NUMBER OF TERMS'),disp(' ')

s=['N = get(HH,''Value'');' s1];

st='1';for j=2:Nmax,st=[st '|' int2str(j)];end
HH=uicontrol('style','popupmenu','string',st);
set(HH,'unit','norm','pos',[0.005 0.895 .08 .05],'value',N,'callback',s);