SCILABでの波形生成＆送受信

１）SCILAB関連

S1)　１周期１６点のCOS波形の生成

-->Fs=1/16

-->n=0:Fs:10

-->x=cos(2*%pi*n)
-->plot2d(n,x)
Fsは計算する点のSTEP

nに０から１０STEP=Fsの値を設定（横ベクトル）

%pi　は円周率　nと同一サイズのxベクトルが計算される

横軸をｎ、縦軸をｘで２次元プロット

S2)　位相を０、πの２種類で、位相変調PSK

-->Fs=1/16

-->n=0:Fs:8-Fs

-->zero=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]

-->one=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]

-->phi=%pi*[zero one zero one one one zero zero]

-->subplot(2,1,1)

-->plot2d(n,phi)

-->subplot(2,1,2)

-->x=cos(2*2*%pi*n+phi)
-->plot2d(n,x) 　
　

１シンボル１６点なので、１６回コピー

　

８シンボル分の位相情報を生成

たて２個のプロットの上側１番に

phiを２次元プロットする。

たて２個のプロットの２番目を指定

１シンボル＝１秒として、２Hzの波に対して位相変調している

結果をプロット

S3)　S2)と同様の変調を複素振幅Xで実現する場合

-->Fs=1/16

-->n=0:Fs:8-Fs

-->zero=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]

-->one=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]

-->phi=%pi*[zero one zero one one one zero zero]

-->X=exp(%i*phi)

-->x=X　.*　exp(%i*2*2*%pi*n)
-->param3d(real(x), imag(x), n) 　
　

　

　

ここまでは　S2)と同じ

複素振幅Xの計算、値としては’１’か’－１’である。

exp(%i*2*2*%pi*n) の各要素とXの各要素を乗算している。（ドット＊演算）

３次元曲線の描画

-->scf();

-->subplot(2,1,1)

-->plot2d(n,real(x))

-->subplot(2,1,2)
-->plot2d(n,imag(x)) 新グラフィックスウインドウ生成
上段

ｘのREAL成分の描画（実信号）

下段

ｘのIMAG成分の描画（参考に！）

S4)　複素振幅でQPSKを実現する場合　　

-->zero=(1/sqrt(2))*(1 + %i)*ones(1,16)

-->one=(1/sqrt(2))*(-1 + %i)*ones(1,16)

-->two=(1/sqrt(2))*(-1 - %i)*ones(1,16)

-->three=(1/sqrt(2))*(1 - %i)*ones(1,16)

-->X=[zero one two three]

-->plot2d(real(X), imag(X), style=-1)
-->square(-1,-1,1,1) 　
　

　

コンスタレーションの４点の複素数を生成。１シンボル１６点としている。

４シンボルを連結

コンスタレーションをプロット
style = 0 , -1, -2, -3, ...で点の形が変わる。

　

-->scf();

-->n=0:1/16:4-1/16;

-->x=X .* exp(%i*2*%pi*n)
-->param3d(real(x), imag(x),n) 　
　

複素振幅 X と回転を乗算する。乗算することで、変調できる。

回転を３次元で表示

--> scf();

-->subplot(3,1,1)

-->plot2d(n, real(X))

-->subplot(3,1,2)

-->plot2d(n, imag(X))

-->subplot(3,1,3)
-->plot2d(n, real(x)) 新グラフィックスウインドウ生成
１段目

複素振幅XのREAL成分

２段目

複素振幅XのIMAG成分

３段目

複素変調信号のREAL部＝実際の送信波形