MATLAB 插值与拟合

MATLAB 插值与拟合一.1维插值函数1.通用接口:进行1维插值:vq=interp1([x,]v,xq[,method,extrapolation)#参数说明: x,v:分别指定样本点的x,y值;x默认为[1,2…length(v)](v为向量)或[1,2…size(v,1)](v为高维数组) #注意:要求x是单调的 xq:指定查询点(即求该点的插值) method:指定插值方法;可为’linear'(默认值;线性插值)/’nearest'(最近邻点插值)/’next'(下1个邻点

大家好,又见面了,我是你们的朋友全栈君。

一.插值
1.1维插值函数
(1)通用接口:

进行1维插值:vq=interp1([x,]v,xq[,method,extrapolation)
  #参数说明:
  	x,v:分别指定样本点的x,y值;x默认为[1,2...length(v)](v为向量)或[1,2...size(v,1)](v为高维数组)
  	  #注意:要求x是单调的
  	xq:指定查询点(即求该点的插值)
  	method:指定插值方法;可为'linear'(默认值;分段线性插值)/'nearest'(最近邻点插值)/'next'(下1个邻点插值)/'previous'(上1个邻点插值)/'pchip'(保形分段3次插值)/'cubic'(保形分段3次插值)/'v5cubic'(用于MATLAB5的3次卷积)/'makima'(修正了Akima公式的3次埃尔米特插值)/'spline'(3次样条插值)
  	  #注意:x等距时可使用快速插值,此时指定method为'*linear'等
  	extrapolation:指定外插策略;可为'extrap'(使用method计算)/标量值(返回该值)/NaN(默认值;返回NaN)
  	  #即当xq落在x外,如何处理
  	vq:返回插值

#实例:
>> x=[1,2,3,4];
>> v=[1,-1,3,9];
>> interp1(x,v,3.5)

ans =

     6

>> v=[1,1,1;-1,-1,-1;3,3,3;9,9,9];
>> interp1(v,3.5)

ans =

     6     6     6

(2)3次样条插值:

进行3次样条插值:s=spline(x,y[,xq])
  #参数说明:
  	x,y:分别指定样本点的x,y值

######################################################################################################################

csape()

(3)分段3次埃尔米特插值(保型分段3次插值):

pchip()

(4)分段多项式插值:

mkpp()

2.2维插值
(1)通用接口:

interp2()

二.拟合
1.1元多项式曲线拟合(Polynomial Curve Fitting):

进行1元多项式曲线拟合:[p,S,mu] = polyfit(x,y,n)
  #参数说明:
  	x,y:分别指定数据的x,y值
  	n:指定多项式的次数
  	p:返回按降幂排列的多项式系数
  	S:返回1个结构体,可用作polyval()的输入来获取误差估计值
  	mu:返回1个二元素向量,mu(1)为mean(x),mu(2)为std(x)
  	  #要求返回mu时,polyfit()会对x进行Z-Score Standarlization,这会改善多项式和拟合算法的数值属性

#实例:
>> polyfit([1,2,3],[1,4,9],1)

ans =

    4.0000   -3.3333

>> [p,S,mu]=polyfit([1,2,3],[1,4,9],1)

p =

    4.0000    4.6667


S = 

  包含以下字段的 struct:

        R: [2×2 double]
       df: 1
    normr: 0.8165


mu =

     2
     1

2.多元线性回归:

进行多元线性回归:[b,bint,r,rint,stats] = regress(y,X[,alpha])
  #设多元线性方程为y=b1*x1+b2*x2+...+bn*xn
  #参数说明:
  	y:指定数据点的y坐标;为column array
  	X:指定数据点的x坐标;为matrix(每列对应1个变量)
  	  #y,X的行数应相等;X(:,1)对应x1,X(:,2)对应x2...;若方程有常数项,应有某列全为1
  	alpha:指定置信水平;默认为0.05
  	  #将使用100*(1-alpha)%置信水平
  	b:返回超平面的各个系数构成的向量
  	  #b(1)为b1,b(2)为b2...
  	bint:返回1个矩阵,其中包含系数估计值的95%置信区间
  	r:返回由残差组成的向量
  	rint:返回1个矩阵,其中包含可用于诊断离群值的区间
  	stats:返回1个向量,其中包含R2统计量/F统计量及其p值,以及误差方差的估计值
  	  #X应包含1个全为1组成的列,以便正确计算模型统计量

3.曲线拟合工具箱(Curve Fitting Toolbox):

>> cftool%见下图

在这里插入图片描述

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/136219.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • django restful API 代码自动生成_restful接口设计规范

    django restful API 代码自动生成_restful接口设计规范restful接口规范什么是接口规范?接口规范就是为了采用不同的后台语言,也能使用同样的接口获取到同样的数据。如何写接口:接口规范是规范化书写接口的,写接口要写url、响应数据​注:如果将请求参

  • onedrive免费扩容25t_onedrive怎么免费扩容1T

    onedrive免费扩容25t_onedrive怎么免费扩容1TOneDrive存储我们都知道没有开office365,自己onedrive的储存空间只有5GB,onenote做笔记以及用onedrive同步文档空间不够,但是又不想开office365;所以在网上看到别人说onedrive可以推荐别人注册,可以扩容10GB;加起来一共15GB,用来做笔记完全够用;或许有大佬会说可以弄到Office教育版的微软账号,有5T或1T的存储空间,但是这个会涉及到账号里面文档的安全性,这种账号是属于教育机构的,全局管理员可以有权查看里面储存的文件并且有权删去账号,这样的

  • 常用jar包下载地址

    常用jar包下载地址http://www.java2s.com/Code/Jar/CatalogJar.htm

  • matlab 定位_matlab读取mp3音频文件

    matlab 定位_matlab读取mp3音频文件PAGEPAGE1毕业设计(论文)基于MATLAB的声源定位系统专业年级07级应用物理学学号姓名指导教师评阅人二○一一年六月中国南京河海大学本科毕业设计(论文)任务书(理工科类)Ⅰ、毕业设计(论文)题目:基于matlab系统声源识别与实时定位初步研究Ⅱ、毕业设计(论文)工作内容(从综合运用知识、研究方案的设计、研究方法和手段的运用、应用文献资料…

  • pycharm突然很卡_pycharm启动慢

    pycharm突然很卡_pycharm启动慢使用pycharm时常出现 theIDEisrunninglowonmemory的问题,表示pycharm这款IDE使用内存不足,需要在系统内存充足的情况下扩充IDEmemory。首先,打开File—>Appearance—>WindowsOptions—>选中showmemmoryindicator—>OK右下角会出现 …

  • windows环境下tomcat部署war包[通俗易懂]

    windows环境下tomcat部署war包[通俗易懂]1.从Eclipse中将项目打包成war包,Runas>>MavenInstall,然后将war包放在tomcat\webapps下,如图所示2.然后运行tomcat\bin目录下的startup.bat(运行后悔在webapps下生成相应名称的目录)3.运行成功,使用Postman测试成功…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号