1条回答 默认 最新
- 关注
码龄 粉丝数 原力等级 --
- 被采纳
- 被点赞
- 采纳率
CSDN专家-黄老师 2021-11-20 09:29最佳回答 专家已采纳插值基函数
n 次代数插值问题的解可表示为pn(x)=sumlimitsnk=0lk(x)yk,
MATLAB 程序
function yi = LagInterp(x,y,xi) % Lagrange 插值多项式,调用格式为 % yi = LagInterp(x,y,xi) % 其中 % x 为插值节点,y为节点处函数值, % xi 为为估计函数自变量,yi 为xi处函数估计值 % n = length(x); m = length(xi); p = zeros(n,m); for k = 1:n t = ones(n,m); for j = 1:n if j~=k if abs(x(k) - x(j))<eps error('% 输入的插值节点必须互异!'); end t(j,:) = (xi - x(j))/(x(k) - x(j)); end p(k,:) = prod(t); end yi = y*p; end x = [1.1275, 1.1503, 1.1735, 1.1972]; y = [0.11910, 0.13954, 0.15932, 0.17903]; xi = [1.1300 1.1600]; yi = LagInterp(x,y,xi)
采纳该答案 已采纳该答案 专家已采纳评论解决 1 无用打赏举报微信扫一扫
分享评论登录 后可回复...
报告相同问题?
提交
相关推荐 更多相似问题
- 2021-11-20 00:08回答 1 已采纳 插值基函数n 次代数插值问题的解可表示为 pn(x)=sumlimitsnk=0lk(x)yk, MATLAB 程序 function yi = LagInterp(x,y,xi) % La
- 2018-03-28 08:18回答 1 已采纳 第一步 将function yy=Lagrange(x,y,xi)m=length(x);n=length(y);if m ~= n , error('向量x与y的长度必须一致'); end; s=
- 2021-10-09 14:09回答 1 已采纳 你好相关问题请参考我之前的回答: 利用Matlab解决以下实验问题-编程语言-CSDN问答 CSDN问答为您找到利用Matlab解决以下实验
- 2021-05-29 21:27拉格朗日插值多项式拟合一组点 LAGRANG(X,Y,N,XX) 其中 X 和 Y 是定义一组 N 个点的行向量,使用拉格朗日方法在 X 中找到通过这些点的 N 阶多项式。
- 2021-06-01 00:22拉格朗日插值多项式
- 2021-04-01 10:43@二十三的博客 function yy = lag(x,y,xx) m=length(x); n=length(y); if m~=n, error(’’) end s=0; for i=1:n t=ones(1,length(xx)) ; for j=1:n if j~=i, t=t.(xx-x(j))/(x(i)-x(j)); end end s=s+ty(i); ...
- 2021-10-06 23:54回答 1 已采纳 我们观察aa,实际描述意思就是[A, G'; G, 0]C=[b;1]求C这个向量吧,光看最后一行的话,就是说GC=1举个例子,如果我们要拟合的为三次多项式,多项式表达式为c0+c1x+c2x^2+c
- 2021-10-26 02:07回答 2 已采纳 你好,xy构成的多项式也可以写,但是比较繁琐,而且求解很难满足最后一个条件,全局大于0,下面是我的一维拓展方案:主函数: %% 为了得到分布,假设了一个协方差矩阵 mu=[20,180];%数学期望
- 2022-03-22 10:26回答 1 已采纳 这个你先把自由度、自由度一阶导分别用theta、omega表示,然后用subs替换,我举个很简单的例子吧,比方说单摆的自由摆动方程的推导 代码实现如下 很容易理解吧,有帮助的话,两个回答都采纳了吧
- 2020-12-02 08:06实验四用MATLAB实现拉格朗日插值、分段线性插值.pptx
- 2021-05-30 02:52使用拉格朗日多项式插值法近似点定义函数
- 2011-04-11 22:25求已知数据点的拉格朗日插值多项式 matlab源代码
- 2021-11-16 16:55回答 1 已采纳 你好,你给出的代码是拉格朗日(Lagrange)插值的代码,不是拟合的代码,拟合建议用polyfit函数,因为你如果用拉格朗日多项式插值的话,虽然可以保证多项式经过每一个点,但是当插值点个数过多的话,
- 2021-10-01 23:27回答 3 已采纳 你好,之前也提醒过题主用fmincon函数,现在把程序发一下,我是按照fmincon求解的(额外加了非线性约束保证函数值全体大于0),但是初值选取使用的是拉格朗日乘子法加最小二乘法【这样计算比较稳定】
- 2021-07-02 15:52回答 5 已采纳 错误原因:输出的数据太多了。或者参数不对。
- 2021-04-20 02:00weixin_39682477的博客 %拉格朗日插值多项式 利用矩阵求解x=1:0.2:3;%已知数据点x坐标向量:xy=sin(x);%已知数据点x坐标向量:yx1=1.1:0.2:3.1;%插值点的x坐标:x1L=zeros(11,11);%另L矩阵为0for i=1:11A=ones(10,1);%另A矩阵为10行1列的...
- 2021-06-01 21:33第一个脚本 barylag.m 对一组给定数据执行重心拉格朗日插值。 这种方法遵循 LN Trefethen 的一篇论文(脚本注释中的引用),并且比之前的脚本快得多。 此外,它已被矢量化以进一步缩短计算时间。 第二个脚本 ...
- 2020-05-15 21:02拉格朗日插值及其画图 1
- 2021-10-30 16:46回答 1 已采纳 你好同学,建议从积分得到的非线性方程考虑,或者先由最小二乘得到正数表达式,然后用积分考虑改变前面的系数除以一个比例就行(实质是改变a0)。方法1 非线性方程组方法: x = [0:1:40]; % x
- 2020-10-11 16:04寻欢公子的博客 一 ![在这里插入图片描述](https://img-blog.csdnimg.cn/20201011160101791.JPG?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80...color_FFFF...
- 没有解决我的问题, 去提问