Question

各位matlab大神给我解读一下这些程序，尽量具体点，谢啦！！！有奖励的

fc=2*10^9;
v=110;
c=300*10^6;
fm=fc*(v*10^3/3600)/c;
N=128;
gap=2*fm/(N-1);
T=1/gap;
sf0=1.5/(pi*fm);
for n=（1:(N-2)/2）
sf(n) = 1.5/(pi*fm*sqrt(1-(n*gap/fm)^2));
end
SEf = [fliplr(sf),sf0,sf];
figure(1);
plot(SEf)，title('多普勒功率谱')， xlabel('f')， ylabel('size')，
grid;
x1=randn(1,N);

x2=randn(1,N);
x=x1+1i*x2;
pha=2*pi*rand(1,N);
i= linspace(1,N,N);
sigma=sqrt(var(x)/2);
syms t;
tt=linspace(0,0.35,N);
ci=sigma*(2/N).^0.5;
fi=fm.*sin(pi/(2*N)*(i-1/2));
ph=2.*pi.*fi*t+pha;
Tc1=ci.*cos(ph);
Ts1=ci.*sin(ph);
Tc2=sum(Tc1);
Ts2=sum(Ts1);
Tc=subs(Tc2,t,tt);
Ts=subs(Ts2,t,tt);
st=Tc.*cos(2*pi.*tt)-Ts.*sin(2*pi.*tt);
rt = sqrt((real(Tc).^2 + real(Ts).^2));
% rt=sqrt(Tc.^2+Ts.^2)

rt_db=20*log10(rt);
st_db=20*log10(st);
figure(2);
plot(tt,rt_db)，axis([0 20，0 100]) ,title('瑞利衰落信道')，xlabel('t')，ylabel('dB')，grid
figure(3)
n=0:0.1:10;
r=sqrt(sigma*(x1.^2+x2.^2));
h=hist(r,n);
fr_approx=h/(0.1*sum(h));
pijun=sum(r)/N;
junfanghe=(r-pijun).^2;
junfang=sum(junfanghe)/N;
u=0;
cdf=raylcdf(n,sigma);
subplot(3,1,1);
plot(n,cdf)，title('CDF');
pdf=raylpdf(n,sigma);

subplot(3,1,2);
plot(n,pdf)，title('PDF');
hold on;
plot(n,fr_approx,'ko')；
axis([0 8，0 1]) ；
wucha=fr_approx-pdf; subplot(3,1,3);
plot(n,wucha)，title('误差');
R=raylrnd(sigma,1,1000);

E=mean(R);
D=cov(R);
figure(4)
plot(Tc,'.'),title('时域高斯信号(Tc)'),xlabel('N'),ylabel('正交信号'), grid;
%时域高斯信号Ts
figure(5)
plot(Ts,'.')，title('时域高斯信号(Ts)')， xlabel('N')， ylabel('正交信号')，grid;
我希望各位大神能给我解释一下，每行程序代码是什么意思？qq：929804034

Answer 1

不知到你想问什么，是程序代码看不懂？还是关于matlab瑞利信道仿真相关

追问

不好意思，我没说清楚，我希望有人能给我解释一下，每行程序代码是什么意思？

追答

每一行都不懂吗？这个只要学一下matlab语法即可啊

追问

基本的懂点，刚开始学，大神帮帮忙吧！！！

追答

我先解释一段、其余的应该都是差不多~
遇到不懂的 可以 在matlab中 输入 help 函数名 需求帮助
========================================================================
fc=2*10^9; // fc等于2乘以10的9次方
v=110; //V等于110
c=300*10^6; //C等于300乘以10的6次方
fm=fc*(v*10^3/3600)/c; //fm等于fc×(v×10^3/3600)÷c
N=128; //N等于128
gap=2*fm/(N-1); //gap等于2×fm÷(N-1);
T=1/gap; //T=1÷gap
sf0=1.5/(pi*fm);//sf0=1.5÷(pi×fm);其中 PI 应该是常量3.1415926...

for n=（1:(N-2)/2） //此处是一个循环 n从1 到 (N-2)/2
sf(n) = 1.5/(pi*fm*sqrt(1-(n*gap/fm)^2)); //sf应该是数组，不断的将右边值赋值给 sf(n)
end //循环结束
SEf = [fliplr(sf),sf0,sf]; //定义sef fliplr 表示翻转、即将sf数组中数值翻转 ，//假设 sf中数值为1,2,3 那么 sef=[3，2，1，sf0，1，2，3]
figure(1);
plot(SEf)，title('多普勒功率谱')， xlabel('f')， ylabel('size')，grid; //此处表示绘图 ，数据位 SEF，标题 多普勒功率谱 ，横坐标 名为 f，纵坐标名为 size gird 表示图表中显示网格

x1=randn(1,N); //取 1到N中随机数赋值给 x1

x2=randn(1,N); //取 1到N中随机数赋值给 x2

x=x1+1i*x2; //1i 复数单位

pha=2*pi*rand(1,N);
i= linspace(1,N,N); //linspace此函数用来生成1到n之间的等间距的数组 即 1,2，3...n

sigma=sqrt(var(x)/2);//sqrt 平方根
syms t; //相当于 t=syms(‘t’)

追问

哇哇！厉害呀！好人啊！再给翻译翻译吧，我也知道，help，只是本人英语太渣，看不懂啊，谢啦！！！！！