![](https://img-blog.csdnimg.cn/20201014180756930.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数值分析
暴杀汉堡
赢麻了
展开
-
matlab上的简单拉格朗日插值代码
备注: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-1 a=1;原创 2021-06-03 09:05:35 · 639 阅读 · 1 评论 -
牛顿插值的简单matlab代码
备注: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); % 差分表原创 2021-06-03 08:54:27 · 1485 阅读 · 0 评论 -
matlab上简单牛顿-柯特斯公式实现
备注: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]/2 Cotescoeff(2,:)=[1,4,1,0,0,0,0,0]/6 Cot原创 2021-06-03 08:49:25 · 4642 阅读 · 2 评论