www.cn-17.net 发表于 2012-5-28 13:59

要什么样的？

我这里有一个调制部分的  输出了调制信号波形和其频谱 你看能在此基础上加上解调信号的程序吗  还要求输出误码率  (4ppm就是单脉冲位置的调制与解调  及一帧有四个时隙调制出的波形是一帧李只有一个时隙有脉冲(显示为1))  真的非常感谢

%######################################%

% 无线光通信系统中光调制编码

% --------*********-------*********-----

% Function:

% Programme  name:

% Programmer:

% The guidance of teacher: JiaKejun

% Date: 04/05/2012

%**************************************%

% Define variables:

%

%

%

%

%

%

%

%

%

%

%

%

%

%

%*******************************%

%--清理内存--%

clc;

clear all;

close all;

%*************** 编码程序  ***************%

%------------输入维数与编码序列-----------%

L=input('请输入编码维数L的值:');           % 例如： 3

% P=input('请以行矩阵形式输入OOK序列码:'); % 例如：[1 0 0 0 1 1 1 0 1]

clc;

disp('编码维数L的值为:');

disp(L);

P=randn(1,8)>0;

%-----------输出OOK序列------------%

c=1;

% u=numel(P);      %求输入元素的个数

M=floor(8/L);     %求编码的组数

n=M*2^L;          %求OOK转换为PPM方式时的位码数

A(L,M)=zeros;    %A设初值

for i=1:M        %循环求OOK转换为PPM方式的编码

for j=1:L

if c<=M*L

A(j,i)=P(1,c);

c=c+1;

end

end

end

disp('传输的OOK序列为:A=');

disp(A);

t1=0:0.01:M*L;     %设置抽样频率

%----------输出OOK编码图像------------%

figure(1)

for i=1:M

for j=1:L

if A(j,i)==1

date_modulated1=rectpuls(t1-((i-1)*L+j-0.5));

subplot(2,1,1),title('OOK序列'),xlabel('TIME'),v=line(t1,date_modulated1);

grid on,set(gca,'xtick',(1:M*L),'ytick',[0.5,1]),hold on;

set(v,'Color',[0 0 1],'LineWidth',2),axis([0,M*L,0,1.1]);

end

end

end

%----------输出OOK编码图像中的分组分隔线------------%

for i=1:M

t1=L*i;

y1=0:0.2:1;

v1=plot(t1,y1,'r*');

end

for i=1:M

for j=1:L

if A(j,i)==1

y=fft(date_modulated1);

plot(abs(y)),grid on;

subplot(2,1,2),title('OOK频谱');

axis([0,200,0,60]);

end

end

end

figure(2)

signal_length1=length(date_modulated1);

noise1=1*randn(1,signal_length1);

signal_addnoise1=date_modulated1+noise1;

subplot(3,1,1);

plot(signal_addnoise1),grid;

title('加噪后的OOK序列');

axis([0,800,-5,5]);

subplot(3,1,2);

y9= xcorr(signal_addnoise1);

plot(abs(y9)),grid;

title('加噪后的OOK序列的自相关函数');

axis([0,1600,0,1000]);

subplot(3,1,3);

y7=fft(y9);

plot(abs(y7)),grid;

title('加噪后的OOK调制信号频谱');

%---------输出PPM调制序列---------------%

C(1,L)=zeros;      %对C矩阵赋初值

B(2^L,M)=zeros;    %对B矩阵赋初值

b(1,M)=zeros;

for i=1:M

for j=1:L

C(1,j)=A(j,i)*2^(L-j);

b(1,i)=sum(C);

end

end

for i=1:M

for j=1:2^L

if j==b(1,i)+1

B(j,i)=1;

end

end

end

disp('传输的PPM序列为:B=');

disp(B);

t2=0:0.01:n;     %设置抽样频率

%------------输出PPM调制图像-------------%

figure(3)

for i=1:M

for j=1:2^L

if B(j,i)==1

date_modulated=rectpuls(t2-((i-1)*2^L+b(1,i)+0.5));

subplot(2,1,1),title('PPM调制序列'),xlabel('TIME'),v=line(t2,date_modulated);

grid on,set(gca,'xtick',(1:n),'ytick',[0.5,1]),hold on;

set(v,'Color',[0 1 1],'LineWidth',2),axis([0,n,0,1.1]);

end

end

end

%------------输出PPM调制图中的分组分隔线-------------%

for i=1:M

t1=2^L*i;

y1=0:0.2:1;

v1=plot(t1,y1,'r*');

end

for i=1:M

for j=1:2^L

if B(j,i)==1

y1=fft(date_modulated);

plot(abs(y1)),grid on;

subplot(2,1,2),title('PPM频谱');

axis([0,200,0,80]);

end

end

end

figure(4)

signal_length=length(date_modulated);

noise=1*randn(1,signal_length);

signal_addnoise=date_modulated+noise;

subplot(3,1,1);

plot(signal_addnoise),grid;

title('加噪后的PPM调制信号');

axis([0,1000,-4,4]);

subplot(3,1,2);

y5= xcorr(signal_addnoise);

plot(abs(y5)),grid;

title('加噪后的PPM调制信号的自相关函数');

subplot(3,1,3);

y8=fft(y5);

plot(abs(y8)),grid;

title('加噪后的PPM调制信号频谱');