p1=19，p2=24。用matlab编写。第二题中求出后用h分之一来做，即保证分母阶数比分子高

中午十二点前交，急

B计算以下系统的总传递函数H(s)，分别求出冲激与阶跃响应，以及激励为f(t)=8(sint+cost)*u(t)的系统响应的波形图(0≤t≤p2s,步长0.01s)。画出零极点图并由其判断该系统是稳定的吗

function [m]=intgauss(n)    %  n代表所求节点的个数
syms x
for i=1:n
    for j=1:n
        y(i,j)=int(log(x)*x^(i-1)*x^(n-j),0,1);   % 积分的权函数不同则log(x)会改变
    end
end 
y                                 % y表示权函数与正交函数与w(x)的积分值所构成的方阵
for i=1:n
   h(i)=-int(log(x)*x^(i-1)*x^n,0,1);
end
h=h'                              % h表示权函数与正交函数与w(x)的积分值所构成列向量
m=y^(-1)*h;                       % m是系数向量
w=x^n;
for i=1:n
    w=w+m(i)*x^(n-i);
end
x=solve(w==0)                     % x表示的是节点位置
for i=1:n
    for j=1:n
        l(i,j)=x(j)^(i-1);
    end
   z(i)= y(i,n);
end
z=z';
A=l^(-1)
A=l^(-1)*z                       % A代表系数
for i=1:n
   X = sprintf('%f*f(%f)',A(i),x(i));
   disp(X)
end