matlab 多项式拟合程序,科学网—数值分析--多项式拟合程序matlab程序 - 殷春武的博文...

%%%程序编写者  西北工业大学自动化学院    Email: yincwxa2013@mail.nwpu.edu.cn

%%  All rights reserved

clear

clc

x=input('输入节点坐标x=')

y=input('输入节点坐标函数值y=')

w=input('输入权重向量w=')

n=input('输入拟合最高次数n=')

for i=1:n+1

for j=1:n+1

p(i,j)=(w.*(x.^(i-1)))*(x.^(j-1))';

end

end

for i=1:n+1

q(i)=(w.*y)*(x.^(i-1))';

end

C=inv(p)*q'

f=poly2sym(C)

%%%

输入节点坐标x=[1 1 2 3 4 ]

x =

1     1     2     3     4

输入节点坐标函数值y=[2 3 5 6 9]

y =

2     3     5     6     9

输入权重向量w=2

w =

2

输入拟合最高次数n=2

n =

2

C =

0.9231

1.5641

0.1026

f =

32477881928153/35184372088832*x^2+1719748956265/1099511627776*x+7217307095145/70368744177664

>>

转载本文请联系原作者获取授权,同时请注明本文来自殷春武科学网博客。

链接地址:http://blog.sciencenet.cn/blog-748007-751138.html

上一篇:数值分析--雅克比迭代法、高斯赛德尔迭代法、超松弛迭代法matlab

下一篇:数值分析--非线性方程组不动点迭代法matlab程序

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
拟合功能说明: 对于形如y=a1*x^n+a2*x^(n-1)+...+an*x+a(n+1)的函数,若a1,a2,...,a(n+1)未知,则可以由n+1对不同的(x,y)值求出这些a的值。更进一步的说,给你一条X-Y轴平面里的曲线,或是,已知一系列的(x,y)值,要你用一个近似的一元函数拟合它,该怎么办?这个时候你不知道该一元函数的最高次数是多少,即不知道未知数系数a1,a2,...,a(n+1)的个数。在这里显然要先假设一个n的值,再来算出这些系数值:有n+1个未知数系数要求,就得有n+1对(x,y)的值,这样可以精确算出唯一的解。但是,如果给你的(x,y)值对数多于n+1对该怎么办(因为你碰到的情况几乎都是这样的)?你算出来的近似函数得兼顾所有的(x,y)取样点,使得误差最小。你不能只取所有(x,y)取样点中的n+1对计算,而扔掉了多出来的那些取样点。该程序的拟合功能就是完成这个任务的。它会计算所有的取样点,最后得出一个误差最小的函数。 执行过程如下: 先读取右边两个编辑框内的(x,y)取样点,再弹出一个对话框让你假设一个一元函数,即最高次数n,然后程序会开始计算。注意,若n设的太大,计算过程可能会很久,会造成程序假死现象。得出的结果会直接以函数的形式显示在下面的长条形编辑框内。你可以通过点击“←画图”和“画图→”来进行比较拟合的近似程度。注意,你设置的n值不是越大,就拟合得越好,也不是越小,拟合得越好,原因你细想一下就会明白。程序所用的计算原理是最小二乘法,所用的计算方法是矩阵运算。你可以百度相关知识。 该程序的另外一个“辨识”功能,用到的也是最小二乘法,传递函数等相关知识是大学工科学生会学到的。 其他不懂的问题请联系作者:hastings1986@163.com。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值