matlab语言 麻烦把每一行注释一下 本人小白 望大虾指教 采纳后还有重谢

matlab语言 麻烦把每一行注释一下 本人小白 望大虾指教 采纳后还有重谢

你也太贪心了吧，程序是别人的吧。大概帮你l捋捋。
% 调相法
clear all;
t=[-1:0.01:7-0.01]; %变量赋值
tt=length(t); %循环的长度
x1=ones(1,800); %定义1*1800全1矩阵x1

for i=1:tt
if (t(i)>=-1 & t(i)<=1) | (t(i)>=5& t(i)<=7); %如果-1= x1(i)=1; %矩阵x1的对应元素=1
else
x1(i)=-1; %否则，矩阵x1的对应元素=-1
end
end

t1=[0:0.01:8-0.01];
t2=0:0.01:7-0.01;
t3=-1:0.01:7.1-0.01;
t4=0:0.01:8.1-0.01; %定义4个变量
tt1=length(t1); %循环的长度
x2=ones(1,800); %定义1*1800全1矩阵x2

for i=1:tt1
if (t1(i)>=0 & t1(i)<=2) | (t1(i)>=4& t1(i)<=8); %如果0= x2(i)=1; %矩阵x2的对应元素=1
else
x2(i)=-1; %否则，矩阵x2的对应元素=-1
end
end

f=0:0.1:1; %定义频率变量f
xrc=0.5+0.5*cos(pi*f); %给函数xrc赋值
y1=conv(x1,xrc)/5.5; %计算x1和xrc的卷积
y2=conv(x2,xrc)/5.5; %计算x2和xrc的卷积

n0=randn(size(t2)); %产生与t2一样大的随机矩阵n0
f1=1;
i=x1.*cos(2*pi*f1*t); %计算i分量
q=x2.*sin(2*pi*f1*t1); %计算q分量
I=i(101:800);
Q=q(1:700); %截取部分数据作为I、Q
QPSK=sqrt(1/2).*I+sqrt(1/2).*Q; %计算正交相移键控
QPSK_n=(sqrt(1/2).*I+sqrt(1/2).*Q)+n0; %叠加前面的随机数n0

n1=randn(size(t2)); %产生与t2一样大的随机矩阵n1
i_rc=y1.*cos(2*pi*f1*t3); %计算i分量
q_rc=y2.*sin(2*pi*f1*t4); %计算q分量
I_rc=i_rc(101:800);
Q_rc=q_rc(1:700);
QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);
QPSK_rc_n1=QPSK_rc+n1; %叠加前面的随机数n1

figure(1)
subplot(4,1,1); %图形1有4个子图
plot(t3,i_rc); %画变量t3的i分量序列
axis([-1 8 -1 1]);
ylabel('a序列');
subplot(4,1,2); %子图2
plot(t4,q_rc); %画变量t4的q分量序列
axis([-1 8 -1 1]);
ylabel('b序列');
subplot(4,1,3); %子图3
plot(t2,QPSK_rc); %画变量t2的QPSK
axis([-1 8 -1 1]);
ylabel('合成序列');
subplot(4,1,4); %子图4
plot(t2,QPSK_rc_n1); %画变量t2的A加入噪声的QPSK
axis([-1 8 -1 1]);
ylabel('加入噪声');

后悔