作者:Z-HE
链接:https://zhuanlan.zhihu.com/p/36103034
来源:知乎
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

1) polyfit

代码例子如下,拟合一个3次曲线,并画图。

x = 0:1:9;

y = [0 2 4 6 8 20 12 14 16 18]

A=polyfit(x,y,3);

z=polyval(A,x);

plot(x,y,'r*',x,z,'b')

1) lsqcurvefit nlinfit

使用lsqcurvefit(最小二乘拟合)或nlinfit。下面的例子,是拟合一个3次曲线。这两个函数还可以拟合指数函数、三角函数等。

x=[0 1 2 3 4 5 6 7 8 9];

y=[0 1 2 3 4 5+10 6 7 8 9]; %这里,有一个噪音

w=[1 1 1 1 1 0.1 1 1 1 1]; %我们知道噪音的权重很低

a0=[1 1 1 1];

f=@(a,x)a(1)+a(2)*x+a(3)*x.*x+a(4)*x.*x.*x;

f1 = lsqcurvefit(f,a,x,y)

f2 = nlinfit(x,y,f,a0)

y1=f(f1,x)

y2=f(f2,x)

plot(x,y,'r*',x,y1,'b',x,y2,'g')

2)cftool

注意,上述拟合,是无权重拟合。w没有起作用,会发现拟合的曲线被噪音影响较大。下面将采用cftool拟合来降低噪音:

命令行中输入:cftool,回车,就打开了cftool窗口。

在窗口中,将X data选为x,将Y data选为y,此时就自动拟合为一个直线。注意,重要的地方来了:将Weights选择为w,会发现直线更加贴近正常点!

下列参考文章中有所有拟合类型的解释 matlab cftool用法及其菜单 [李园7舍_404]

3)使用fit 注意,这才是最重要的拟合!cftool能实现的,这里全能实现。

x=[0;1;2;3;4;5;6;7;8;9];

y=[0;1;2;3;4;5+10;6;7;8;9]; %这里,有一个噪音

w=[1;1;1;1;1;0.1;1;1;1;1]; %我们知道噪音的权重很低

a0=[1 1];

plot(x,y,'o')

fo = fitoptions('Method','NonlinearLeastSquares',...

'StartPoint',a0);

ft = fittype('a+x*b','options',fo);

[curve,gof] = fit(x,y,ft,'Weights',w)

hold on

plot(curve,'m')

legend('Data','curve')

hold off

参见:

Fit curve or surface to data

Fit type for curve and surface fitting

Create or modify fit options object

4) fit使用SmoothingSpline

x=[0;1;2;3;4;5;6;7;8;9];

y=[0;1;2;3;4;5+10;6;7;8;9]; %这里,有一个噪音

w=[1;1;1;1;1;0.2;1;1;1;1]; %我们知道噪音的权重很低

a0=[1 1];

plot(x,y,'o')

fo = fitoptions('method','SmoothingSpline','SmoothingParam',0.07,'Weights',w);

ft = fittype('smoothingspline');

f = fit(x,y,ft,fo)

y1= f(x)

hold on

plot(x,y1,'r')

legend('Data','f')

hold off

可惜,拟合代码不能转C++

1)写一个文件 fixSmoothingSpline.m,内容如下,是做一个加权样条拟合的,

function y1 = fixSmoothingSpline(x, y, w)%#codegen

f = fit(x,y,'SmoothingSpline','SmoothingParam',0.07,'Weights',w);

y1= f(x);

2)在命令行中输入coder,调出coder窗口,

3)可是,因为上述文件中带有fit字样,不支持转C。

Matlab各种拟合的更多相关文章

  1. MATLAB中拟合算法刚入门

    %%%1.拟合问题:(做预测,主要使用的范围是样本比较小,拟合效果会好,样本比较多,拟合的效果就不是很好) 1.应用预测的场景:已经知道10年的样本,预测第11年以内的数据 2.用拟合的到关系式:样本 ...

  2. matlab切比雪夫拟合

    matlab中没有切比雪夫拟合的现成算法,这里把我程序中的这部分抽出来,说一下. 1.首先是切比雪夫计算式 function [ res ] = ChebyShev(num,i) res=; else ...

  3. [matlab] 1.拟合

    x = [1 2 3 4 5 6 7 8 9 ]; y = [9 7 6 3 -1 2 5 7 20]; p=polyfit(x,y,3); %数字代表拟合函数的阶数 xi=0:0.01:10; yi ...

  4. matlab函数拟合

    1 函数拟合 函数拟合在工程(如采样校正)和数据分析(如隶属函数确定)中都是非常有用的工具.我这里将函数拟合分为三类:分别是多项式拟合,已知函数类型的拟合和未知函数类型的拟合.matlab中关于函数的 ...

  5. matlab的拟合函数polyfit()函数

    matlab的多项式拟合: polyfit()函数 功能:在最小二乘法意义之上,求解Y关于X的最佳的N次多项式函数. clc;clear; close all; x=[ ]; y=[2.7 7.4 2 ...

  6. Matlab 多项式拟合、稳健滤波等实用函数

    Function summary http://www.biomecardio.com/matlab/index.html clinspace Curvilinearly spaced points ...

  7. Matlab多项式拟合測试

    x=0:0.2:4; %生成等差数列 rnd=rand(1,size(x,2))*5; %生成一组随机数 y=x.*x.*x+x.*x+6+rnd; %生成y=x^3+x^2+6函数在垂直方向5个尺度 ...

  8. 【转】Matlab多项式拟合

    转:https://blog.csdn.net/hwecc/article/details/80308397 例: x = [0.33, 1.12, 1.41, 1.71, 2.19] y = [0. ...

  9. matlab多项式拟合以及指定函数拟合

    clc;clear all;close all;%% 多项式拟合指令:% X = [1 2 3 4 5 6 7 8 9 ];% Y = [9 7 6 3 -1 2 5 7 20]; % P= poly ...

随机推荐

  1. CSS新特性之动画

    动画是CSS3中具有颠覆性的特征之一,可以通过设置多个节点来精确控制一个或者一组动画,常常用来实现复杂的动画效果.相比较过度,动画可以实现更多变化,更多控制,连续自动(不需要鼠标经过和鼠标离开来控制) ...

  2. JavaScript基础2

    本节内容: 1:ECMA对象 2:string对象 3:Array对象 4:BOM对象 5:DOM对象之节点 6:DOM对象之EVENT事件 7:DOM节点的增删改查 8: 实例练习 http://w ...

  3. 201871010133-赵永军《面向对象程序设计(java)》第一周学习总结

    <面向对象程序设计(java)>第一周学习总结 正文开头: 项目 内容 这个作业属于哪个课程 https://www.cnblogs.com/nwnu-daizh/ 这个作业的要求在哪里 ...

  4. 攻防世界Web_php_unserialize

    本文借鉴 https://blog.csdn.net/qq_40884727/article/details/101162105 打开页面得到源码 <?php class Demo { priv ...

  5. Leetcode103_Binary-tree-zigzag-level-order-traversal

    原题: https://leetcode-cn.com/problems/binary-tree-zigzag-level-order-traversal/ 思路,如果无从下手,先要了解二叉树的广度优 ...

  6. Linux 学习记录一(安装、基本文件操作).

    Linux 名字的由来,是当时作者将初版的 Linux 发布在网上,供别人下载完善,而那个核心文件夹就叫 Linux,就这么叫着了.而为什么 Linux 的吉祥物是一只企鹅呢?是因为当时大家要发行稳定 ...

  7. IT兄弟连 Java语法教程 逻辑运算符

    表8中显示的布尔逻辑运算符只能操作布尔类型的操作数,所有的二元逻辑运算符都可以组合两个布尔值,得到的结果为布尔类型. 表8  布尔逻辑运算符 布尔逻辑运算符”&“.”|“以及”^“,都会布尔值 ...

  8. gitlab的安装配置与简单使用

    安装 gitlab,建议系统内存 6G ,不然会报错. 一.如何安装 gitlab 下载 gitlab 的 RPM 包 https://packages.gitlab.com/gitlab/gitla ...

  9. springcloud微服务多节点高性能、高可用、高并发部署

    1. 共有三个服务 discovery服务,domain服务,gateway服务. discovery服务是用来注册其他服务的,作为服务治理用. domain服务是主业务服务. gateway服务是所 ...

  10. Linux帮助——常用命令

    Linux帮助——常用命令 摘要:本文主要学习了Linux系统中常用的一些命令. uname命令 uname命令可以显示电脑以及操作系统的相关信息. 基本语法 uname [选项] 选项说明 -a:显 ...