ANGLE

最新推荐文章于 2019-09-25 08:53:37 发布

不断的进步

最新推荐文章于 2019-09-25 08:53:37 发布

阅读量239

点赞数

CC 4.0 BY-SA版权

分类专栏：基于matable算法

本文链接：https://blog.youkuaiyun.com/qq_43359781/article/details/90613667

基于matable算法专栏收录该内容

2 篇文章

订阅专栏

clf;
%我的
%fm=100;f0=2.4e10;B=2.4e8;FR=320R;FD=160v;
Fs = 10240;            % Sampling frequency                    
Ts = 1/Fs;             % Sampling period       
N = 512;             % Length of signal
n=0:N-1;
%t1 = (0:1:N-1)*Ts;        % Time vector
t1=n/Fs;
%y = wgn(m, n, p);
% 遮挡临近目标
S = sqrt(2)*sqrt(5)*sin(2*pi*5000*t1+pi/4) ;
NOISE=wgn(1, 512, 5);
X = S+NOISE;
sigPower = sum(abs(S).^2)/length(S);           %求出信号功率
noisePower=sum(abs(X-S).^2)/length(X-S); %求出噪声功率
SNR=10*log10(sigPower/noisePower) ;

Y = fft(X); %求得Fourier变换后的振幅
P2=angle(Y)/pi*180;
P2=P2(1:256);
mag1=abs(2*Y/N);
P1 = mag1(1:N/2);
%P1(2:end-1) = 2*P1(2:end-1);
f1=(1:N/2)*Fs/N;    %频率序列
%figure(1);plot(1000*t1,X);title('');xlabel('t (milliseconds)');ylabel('X(t)');
figure(2);plot(f1,P1,'r');  xlabel('频率/Hz');ylabel('dbv');title('： FFT','color','r');
figure(3);stem(P2);