ch2example5prg1.m

一些有用的MATLAB通信信号处理的源代码程序

% ch2example5prg1.m
dt=0.01;            % 仿真步进
T=15;               % 仿真时间长度
t=0:dt:T;           % 仿真计算时间序列
g=9.8;              % 重力加速度
L=1;                % 摆线长度
m=10;               % 摆锤质量
k=5;                % 空气阻力比例系数
theta0=3.1;         % 初始摆角设置
v0=0;               % 初始摆速设置
x0=[v0;theta0];     % 初始状态赋值
par=[g;k;m;L];      % 系统参数赋值
% 以欧拉算法计算并作图
[t_out, x_out]=eulerode('pendulumstateeq',t', x0, [], par);
plot(t_out, x_out(:,1),'-.k');hold on;
% 以ode45算法计算并作图对比
[t_out, x_out]=ode45('pendulumstateeq',t', x0, [], par);
plot(t_out, x_out(:,1),'-k');
xlabel('时间 t ');ylabel('线速度 m/s');
legend('欧拉算法','ode45算法');