📄 siga.m
字号:
%if(freq1<0.05)freq1=0.05;end
%if(freq>0.95)freq1=0.95;end
% if(amp1>1.0)amp1=1.0;end
if (val == 1),
%if t==0 f=1;
%else
f=amp1*sin(freq1*t)./(t*pi);
% end;
elseif (val == 2), % square wave
for i=1:N
if abs(i-N/2)<=N/(6*freq1)
f(i)=amp1;
else
f(i)=0;
end
end
%tt=freq1*t*2*pi;
% tmp=rem(tt,2*pi);
%f=amp1*(2*rem((tt<0)+(tmp>pi | tmp<-pi)+1,2)-1);
elseif (val == 3), % sawtooth
for i=1:N
if i==N/2, f(i)=1;
else f(i)=0;
end;
end
%tt=freq1*t*2*pi;
% f=amp1*((tt < 0) + rem(tt,2*pi)/2/pi - .5)*2;
elseif (val == 4), % self-defined waveform
f=eval(fun_str);
end;
set(time_line,'YData',f);
F=fft(f,2*M);
F=F(1:M);
%w=(0:M-1)*ws/(2*M);
FF=abs(F);
% ind=find(FF<min_dB);
% FF(ind)=NaN*ones(size(ind)); % put NaN's in where
% min_dB shows up - this is to work around no clipping in xor mode
% freq_line=plot(w1(1:N/2),FF(1:N/2),w2(1:N/2),FF(1:N/2),'EraseMode','xor');
% axis([-100 100 0 40]);
set(freq_line,'YData',FF);
set(freq_field,'String',num2str(freq1));
elseif strcmp(action,'up'),
pt=get(gca,'currentpoint');
x=pt(1,1);
y=pt(1,2);
set(gcf,'WindowButtonMotionFcn','');
set(gcf,'WindowButtonUpFcn','');
% u=get(gcf,'userdata');
freq=siga_DAT(1);
amp=siga_DAT(2);
fixed_x=siga_DAT(6);
fixed_y=siga_DAT(7);
amp1=y/fixed_y*amp;
if (abs(amp1)>1.0),
amp1=1.0*sign(amp1);
end;
if (abs(amp1)<0.05),
amp1=0.05*sign(amp1);
end;
freq1=fixed_x/x*freq;
set(siga_DAT(8),'linestyle','-','markersize',8);
siga_DAT(1)=freq1; % set amplitude and frequency
siga_DAT(2)=amp1;
% set(gcf,'userdata',u);
siga('redraw');
elseif strcmp(action,'done'),
% close the figure window that is showing the window fnction:
% u = get(gcf,'userdata');
show_window=findobj('Tag','window_window');
if ( ~isempty(show_window) & siga_DAT(12)~=-1 )
close(siga_DAT(12));
end;
close(gcf);
clear global siga_DAT
clear global ADDIT_DAT
clear global fun_str
clear show_window
elseif strcmp(action,'redraw'),
% recomputes time and frequency waveforms and updates display
% u = get(gcf,'userdata');
freq=siga_DAT(1);
amp=siga_DAT(2);
N=siga_DAT(3);
M=siga_DAT(4);
min_dB=siga_DAT(5);
time_line=siga_DAT(8);
freq_line=siga_DAT(9);
freq_field=siga_DAT(10);
popup=siga_DAT(11);
t=siga_DAT(14:14+N-1);
window=siga_DAT(14+N:14+N+N-1);
val = get(popup,'Value');
% if(freq<0.05)freq=0.05;end
%if(freq>10)freq=10;end
% if(amp>1.0)amp=1.0;end
if (val == 1),
%if t==0,f=1
% else
f=amp*sin(freq*t)./(t*pi);
% end;
elseif (val == 2), % square wave
for i=1:N
if abs(i-N/2)<=N/(6*freq)
f(i)=amp;
else
f(i)=0;
end
end
%tt=freq*t*2*pi;
% tmp=rem(tt,2*pi);
% f=amp*(2*rem((tt<0)+(tmp>pi | tmp<-pi)+1,2)-1);
elseif (val == 3), % sawtooth
for i=1:N
if i==N/2, f(i)=1;
else f(i)=0;
end
end
%tt=freq*t*2*pi;
%f=amp*((tt < 0) + rem(tt,2*pi)/2/pi - .5)*2;
elseif (val == 4), % self-defined waveform
f=eval(fun_str);
end;
set(time_line,'YData',f);
F=fft(f,2*M);
F=F(1:M);
%w=(0:M-1)*ws/(2*M);
FF=abs(F);
% ind=find(FF<min_dB);
% FF(ind)=NaN*ones(size(ind)); % put NaN's in where
% min_dB shows up - this is to work around no clipping in xor mode
% freq_line=plot(w1(1:N/2),FF(1:N/2),w2(1:N/2),FF(1:N/2),'EraseMode','xor');
% axis([-100 100 0 40]);
set(freq_line,'YData',FF);
set(freq_field,'String',num2str(freq));
drawnow;
elseif strcmp(action,'setwindow'),
% u = get(gcf,'userdata');
winHndl = ADDIT_DAT;
in1 = get(winHndl,'Value');
in2 = 30;
N=siga_DAT(3);
show_window=siga_DAT(14+N+N);
if (in1==1),
window = boxcar(N);
elseif (in1==2),
window = triang(N);
elseif (in1==3),
window = hanning(N);
elseif (in1==4),
window = hamming(N);
elseif (in1==5),
window = chebwin(N,30);
elseif (in1==6),
window = kaiser(N,4);
end;
siga_DAT(14+N:14+N+N-1)=window;
% set(gcf,'userdata',u);
siga('redraw');
if get(show_window,'value'),
siga('showwin');
end;
elseif strcmp(action,'showwin'),
% u=get(gcf,'userdata');
oldfig=gcf;
N=siga_DAT(3);
t=siga_DAT(14:14+N-1);
window=siga_DAT(14+N:14+N+N-1);
if (siga_DAT(12)==-1),
siga_DAT(12)=figure('NumberTitle','off', ...
'Name','窗口', 'menubar','none',...
'Tag','window_window');
axes('Position',[.15 .62 .8 .3]);
line1=plot(t,window);
title('窗口函数');
xlabel('时间（秒）');
grid;
ylabel('窗口');
axes('Position',[.15 .2 .8 .3]);
W=fft(window,1024);
line2=plot((0:(1/1024):(.5-(1/1024)))*N,abs(W(1:512)));
set(gca,'xlim',[0 N/2]);
xlabel('频率 (Hz)');
ylabel('幅度 (dB)');
grid;
windclose=uicontrol('Style','Pushbutton',...
'Position',[.85 .02 .12 .08],...
'Units','normalized',...
'String','关闭',...
'Callback',[...
'THEHAND=get(gcf,''userdata'');',...
'close;',...
'figure(THEHAND(1));',...
'global siga_DAT;',...
'siga_DAT(12)=-1;',...
'clear THEHAND siga_DAT;']);
set(gcf,'userdata',[oldfig line1 line2]);
% set(oldfig,'userdata',u);
clear line1 line2 W;
else
show_window=findobj('Tag','window_window');
if isempty(show_window)
siga_DAT(12)=-1;
siga('showwin');
else
figure(siga_DAT(12));
drawnow
lines=get(gcf,'userdata');
set(lines(2),'ydata',window);
W=fft(window,1024);
set(lines(3),'ydata',abs(W(1:512)));
figure(siga_DAT(12));
drawnow;
clear W lines;
end
clear show_window;
end
clear oldfig N t window;
elseif strcmp(action,'delshow'),
N=siga_DAT(3);
show_window=siga_DAT(14+N+N);
if ( siga_DAT(12)~=-1 & get(show_window,'value')==0),
show_window=findobj('Tag','window_window');
if isempty(show_window)
siga_DAT(12)=-1;
else
figure(show_window);
THEHAND=get(gcf,'userdata');
close;
figure(THEHAND(1));
siga_DAT(12)=-1;
clear THEHAND;
end
end
clear show_window N;
elseif strcmp(action,'setfreq'),
x = str2num(get(siga_DAT(10),'String'));
if isempty(x), % handle the non-numeric case
set(siga_DAT(10),'String',num2str(siga_DAT(1)));
else
siga_DAT(1)=x;
end;
elseif strcmp(action,'info'),
ttlStr = '连续信号傅立叶变换';
hlpStr1= str2mat(...
' ',...
' 我们看到，在演示窗口中，上方的图给出了一个 ',...
' 连续时间非周期信号的时域波形，用鼠标点击并 ',...
' 拖动信号上的任意一点就可改变该信号的幅值及 ',...
' 压缩或展宽信号的波形；也可通过用鼠标单击“带',...
' 限宽度”的文本框来压缩或展宽信号的波形。 ',...
' 下方的图给出了该连续信号经傅立叶变换后的频谱',...
' 密度的图形即信号的频域波形。 ' ,...
' 频域信号是通过连续信号的采样点经快速傅立叶 ',...
' 变换算法(FFT)而得到的。采样间隔由傅立叶变换 ',...
' 的点数决定。 ' ,...
' ' ,...
' 注意：下图中，频率范围是从0 到100 Hz。处于 ',...
' 负频域的 DFT 是正频域 DFT 的镜象，故在此只 ',...
' 画出了正频域的 DFT。抽样频率是 200Hz，这意 ',...
' 味着“奈奎斯特”频率是100Hz。在高于“奈奎斯 ',...
' 特”频率的频率点处的 DFT 与在低频（负频域） ',...
' 处的DFT相同。 ' ,...
' ');
hlpStr2=str2mat(...
' ',...
' 在上图中，点击并拖动信号波形上的任意一点， ',...
' 这样可以设置信号的带限范围与信号幅度。 ',...
' ' ,...
' 用鼠标单击“信号”文本框，就给出一个下拉式菜',...
' 单，从中选中信号即可改变输入的信号。 ' ,...
' ' ,...
' 菜单“窗口”让你选择一个窗口函数，在对信号 ',...
' 进行DFT 以前，这个窗口将与信号波形相乘—— ',...
' 即“开窗截取”。 ' ,...
' ' ,...
' 检查框“显示窗口函数”若被选中，则在显示信 ',...
' 号波形和其DFT的同时，也显示出窗口函数的形状 ' ,...
' ' ,...
' 左边的文本框给出了信号的带限值，可以通过用 ',...
' 鼠标点击此文本框，然后设置带限值的大小，再 ',...
' 回车来改变信号的形状。 ' );
helpwin(ttlStr, hlpStr1, hlpStr2);
end
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -