利用面积法求高阶微分方程系数

clear;clc;
syms s;dt=0.01;xo=1;
fprintf('------------第一部分-------给定参数生成数据----------------\n');
fprintf(' 假设输入 Xo=1\n');
syms s;dt=0.01;xo=1;
fprintf('------------第一部分-------给定参数生成数据----------------\n');
fprintf(' 假设输入 Xo=1\n');
n=input('1. 请给定系统的阶数: n=');
T=zeros(1,n);
fprintf('2. 根据所定阶次，请依次给出时间常数的值：\n');
B=1;
for m=1:n
fprintf(' 2.%d 请给定时间常数的值：T(%d)',m,m);
T(m)=input('=');
B=(T(m)*s+1)*B;
end
L=xo/(s*B);y=ilaplace(L);t=0:500;y=double(subs(y,t));
ti=0:dt:500;yi=zeros(size(ti));yi=interp1(t,y,ti);
fprintf('3. 采 样 其 中 几 个 点 的 数 据 ：
\n------yi(0)=%.4f\n------yi(100)=%.4f\n------yi(500)=%.4f\n'...
,yi(1),yi(10000),yi(50000));
%-------------至此，已生成一组给定阶数的响应数据
fprintf('------------------------- 第 二 部 分 ------------------------\n
n=input('1. 请给定系统的阶数: n=');
T=zeros(1,n);
fprintf('2. 根据所定阶次，请依次给出时间常数的值：\n');
B=1;
for m=1:n
fprintf(' 2.%d 请给定时间常数的值：T(%d)',m,m);