a783208746的博客

备注:1.由于是敲实验报告自用，所以不能保证代码的通用性，具体问题还要具体分析2.这里只有函数文件，脚本文件还请读者自行建立function s=l2(m,x,y,n)#x为自变量矩阵，y为因变量矩阵，如x=[0.4,0.5,0.6,0.7,0.8];#y=ln(x),y=[-0.916291,-0.693147,-0.510826,-0.356675,-0.223144]#n为x矩阵的长度(这里是5),m为预测值的自变量s=0;w=zeros(1,n);for k=2:n-1a=1;

备注:1.由于是敲实验报告自用，所以不能保证代码的通用性，具体问题还要具体分析2.这里只有函数文件，脚本文件还请读者自行建立function y=n1(x,y,n,x1)#参数说明，x为自变量矩阵，y为因变量矩阵，如x=[0.4,0.5,0.6,0.7,0.8];y=ln(x),y=[-0.916291,-0.693147,-0.510826,-0.356675,-0.223144];n为x矩阵的长度(这里是5),x1为预测值的自变量Y= zeros(n);b=zeros(n,1);% 差分表

备注:1.由于是敲实验报告自用，所以不能保证代码的通用性，具体问题还要具体分析2.牛顿-柯特斯公式最多有7阶并且系数是确定了，这里是直接将柯特斯系数存入了矩阵3.这里只有函数文件，脚本文件还请读者自行建立function y=NewtonCotes(fun,a,b,n)sum=0;% 生成柯特斯系数矩阵Cotescoeff=zeros(7,8);Cotescoeff(1,:)=[1,1,0,0,0,0,0,0]/2Cotescoeff(2,:)=[1,4,1,0,0,0,0,0]/6Cot