请matlab大神帮我看看我的程序吧

....

figure(2)

B2 = repmat(B1,1,1000)

mesh(x,y,B2)

(单缝)

4.

[x,y]=meshgrid(xs,xs)%产生网格

.....

figure(2)

B2 = repmat(Br,500,1)

mesh(x,y,B2)%三维强度图像

clear

close all

Fd=1%消息序列的采样速率

Fs=3*Fd%以调信号的采样速率

M=4%4比特调制，qpsk

SNR=0:12%去信噪比从0到12

SNR1=0:0.1:12%步进为0.1

for n=1:length(SNR)%对每一snr计算一次

Eb_N0=10^(SNR(n)/10)%换算10进制

sgma=sqrt(1/(4*Eb_N0))%噪声信号的平方差开根号

x=randint(10000,1,M)%产生10000个输入符号

y=dmodce(x,Fd,Fs,'psk',M)%调制

ynoise=y+sqrt(Fs/Fd)*sgma*(randn(length(y),1)+1i*randn(length(y),1))%经过信道加噪声

z=ddemodce(ynoise,Fd,Fs,'psk',M)%解调

[numbers,pm(n)]=symerr(x,z)&求误符号率

end

for m=1:length(SNR1)%此循环用于画出理想的误符号率曲线

Eb_N0=10^(SNR1(m)/10)

pm1(m)=erfc(sqrt(log2(M)*Eb_N0)*sin(pi/M))

end

semilogy(SNR,pm,'-*',SNR1,pm1)%纵坐标取对数

xlabel('信噪比SNR/db')

ylabel('误码率(pe)')

legend('仿真误码率曲线','理论误码率曲线')

这个net是个结构体类型,其中结构体中的变量LW又是个元胞数组。

net.LW{2,1}=net.LW{2,1}*0.01这句话是把结构体net中的lw元胞数组的第2行第一个元素变为原来的0.01倍

下面一句意思类似。

matlab元胞数组：http://www.3800hk.com/Article/jxdz/matlab/matlabjc/2005-08-06/Article_35366.html

matlab结构体：http://hi.baidu.com/%C6%BF%D6%D0%B5%C4%C5%AE%CE%D7/blog/item/91939438af235ef43a87ce08.html

---