matlab波形逐渐增大,matlaB频率渐变的信号产生方法

close all;

% clear all;

SampFreq=1500;

t=-8:1/SampFreq:8;

N1=length(t);

% Sig=sin(2*pi*(5*t+4*t.^2));

Sig=sin(2*pi*(5*t+0.2*t.^2+0.2*t.^3+0.2*t.^4+0.2*t.^5))+2*sin(2*pi*2*t)+sin(2*pi*6*t);

Sig2=sin(2*pi*5*t)+2*sin(2*pi*4*t)+sin(2*pi*6*t);

sq=square(2*pi*5*t);

Sig=Sig+Sig2;

% plot(t,Sig2,'k',t,sq,'b');

plot(t,Sig,'k',t,sq,'b');

xlim([2,4])

[tf_b,tf_a]=tf(Hd);%得到传递函数;FDATool,Design a FIR Filter

tf_d1=filter(tf_b,tf_a,Sig);

tf_d2=filter(tf_b,tf_a,Sig2);

subplot(212); plot(t,tf_d1,'k-o',t,tf_d2,'b-o',t,sq,'b');

title('FIR时域');

xlabel('t/S');

ylabel('y');

grid on;

% plot(t,Sig,'k.',t,Sig2,'g.');

hold on;

plot(Sig2);

% Sig=Sig+sin(2*pi*(45*t-4*t.^2));

y=abs(fft(Sig));

y1=y*2/N1;

tf=SampFreq*(0:N1-1)/N1

subplot(211);

plot(t,Sig,'k-.');

subplot(212);

plot(tf,y1)