function [y_original]=QPSK()%写错名字啦!!
clc;
clear;
delta=2*pi/32;
x=0:0.01*pi:2*pi-0.01*pi;%自变量
original_phase=0:delta:2*pi-delta;%几种可能出现的相位情况
size_phase = length(original_phase);
% figure(1);
% plot(x,y_sin,'r');
% hold on
% plot(x,y_cos,'g');
size_x = length(x);%抽样点个数
y_original=zeros(1,size_phase*(size_x));%信号a
for i=1:size_phase
for j=1:size_x
y_original(j+(i-1)*size_x)=sin(x(j)+original_phase(i));%原始信号
end;
end;
x1=0:0.01*pi:size_phase*2*pi-0.01*pi;
y_sin=sin(x1);%正弦信号
y_cos=cos(x1);%余弦信号
length(x)
length(y_original((1:size_x)+800))
find(y_original==0)
figure(1);
% plot(x1,y_original,'*b',x1,y_sin,'.r');
plot(x1,y_original,'*b');
figure(2);
y_y_sin=y_original.*y_sin;
for i=1:size_phase
plot(x1,y_y_sin,'*');
%axis在plot之后有效限制坐标,放在前面没有作用
axis( [0 i*2*pi -1 1] );%限制横纵坐标时必须有[]
% pause;
end;
figure(3);
y_y_cos=y_original.*y_cos;
for i=1:size_phase
plot(x1,y_y_cos,'.');
%axis在plot之后有效限制坐标,放在前面没有作用
axis( [0 i*2*pi -1 1] );%限制横纵坐标时必须有[]
pause;
end;