フラット フェージング チャネル

変調次数を 4、変調後のシンボル レートを 45 bps、周波数分離を 200 Hz に、それぞれ設定します。

M = 4;           % Modulation order
symbolRate = 45; % Symbol rate (bps)
freqSep = 200;   % Frequency separation (Hz)

シンボルあたりのサンプル数パラメーター sampPerSym を、変調次数、周波数分離およびシンボル レートの関数として計算します。出力信号のエイリアシングを回避するには、sampPerSym と symbolRate の積が M と freqSep の積より大きくなければなりません。FSK 出力信号のサンプルレートを計算します。

sampPerSym = ceil(M*freqSep/symbolRate);
fsamp = sampPerSym*symbolRate;

FSK 変調器を作成します。

fskMod = comm.FSKModulator(M, ...
    FrequencySeparation=freqSep, ...
    SamplesPerSymbol=sampPerSym, ...
    SymbolRate=symbolRate);

3 パス フェージング チャネルのパス遅延と平均パス ゲインを設定します。

pathDelays = [0 3 10]*1e-6; % Discrete delays (s)
avgPathGains = [0 -3 -6];   % Average path gains (dB)

慣例により、最初のパスの遅延は一般にゼロに設定されます。それ以降のパスの場合、1 マイクロ秒の遅延はパス長では 300 m の違いに相当します。チャネルの平均遅延プロファイルは、パス遅延とパス ゲインの影響を受けます。

定義されたパラメーターを使用してレイリー チャネルを作成します。Visualization プロパティを設定して、インパルス応答と周波数応答を表示します。

channel = comm.RayleighChannel( ...
    SampleRate=fsamp, ...
    PathDelays=pathDelays, ...
    AveragePathGains=avgPathGains, ...
    MaximumDopplerShift=0.01, ...
    Visualization='Impulse and frequency responses', ...
    SamplesToDisplay='10%');

ランダム データ シンボルを生成して FSK 変調を適用します。

data = randi([0 3],2000,1);
modSig = fskMod(data);

FSK 変調信号とフラット フェージング チャネル フィルター処理後の受信信号をプロットするスペクトル アナライザー オブジェクトを作成します。

sa1 = spectrumAnalyzer( ...
    SampleRate=fsamp, ...
    ChannelNames=["FSK modulated signal", ...
        "Flat fading channel"]);

変調信号は、200 Hz で区切られ、それぞれが約 20 dBm のピーク強度をもつ 4 つのトーンで構成されています。

信号をレイリー フェージング チャネルを通して渡し、25 dB の S/N 比をもつ AWGN を適用します。

snrdB = 25;
raylFadedData = channel(modSig);

rxSig = awgn(raylFadedData,snrdB);

インパルス応答と周波数応答から、チャネルがフラットであるかのように動作していることが分かります。このようにフラットな応答になるのは、信号帯域幅 800 Hz がコヒーレンス帯域幅 50 kHz よりもはるかに小さいためです。

チャネル フィルター処理の前と後で、FSK 変調信号の信号スペクトルをプロットします。FSK 信号で構成される 4 つのトーンは、同一の周波数分離とピーク強度レベルを互いに相対するように維持します。絶対ピーク強度レベルは、フェージング チャネルによって低下しています。

sa1(modSig,rxSig)
release(sa1)

周波数選択性フェージング

シンボル レートを 45 kbps、周波数分離を 200 kHz に上げます。シンボルあたりのサンプル数とサンプルレート パラメーターの新しい値を計算します。FSK 変調器オブジェクトを解放し、その構成を更新します。

symbolRate = 45e3;
freqSep = 200e3;

release(fskMod)
fskMod.SymbolRate = symbolRate;
fskMod.FrequencySeparation = freqSep;

FSK 変調信号と周波数選択性フェージング チャネル フィルター処理後の受信信号をプロットするスペクトル アナライザー オブジェクトを作成します。

sampPerSym = ceil(M*freqSep/symbolRate);
sa2 = spectrumAnalyzer( ...
    SampleRate=sampPerSym*symbolRate, ...
    ChannelNames=["FSK modulated signal", ...
        "Frequecy-selective fading channel"]);

送信データに FSK 変調を適用します。

modSig = fskMod(data);

チャネル オブジェクトを解放し、チャネルのサンプル レート プロパティを更新します。レイリー フェージング チャネルを通して信号を渡し、AWGN を適用します。インパルス応答と周波数応答から、マルチパス フェージングが周波数選択的であることが分かります。

release(channel)
fsamp = sampPerSym*symbolRate;
channel.SampleRate = fsamp;

rxSig = awgn(channel(modSig),25);

チャネル フィルター処理の前と後で、FSK 変調信号の信号スペクトルをプロットします。スペクトルの形状はフラット フェージングの場合と同じですが、今度は 4 つのトーンが 200 kHz で区切られています。4 つのトーンはまだ識別可能ですが、周波数選択性フェージングによって、トーンの相対的なピーク強度レベルが異なっています。800 kHz の信号帯域幅が、50 kHz のコヒーレンス帯域幅よりも大きくなっています。

sa2(modSig,rxSig)
release(sa2)