clear all
clear all
f0 = 0.2e6; % 载波频率
fs = 2e6; % 采样频率
b = 400; % 码元速率
l = 0.1; % 波形的时间长度
T = fs/b; % 基带码的宽度
N = l*b; % 产生所需基带码的长度
L = N*T; % 所需数据的总长度
% xx = round(rand(1,N)); % 产生随机的基带码流
% figure(1)
% stem(xx)
for i = 1:N;
xx(i) = mod(i,2);
end
fp=fopen('code010.dat','wb'); % 写文件
for i=1:length(xx)
fprintf(fp,'%5.5d\n',xx(i));
end
for i = 1:N; % 将基带码流展宽
yy((1+(i-1)*T):i*T) = xx(i);
end
% figure(2)
% plot(yy)
theta = 2*pi*f0; % 载波角频率
t = (0:L-1)/fs;
zz = sin(theta*t);
zz = yy.*zz;
zz = zz+1;
zz = round(zz/max(zz)*(2^14-1));
figure(3)
plot(zz)
% 展宽插值
% for ii=1:length(zz)
% zzz((1+(ii-1)*4):ii*4) = zz(ii);
% end
fd=fopen('ask010.dat','wb'); % 写文件
for i=1:length(zz)
fprintf(fd,'%5.5d\n',zz(i));
end
fclose('all');
msgbox('ask.dat,write ok!');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 产生噪声，控制信噪比s/n
%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% nn = randn(1,L); % 产生正态分布的随机数
% nn = nn-min(nn); % 将噪声数据变成正值
% nn = round(nn/max(nn)*(2^14-1));
% figure(4)
% plot(nn);
%
% fd=fopen('noise.dat','wb'); % 写文件
% for i=1:length(nn)
% fprintf(fd,'%5.5d\n',nn(i));
% end
% fclose('all');
% msgbox('noise.dat,write ok!');
评论0