电力系统仿真课程论文

subplot(2,1,2); stem(n,imag(x)); xlabel('时间(k)'); ylabel('幅值f(k)'); title('虚部');

运行结果为：

3.2用MATLAB实现信号的时域运算

信号的时域运算包括信号的相加、相乘、数乘、微分、积分。

（1）信号的相加和相乘：已知信号f1(t)和f2(t)，信号相加和相乘记为f(t)?f(t)?f2(t)；1f(t)?f(t)?f2(t)。 1（2）信号的微分和积分：对于连续时间信号，其微分运算是用diff函数来完成的，其语句格式为：diff(function,’variable’,n)，其中function表示需要进行求导运算的信号，或者被赋值的符号表达式；variable为求导运算的独立变量；n为求导的阶数，默认值为求一阶导数。连续信号的积分运算用int函数来完成，语句格式为：diff(function,’variable’,a,b)，其中function表示需要进行被积信号，或者被赋值的符号表达式；variable为求导运算的独立变量；a,b为积分上、下限，a和b省略时为求不定积分。

3.2相加

例1：t=0:0.001:6; y1=2*rectpuls(t-3,2); y2=2*rectpuls(t-3,6);

f1=y1+y2; plot(t,f1); grid on;

axis([0 7 0 5]);

函数f1(t)由两个矩形脉冲信号相加。运行程序，得如下图波形。

3.3用MATLAB实现信号的时域变换（参数变化，分析波形变化） 信号的时域变换包括信号的反转、使移、展缩、倒相

信号的平移包含信号的左移与右移，信号的翻转包含信号的倒相与折叠，平移和翻转信号不会改变信号f(t)的面积和能量。信号的尺度变换是对信号f(t)在时间轴上的变化，可使信号压缩或扩展。f(at)将原波形压缩a倍，f(t/a)将原波形扩大a倍。 设f(t)=(1+t/2)*[&(t+2)-&(t-2)],用MATLAB来求f(t+2),f(t-2),f(-t),f(2t),并绘出其时域波形。 程序如下： syms t

f=sym('(t/2+1)*(heaviside(t+2)-heaviside(t-2))') subplot(2,3,1),ezplot(f,[-3,3]) title('f(t)') y1=subs(f,t,t+2)

subplot(2,3,2),ezplot(y1,[-5,1]) title('f(t+2)') y2=subs(f,t,t-2)

subplot(2,3,3),ezplot(y2,[-1,5]) title('f(t-2)') y3=subs(f,t,-t)

subplot(2,3,4),ezplot(y3,[-3,3]) title('f(-t)') y4=subs(f,t,2*t)

subplot(2,3,5),ezplot(y4,[-2,2]) title('f(2t)') 运行结果：

3.4用MATLAB实现信号简单的时域分解

信号的交直流分解即将信号分解成直流分量和交流分量两部分之和，其中直流分量定义为 ：

程序如下：

t=-1:0.001:1; f=sin(2*pi*t)+2; g=mean(f); h=f-g;

subplot(311); plot(t,f);

axis([-1 1 0.5 3.5]); subplot(312); plot(t,g);

axis([-1 1 1.5 2.5]); subplot(313); plot(t,h);

axis([-1 1 -1.5 1.5]); 运行结果为：

3.5用MATLAB实现连续时间系统的卷积积分的仿真波形