デジタル信号シミュレーション実験-1つの離散時間信号とシステム時間領域分析の実験

1.実験の目的

(1)離散時間信号の生成と基本操作に
精通している(2)離散時間システムの時間領域特性に精通している
(3)畳み込み法を使用してシステムの時間領域特性を観察および分析する

2.実験原理

(1)典型的な離散時間信号
ここに写真の説明を挿入
(2)シーケンスの基本操作
ここに写真の説明を挿入
ここに写真の説明を挿入
(3)線形畳み込み
ここに写真の説明を挿入
(4)私たちが主に研究する線形時間不変離散時間システムは次の形式を使用します
ここに写真の説明を挿入

3.実験内容

(1)Matlabを使用して、典型的な離散時間信号を生成し、それらのグラフを描画します。
ここに写真の説明を挿入
(2)Matlabを適用して、2つの有限長シーケンスの線形畳み込みを計算します。
(3)差分方程式で表される因果線形時間不変離散時間システムは、
ここに写真の説明を挿入
(4)入力信号が
ここに写真の説明を挿入
(5)ECG(パルス、EEG)信号の負荷である場合、その時間領域波形を描画します。

4.実験レポートの要件

(1)実験報告書には、実験の目的と実験原理の要点を簡潔に記載してください。
(2)実験中に記録された各信号の時間領域波形を実験レポートに添付し、結果のグラフを分析し、各信号のパラメーター変更がその時間領域特性に与える影響を説明します。
(3)実験の主な結論を要約します。

Matlabプログラム1_a:

Problem1_a
clear
N=input('Type in the length of sequence=');%% 输入一个N
n=-(N-1):1:N-1;
x1=[zeros(1,N-1),1,zeros(1,N-1)];%%zeros()零矩阵
stem(n,x1);%%绘制火柴梗,产生离散信号
xlabel('Time index n');
ylabel('Amplitude');
title('unit sample sequence LEI');

操作結果1_a:

ここに写真の説明を挿入

Matlabプログラム1_b:

Problem1_b
clear
N=input('Type in the length of sequence=');%% 输入一个N
n=-(N-1):1:N-1;
x1=[ones(1,N-1),0,ones(1,N-1)];%%ones()1矩阵
stem(n,x1);%%绘制火柴梗,产生离散信号
xlabel('Time index n');
ylabel('Amplitude');
title('unit step sequence LEI');

操作結果1_b:

ここに写真の説明を挿入

Matlabプログラム1_c:

Problem1_c
clear
N=input('Type in the length of sequence=');%% 输入一个N
n=0:1:N-1;
x1=sin(pi/6*n);%%sin(pi/6)
stem(n,x1);%%绘制火柴梗,产生离散信号
xlabel('Time index n');
ylabel('Amplitude');
title('sinusoidal sequence LEI');

操作結果1_c:

ここに写真の説明を挿入

Matlabプログラム2:

Problem2
x=[0 1 2 3 4 5];%%任意有限序列
y=[5 4 3 2 1 0];
z=conv(x,y)%%计算线性卷积
stem(y)

運用結果2:

ここに写真の説明を挿入

Matlabプログラム3:

Problem3
N=41;
a=[0.9,-0.45,0.35,0.002];
b=[1,0.71,-0.46,-0.62];
x1=[1 zeros(1,N-1)];%%ones()1矩阵
x2=ones(1,N);%%ones()1矩阵
k=0:1:N-1;
h=filter(a,b,x1);%%实现差分方程的仿真
y=filter(a,b,x2);
subplot(2,1,1);
stem(k,h,'.');%%绘制火柴梗,产生离散信号
xlabel('n');
ylabel('unit sample sequence');
title('Made by LEI');
subplot(2,1,2);
stem(k,y,'.');%%绘制火柴梗,产生离散信号
xlabel('n');ylabel('unit step sequence');grid on;
Problem4
n=40;
k=0:1:n-1;
num=[0.9 -0.45 0.35 0.002];
den=[1 0.71 -0.46 -0.62];
y1=impz(num,den,n);%%系统的冲激响应
x=[ones(1,n)];
y2=filter(num,den,x);%%实现差分方程的仿真
figure(1)
subplot(211)
stem(k,y1);%%绘制火柴梗,产生离散信号
xlabel('Time index n');ylabel('Amplitude');
title('unit sample response LEI');
subplot(212)
stem(k,y2);%%绘制火柴梗,产生离散信号
xlabel('Time index n');ylabel('Amplitude');
title('unit step response LEI');
x1=[1 2 0 -0.5];
y=conv(y1,x1)   %%计算卷积
figure(2)
stem(y)%%绘制火柴梗,产生离散信号
xlabel('Time index n');ylabel('y[n]');
title('time domain waveform LEI');

運用結果3:

ここに写真の説明を挿入
ここに写真の説明を挿入
能力に限りがあり、説明がわかりません。何か問題が発生した場合は、メッセージやプライベートメッセージを残すことができます。プログラムファイルはパッケージ化され、後で誰もが学習して使用できるようにアップロードされます。

この記事は皆様のお役に立てれば幸いです。もちろん、上記に問題がある場合は訂正してください。

決定の高さを共有し、ギャップを開くことを学ぶ

おすすめ

転載: blog.csdn.net/qq_42078934/article/details/109251683