求教大哥们～～帮小弟解答下matlab通信仿真关于GUI的问题

cyh5000 · 发表于 2006-6-12 20:41

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

<v:shapetype>
关于最佳基带系统的误码率分析模型
<v:shapetype><v:stroke joinstyle="miter"></v:stroke><v:formulas><v:f eqn="if lineDrawn pixelLineWidth 0"></v:f><v:f eqn="sum @0 1 0"></v:f><v:f eqn="sum 0 0 @1"></v:f><v:f eqn="prod @2 1 2"></v:f><v:f eqn="prod @3 21600 pixelWidth"></v:f><v:f eqn="prod @3 21600 pixelHeight"></v:f><v:f eqn="sum @0 0 1"></v:f><v:f eqn="prod @6 1 2"></v:f><v:f eqn="prod @7 21600 pixelWidth"></v:f><v:f eqn="sum @8 21600 0"></v:f><v:f eqn="prod @7 21600 pixelHeight"></v:f><v:f eqn="sum @10 21600 0"></v:f></v:formulas><v:path connecttype="rect" gradientshapeok="t" extrusionok="f"></v:path><lock aspectratio="t" v:ext="edit"></lock></v:shapetype>
(一)，用到的函数
1、f2t函数
 function x=f2t(X)
global dt df t f T N
%x=f2t(X)
%x为时域的取样值矢量
%X为x的傅氏变换
%X与x长度相同并为2的整幂
%本函数需要一个全局变量dt(时域取样间隔)
X=[X(N/2+1:N),X(1:N/2)];
x=ifft(X)/dt;
end
2、t<st1:chmetcnv w:st="on" TCSC="0" NumberType="1" Negative="False" HasSpace="False" SourceValue="2" UnitName="F">2f</st1:chmetcnv>函数
function X=t<st1:chmetcnv w:st="on" TCSC="0" NumberType="1" Negative="False" HasSpace="False" SourceValue="2" UnitName="F">2f</st1:chmetcnv>(x)
global dt df N t f T
%X=t<st1:chmetcnv w:st="on" TCSC="0" NumberType="1" Negative="False" HasSpace="False" SourceValue="2" UnitName="F">2f</st1:chmetcnv>(x)
%x为时0域的取样值矢量
%X为x的傅氏变换
%X与x长度相同,并为2的整幂。
%本函数需要一个全局变量dt(时域取样间隔)
H=fft(x);
X=[H(N/2+1:N),H(1:N/2)]*dt;
end
(二),主程序
%本程序中时间单位是us
%频率单位为MHZ
%码速率单位是Mb/s
global dt df t f N T
close all
clear Eb_N0 Pe
 N=2^13; %采样点数
 L=8; %每个码元的采样点数
 M=N/L; %码元数
 Rb=2; %码速率是2MB/s
 Ts=1/Rb; %码元间隔
 dt=Ts/L; %时域采样间隔
 df=1/(N*dt); %频域采样间隔
 T=N*dt; %截短时间
 Bs=N*df/2 %系统带宽
 alpha=1;
 t=linspace(-T/2,T/2,N);
 f=linspace(-Bs,Bs,N)+eps;%eps为了避免f=0;
 figure(1)
 set(1,'position',[10,50,500,500])
 hr1=sin(pi.*t/Ts)./(pi.*t/Ts);
 hr2=cos(alpha*pi.*t/Ts)./(1-(2*alpha.*t/Ts).^2);
 hr=hr1.*hr2; %升余弦脉冲波形
 HR=abs(t2f(hr)); %取模是为了忽略时延
 GT=sqrt(HR);
 GR=GT;
 for loop1=1:10
 Eb_N0(loop1)=(loop1-1)
 eb_n0(loop1)=10^(Eb_N0(loop1)/10);
 Eb=1;
 n0=Eb/eb_n0(loop1);%信道的噪声谱密度
 sita=n0*Bs; %信道中噪声功率
 n_err=0; %误码计数
 for loop2=1:5
 a=sign(randn(1,M));
 imp=zeros(1,N);
 imp(L/2:L:N)=a/dt; %产生冲激序列
 IMP=t2f(imp);
 n_ch=sqrt(sita)*randn(size(t));%信道噪声
 nr=real(f2t(t2f(n_ch).*GR)); %输出噪声
 sr=real(f2t(IMP.*HR))+nr; 
 y=sr(L/2:L:N); 
 aa=sign(y); %判决
 n_err=n_err+length(find(aa~=a));
end
pe(loop1)=n_err/(M*loop2)
figure(1)
 semilogy(Eb_N0,pe,'r');
 eb_n0=10.^(Eb_N0/10);
 hold on
 semilogy(Eb_N0,0.5*erfc(sqrt(eb_n0)));
 axis([0,9,1e-4,1])
title('最佳接收基带接收系统的误码率分析')
 Xlabel('Eb/N0')
 Ylabel('Pe')
End 请教大哥们～怎么将上面的程序通过GUI制作成一个可以通过改变Eb的数值从而改变Pe与Eb/N0图形的分析模型 小弟感激不尽～～请各位大哥们指教指教！！！谢谢 </v:shapetype>

cyh5000 · 发表于 2006-6-13 12:53

或者请高手帮我把上面主程序改成function定义的关于Eb的函数文件。。。 像上面的f2t、t2f函数那样子d ~~ 期待各位大哥指导！先谢谢了

账号		自动登录	找回密码
密码			我要加入