matlab程序求解释……越详细越好……程序如下……若解释的令人满意,分大大地有赏……先上20分……

function y=pipei(x)%构造发送信号波形t1=0:100;--自变量,间隔为1,等价于t1=0:1:100s01=-2+0*t1;s11=2+0*t1;--实际上产生了t1上的两个常量信号 -2和2t2=100:200;--第二个区间段s02=2+0*t2;s12=-2+0*t2;--第二个区间t2上的两个常量信号t3=200:1:400;s03=0*t3;s13=0*t3;--第三个区间t3上的两个零信号t=[t1,t2,t3];--下面这三行将自变量区间和函数空间合并。s0=[s01,s02,s03];s1=[s11,s12,s13];%匹配滤波器系统函数h0=s1;h1=s0;%码元周期 Tb=200; for i=1:length(x)--x是函数的传入参数即pipei(x)中的x,length求其长度if(x(i)==1)for j=1:Tbd((i-1)* Tb+j)=1;--如果x的第i个元素为1,执行此循环--此时每个x元素对应产生一个Tb（200）长的1序列end elsefor j=1:Tb d((i-1)* Tb+j)=0;--否则产生一个Tb（200）长的零序列end end;end;--这个循环结束后,d就是一个200*length(x)长度的数列,每一个x元素对应d中200个值,1对应1,其它数对应0。%对发送二进制信号进行波形映射--后面的太长,具体功能我就不解释了,而且我对编码也不大懂,只给出函数代码表面的含义k=1;for g=1:length(x)--for循环if(x(g)==1)y=s1;elsey=s0;end;%叠加高斯白噪声,信噪比设为10dB ynr=awgn(y,10,'measured');--产生高斯噪声%产生输出信号r0=ifft(fft(ynr,512)。*fft(h0,512),512);--两个信号分别fft之后相乘,对乘积再进行ifft（逆傅里叶变换）,时域上看就是卷积。512指进行fft或ifft的点数。r1=ifft(fft(ynr,512)。*fft(h1,512),512); --同上tt=0:511;figure(k);--绘图subplot(211);--将一个图像分成上下两块,这里绘制第一块plot(t,y,'k',t,ynr,'g');--绘制两条曲线,一个是黑色black,一个是绿色greenaxis([0,400,-3,3]);--设置坐标系的范围：x是从0到400,y轴从-3到3legend('发送波形','叠加高斯白噪声后波形');--给两条曲线添加标记subplot(212);--绘制第二块plot(tt,r0,'b',tt,r1,'r');--一条是蓝色blue,一条是红色redaxis([0,400,-900,900]);--设置坐标轴范围gtext('\leftarrow r0=conv(y,h0)');--在图像上给出文本标记gtext('\leftarrow r1=conv(y,h1)');if r0(Tb)>r1(Tb)m(k)=0;else m(k)=1;end;k=k+1;end;for i=1:length(m)--跟上面的大同小异if(m(i)==1)for j=1:Tbn((i-1)* Tb+j)=1;end elsefor j=1:Tb n((i-1)* Tb+j)=0;end end;end;figure(k);--绘图,与上面的大同小异。subplot(211);plot(d);axis([0,length(d),-0。2,1。2]);title('发送的二进制信号');subplot(212);plot(n);axis([0,length(n),-0。2,1。2]);title('接收判定后的二进制信号');