fsk调制与解调源代码

lf251190685

新建 文本文档 (2).txt (1.74 KB, 下载次数: 5)

2009-12-12 13:56 上传

点击文件名下载附件

cowry

这样的帖子下次请直接粘贴内容，谢谢合作


clear all;

x=randn(10,1)>0;

fl=5000;fh=8000;fs=100000;

ts=1/200;%码元速率200波特

tt=(0:1/fs:ts);

t=[tt;tt+ts;tt+2*ts;tt+3*ts;tt+4*ts;tt+5*ts;tt+6*ts;tt+7*ts;tt+8*ts;tt+9*ts];

y=zeros(10,length(tt));

i=1;

%对该输入信号FSK调制

while i<=10

y(i,:)=x(i)*sin(2*pi*fh*t(i,:))+~x(i)*sin(2*pi*fl*t(i,:));

i=i+1;

end

t=reshape(t',length(tt)*10,1);

y=reshape(y',length(tt)*10,1);

subplot(211);plot(t,y);

title('FSK信号的时域图形');

%该输入信号的频域图形

n=length(y);

r=fft(y)/n;r=fftshift(r);

f=linspace(-fs/2,fs/2,n);

subplot(212);

plot(f,abs(r));

set(gca,'XTick',-fs/2:5000:fs/2);

title('FSK信号的频谱图');





clear all

%二进制随机序列

x=[0 1 1 0 0 1 0 1 0 1];%定义一个二进制10元素数组，作为FSK信号输入数据。
figure(1);
stem(x,'.');
title('二进制随机序列');xlabel('时间');ylabel('幅度');

%FSK信号的调制

f0=1000;f1=2000;fs=8000;ts=1/125;%0对应频率f0，1对应频率f1，采样频率fs，码元速率125波特。
tt=(0:1/fs:ts);
t=[tt;tt+ts;tt+2*ts;tt+3*ts;tt+4*ts;tt+5*ts;tt+6*ts;tt+7*ts;tt+8*ts;tt+9*ts;];
y=zeros(10,length(tt));
i=1;
%开始调制
while i<=10
y(i,:)=x(i)*cos(2*pi*f1*t(i,:))+~x(i)*cos(2*pi*f0*t(i,:));
i=i+1;
end
t=reshape(t',length(tt)*10,1);%阵列重新排列
y=reshape(y',length(tt)*10,1);
figure(2);
plot(t,y);
title('FSK信号时域波形');xlabel('时间');ylabel('幅度');

tt=(0:1/fs:ts);
t=[tt;tt+ts;tt+2*ts;tt+3*ts;tt+4*ts;tt+5*ts;tt+6*ts;tt+7*ts;tt+8*ts;tt+9*ts;];
y=zeros(10,length(tt));
i=1;
%开始调制
while i<=10
y(i,:)=x(i)*cos(2*pi*f1*t(i,:))+~x(i)*cos(2*pi*f0*t(i,:));
i=i+1;
end
t=reshape(t',length(tt)*10,1);%阵列重新排列
y=reshape(y',length(tt)*10,1);
figure(2);
plot(t,y);
title('FSK信号时域波形');xlabel('时间');ylabel('幅度');

gjlin

一直找这个啊，就要这个