%該程序產生三組偽隨機數，模擬隨機編碼過程

close?all;
clear;
clc;
t=0:5:1000;??????%時間變量為1s，步長5ns


x1=rand（1length（t））;???????????????????????%生成一組1*length（t）維的偽隨機數
x1=round（x1）;???????????????????????????????%將x1數組中的元素轉換到與其最近的整數??????????
index1=find（x1==0）;?????????????????????????%找出x1數組中為0的元素
x1（index1）=?-ones（1length（index1））;????????%返回一個index1長度的一維數組


x2=rand（1length（t））;
x2=round（x2）;
index1=find（x2==0）;
x2（index1）=?-ones（1length（index1））;



x3=rand（1length（t））;
x3=round（x3）;
index1=find（x3==0）;
x3（index1）=?-ones（1length（index1））;

y1=zeros（11）;y2=zeros（11）;y3=zeros（11）;?????%生成三個1*1的零矩陣，存放偽隨機值
t1=zeros（11）;t2=zeros（11）;t3=zeros（11）;?????%生成三個1*1的零矩陣，存放時間值

%%將生成的偽隨機值存放到對應的矩陣中
y1=x1（1）;y2=x2（1）;y3=x3（1）;
t1=t（1）;t2=t（1）;t3=t（1）;
k=2;m=2;n=2;

for?i=2:length（t）
????if（x1（i）==x1（i-1））
????????y1（k）=x1（i）;
????????t1（k）=t（i）;
????????k=k+1;
????????
????else
????????y1（k）=x1（i-1）;
????????y1（k+1）=x1（i）;
????????t1（k）=t（i）;
????????t1（k+1）=t（i）+0.01;????%如果該時刻階躍變化，則右移0.01來記錄
????????k=k+2;
????end
?????
????if（x2（i）==x2（i-1））
????????y2（m）=x2（i）;
????????t2（m）=t（i）;
????????m=m+1;
????????
????else
????????y2（m）=x2（i-1）;
????????y2（m+1）=x2（i）;
????????t2（m）=t（i）;
????????t2（m+1）=t（i）+0.01;??
????????m=m+2;
????end
????
????if（x3（i）==x3（i-1））
????????y3（n）=x3（i）;
????????t3（n）=t（i）;
????????n=n+1;
????????
????else
????????y3（n）=x3（i-1）;
????????y3（n+1）=x3（i）;
????????t3（n）=t（i）;
????????t3（n+1）=t（i）+0.01;???
????????n=n+2;
????end
end


%%%畫圖說明產生的三種偽隨機碼
subplot（311）;
plot（t1y1‘y‘）;
grid?on;
axis（[-11001-1.51.5]）;

subplot（312）;
plot（t2y2‘y‘）;
grid?on;
axis（[-11001-1.51.5]）;

subplot（313）;
plot（t3y3‘y‘）;
grid?on;
axis（[-11001-1.51.5]）;
whitebg（‘black‘）;

%將三組偽隨機值存放到新矩陣并傳送到Workspace中，方便以后調用數據
simin1=[t1;y1]‘;
simin2=[t2;y2]‘;
simin3=[t3;y3]‘;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2069??2013-07-21?16:41??GenerateCode.m

?????文件???????1583??2013-07-22?13:12??GPS_signal.m

-----------?---------??----------?-----??----

?????????????????3652????????????????????2