上传者: 38742532
|
上传时间: 2021-09-13 14:14:44
|
文件大小: 23KB
|
文件类型: PDF
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)*xn,0,1);
end
h=h’ % h表示权函数与正交函数与w(x)