满意答案

1位老农民

2013.10.03

采纳率:59%    等级:12

已帮助:5911人

matlab 非线性的拟合有两个命令lsqcurvefit和lsqnonlin。这里用lsqcurvefit(lsqnonlin一样做),先介绍下lsqcurvefit(原理是最小二乘法)

已知数据点:xdata=(xdata1,xdata2,…,xdatan)

ydata=(ydata1,ydata2,…,ydatan)

lsqcurvefit用以求含参量x(向量)的向量值函数

F(x,xdata)=(F(x,xdata1),…,F(x,xdatan))T

中的参变量x(向量),使得1/2*∑[a+b*exp(-0.02*k*t(i))-c(i)]^2最小

1.先定义个函数fun

function y=fun(x,t)

y=(x(1)+x(2)*t.^x(3)).*(1-t).^x(4);

2.调用

clc

clear

tdata=[0.0015 0.0073 0.0125 0.0415 0.0725 0.1030 0.1323 0.1607 0.1879 0.2141 0.2392 0.2632 0.2862 0.3081 0.3291];

ydata=[-2.8988 -0.1635 0.0122 0.0174 0.0181 0.0176 0.0170 0.0164 0.0159 0.0153 0.0146 0.0140 0.0134 0.0129 0.0125];

x0=[1 1 1 1];

x=lsqcurvefit('fun',x0,tdata,ydata)

结果a =-12.4986 b =31.6477 m =0.2062

n =250.3035

*****有一个地方x0=[1 1 1 1];这边是x0迭代的初始值即a,b,m,n迭代的初始值,你可以自己修改一下(根据建的模型,找你认为近似的解或可行解带入可以提高结果精度)初始值的不同对解有点影响,但不大

11分享举报

matlab由x得到y,matlab 已知x、y值以及函数,求未知参数相关推荐

  1. matlab中sine wave出来的波形在示波器中是直线,已知仿真模型和Sine Wave模块的参数设置如下,则示波器的输出波形是...

    已知仿真模型和Sine Wave模块的参数设置如下,则示波器的输出波形是 答:A <UCP500>规定的交单期为"单据必须在提单日后15天内提交" ( ) 答:对 在环 ...

  2. (海伦公式)已知三角形三条边长,求面积

    (海伦公式)已知三角形三条边长,求面积 海伦公式:  S=(△)=√[p(p-a)(p-b)(p-c)]  其中p是三角形的周长的一半p=(a+b+c)/2. ----以下转自百度百科-------- ...

  3. 【C#】已知圆心上的三点求圆心

    已知圆心上的三点:,, ,求圆心: 把三点坐标圆的方程: 得: 把  展开整理得: 令 ; ; , , ,  求 求 把  移项整理得: 把  移项整理得: 还原代入式子: 整理结果 简化: 其中: ...

  4. 小程序转盘抽奖,已知圆心、半径、角度求圆上点的坐标

    因公司的需求,要在小程序语音房内部迭代一个真心话大冒险的功能,也等同于幸运转盘: 业务场景: 发送一个socket,拿到响应的结果,根据返回的数值开始启动幸运转盘,由于转盘上的人数不是固定的,所以要根 ...

  5. equals变量在前面和后面的区别,equals已知(存在实际值)的变量在前的好处

    代码 String a = null;System.out.println(a.equals("a")); //会报错,空指针异常System.out.println(" ...

  6. java正方形矩阵_已知一个NxN的矩阵A,求矩阵中所有边长为m的正方形的子矩阵

    已知一个NxN的矩阵A,求矩阵中所有边长为m的正方形的子矩阵内元素的中位数.(m 例子: // m = 3 n = 5 // 01234 ݂         //  0 11357 ݂         ...

  7. Python实现“已知三角形两个直角边,求斜边”

    用Python实现"已知三角形两个直角边,求斜边" 要求:用户输入两个直角边(数值为浮点类型),若非浮点类型,则提示用户,继续输入. 思路:伪代码描述下步骤 1.-input a ...

  8. python已知两条直角边求斜边,Python实现“已知三角形两个直角边,求斜边”

    用Python实现"已知三角形两个直角边,求斜边" 要求:用户输入两个直角边(数值为浮点类型),若非浮点类型,则提示用户,继续输入. 思路:伪代码描述下步骤 1.-input a ...

  9. 二叉树:已知先序和中序求后序,已知中序和后序求先序

    树的三种遍历方式的遍历顺序: 先序遍历:根.左子树.右子树(特点:第一个元素为根) 中序遍历:左子树.根.右子树(特点:根的两边分别为左子树和右子树) 后序遍历:左子树.右子树.根(特点:最后一个元素 ...

  10. java练习题 有50枚硬币,可能包括4种类型:1元,5角,1角,5分。已知总价值为20元。求各种硬币的数量。

    问题 有50枚硬币,可能包括4种类型:1元,5角,1角,5分.已知总价值为20元.求各种硬币的数量. 解决思路 对付这一类的问题,最简单的方式莫过于暴力穷举,就是用四个循环,一个一个试,结果满足要求输 ...

最新文章

  1. 人工智能是人性的罗夏测试
  2. asp备份恢复mssql数据库的代码
  3. 直播预告 | 清华大学博士生牛艺霖:常识知识获取与常识推理
  4. dm9000 driver 2
  5. 中小企业的当务之急:内部知识管理
  6. MySQL索引的索引长度问题
  7. 如何用Pygame写游戏(二十一)
  8. 理解点击屏幕的事件响应---对UIView的hitTest: withEvent: 方法的理解
  9. HashSet、LinkedHashSet和TreeSet三者区别与联系
  10. GAMIT数据的下载与准备
  11. 2016 年全国房价会呈什么趋势?
  12. linux清除回收站权限错误,在Ubuntu 14.04 中修复无法清空回收站的问题
  13. 简单易懂的现代魔法-递归
  14. Python中的六种基本数据类型
  15. MySQL基础 - 带搜索条件的查询
  16. reactos回顾总结(一)reactos源码编译
  17. 【现成工具】java获取国家法定节假日包含指定月份节假日和周末
  18. 企业所得税汇总纳税信息报告网上办理操作指南
  19. vscode 护眼主题 界面UI配色 语法配色
  20. 在Flutter中使用自定义字体,等宽字体

热门文章

  1. 2006---2009年杭电计算机历年研究生复试---笔试编程
  2. 爬虫代理及ssl验证
  3. form表单的一个页面多个上传按钮实例
  4. python下各种包的安装
  5. (转)AS3 Socket
  6. c++对象的内存布局2--进阶篇---C++ 对象的内存布局(上)
  7. Photoshop阴影与内阴影
  8. [网络安全自学篇] 十三.Wireshark抓包原理(ARP劫持、MAC泛洪)及数据流追踪和图像抓取(二)
  9. [Python图像处理] 十九.图像分割之基于K-Means聚类的区域分割
  10. OpenGL ES之GLSL实现仿抖音“灰度滤镜”和“颠倒滤镜”效果