%=============设置系统参数==============% f1=380e6; %设置波形频率 f2=500e6; f3=300e6; Fs=512e6; %设置采样频率 L=1024; %数据长度 N=64; %数据位宽 Tp=60e-6; %预设的采样时间周期Tp为60us W=Tp*Fs; %=============产生输入信号==============% t=0:1/Fs:(1/Fs)*(L-1); y1=sin(2*pi*f1*t); y2=sin(2*pi*f2*t); y3=sin(2*pi*f3*t); y4=y1+y2+y3; y_n=round(y4*(2^(N-3)-1)); S0=fft(y4); S1=abs(S1); S2=(S1); S2=awgn(S2,10); %=================画图==================% f=0:Fs/W:Fs-Fs/W; figure(1); plot(f/1e6,64*log10(S2)); title('载频为384MHz有用信号带宽为64MHz的带通信号序列频谱'); xlabel('frequency(MHz)'); ylabel('Magnitude(dB)'); grid on %=============写入外部文件==============% fid=fopen('F:\FPGA\radar_channelize_process\xinhao.txt','w'); %把数据写入sin_data.txt文件中,如果没有就创建该文件 for k=1:length(y_n) B_s=dec2bin(y_n(k)+((y_n(k))<0)*2^N,N); for j=1:N if B_s(j)=='1' tb=1; else tb=0; end fprintf(fid,'%d',tb); end fprintf(fid,'\r\n'); end fprintf(fid,';'); fclose(fid); 代码改错
时间: 2023-08-06 13:19:23 浏览: 175
set /p= 命令使用详解
在这段代码中,变量 `S1` 没有被定义过,可能是 `S0` 的笔误。此外,变量 `S2` 在定义后仅仅做了一次加高斯噪声的操作,没有在频域上进行FFT变换,这是不合理的。以下是改正后的代码:
```matlab
%=============设置系统参数==============%
f1 = 380e6; %设置波形频率
f2 = 500e6;
f3 = 300e6;
Fs = 512e6; %设置采样频率
L = 1024; %数据长度
N = 64; %数据位宽
Tp = 60e-6; %预设的采样时间周期Tp为60us
W = Tp * Fs;
%=============产生输入信号==============%
t = 0:1/Fs:(1/Fs)*(L-1);
y1 = sin(2*pi*f1*t);
y2 = sin(2*pi*f2*t);
y3 = sin(2*pi*f3*t);
y4 = y1 + y2 + y3;
y_n = round(y4 * (2^(N-3)-1));
S0 = fft(y4);
S1 = abs(S0); % 将 S0 改为 S1
S2 = awgn(S1, 10); % 在频域上进行高斯噪声操作
%=================画图==================%
f = 0:Fs/W:Fs-Fs/W;
figure(1);
plot(f/1e6, 64*log10(S2));
title('载频为384MHz有用信号带宽为64MHz的带通信号序列频谱');
xlabel('frequency(MHz)');
ylabel('Magnitude(dB)');
grid on
%=============写入外部文件==============%
fid = fopen('F:\FPGA\radar_channelize_process\xinhao.txt', 'w');
% 把数据写入sin_data.txt文件中,如果没有就创建该文件
for k = 1:length(y_n)
B_s = dec2bin(y_n(k) + ((y_n(k))<0)*2^N, N);
for j = 1:N
if B_s(j) == '1'
tb = 1;
else
tb = 0;
end
fprintf(fid, '%d', tb);
end
fprintf(fid, '\r\n');
end
fprintf(fid, ';');
fclose(fid);
```
阅读全文