连续系统的微分方程matlab实例.rar_matlab路面_时域路面_路面_路面时域_路面模型

%% 时域分析 %% %%%%% 数字信号处理与matlab 第七章7.3 124页 对于连续系统，输入信号f(t),冲激响应h(t), %%%%% 零状态响应为y(t),则有y(t)=f(t)*h(t),其中h(t)包含了连续系统的固有特性，与系统 %%%%% 的输入无关。 % 连续系统的微分方程：2y''(t)+y'(t)+8y(t)=f(t)，利用matlab绘制该系统的冲激 % 响应和阶跃响应的波形。 b=[1]; a=[2 1 8]; subplot(1,2,1); impulse(b,a);%绘出由向量a和b表示的连续系统在指定时间范围内的冲激响应h(t)的时域 %波形图，并求出指定时间范围内冲激响应的数值解。 subplot(1,2,2); step(b,a)%绘出由向量a和b表示的连续系统的阶跃响应g(t)在指定时间范围内的波形图 %并求出其数值解。 % %% y''(t)+5y'(t)+6y(t)=3f''(t)+2f(t) c=[1 5 6]; d=[3 2]; impulse(d,c); %% 利用MATLAB求LTI【线性时不变系统（Linear time invariant）】连续系统的响应 %****线性时不变系统就是系统的参数不随时间的变化而变化 % y''(t)+2y'(t)+y(t)=f'(t)+2f(t),若要求当输入信号为f(t)=exp(-2*t)*e(t) %时该系统的零状态响应y(t),命令如下： a=[1 2 1]; b=[1 2]; p=0.5;%定义取样时间间隔 t=0:p:5;%定义时间范围 x=exp(-2*t);%定义输入信号 lsim(b,a,x,t);%对系统输出信号进行仿真 % a和b是描述系统的微分方程系数，表示该系统的两个行向量，x和t则是表示输入信号 % 的行向量，其中，t为输入信号时间范围的向量，x则是输入信号在向量t定义的时间点上 % 的取样值。 %%%%% y=lsim(b,a,x,t)与impulse和step不同的是，该格式并不绘出系统的零状态响应的曲线 %%%%% 而是求出与向量t定义的时间范围相一致的系统零状态响应的数值解。 %% 频域分析 %%%%% 数字信号处理与matlab 第八章 周期信号频域分析及MATLAB实现 137页 %% 利用FFT实现周期信号的频谱分析；161页 %% 第九章 连续时间信号的频域分析及MATLAB实现 168 页 % example 设f(t)=（1/2）*exp(-2*t)*e(t),试画f(t)及其幅频图 % 这里采用符号变量定义符sym,并用符号函数作图函数ezplot()绘图。实现该过程的命令如下： syms t v w x;%生成符号表达式 x=1/2*exp(-2*t)*sym('Heaviside(t)'); F=fourier(x); subplot(211); ezplot(x); subplot(212); ezplot(abs(F));