爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 12462|回复: 18

[秀图] 用matlab画多项式拟合曲线

[复制链接]

新浪微博达人勋

发表于 2014-6-12 15:05:33 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
之前要画多项式拟合曲线,就是当做折线图的趋势线画的。试了grads和Excel最后还是用MATLAB实现了。
分享下,希望有所帮助吧。其实不难的,但还是分享下
>> load('E:\dy500hgt\data\1\txt\wz_yr.txt')
>> t=[1:35]
t =
  Columns 1 through 15
     1     2     3     4     5     6     7     8     9    10    11    12    13    14    15
  Columns 16 through 30
    16    17    18    19    20    21    22    23    24    25    26    27    28    29    30
  Columns 31 through 35
    31    32    33    34    35
>>plot(t,wz_yr,’-r') (表示用红色实线画图)
>> hold on (试原图保留,不至于被后面的覆盖)


>> b=polyfit(t,wz_yr',2) (其中的wz_yr'上的'是似的那个矩阵转置下,做成和t一样的矩阵类型)
b =
   -0.0101    0.3820  129.1424  (二次多项式拟合的系数)
>> y=polyval(b,t,1) (计算对应的拟合值)
y =
  Columns 1 through 5
  129.5144  129.8663  130.1980  130.5096  130.8011
  Columns 6 through 10
  131.0724  131.3237  131.5548  131.7658  131.9567
  Columns 11 through 15
  132.1275  132.2782  132.4087  132.5191  132.6094
  Columns 16 through 20
  132.6796  132.7296  132.7596  132.7694  132.7591
  Columns 21 through 25
  132.7287  132.6781  132.6075  132.5167  132.4058
  Columns 26 through 30
  132.2748  132.1237  131.9524  131.7610  131.5496
  Columns 31 through 35
  131.3179  131.0662  130.7944  130.5024  130.1903
>> plot(t,y,'b')

其中>>后面的才是在matlab命令框中需要输入的。最后放个自己做好的图。 9[TZVMGLEK%KR080T~BIH.jpg

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2014-6-12 20:27:05 | 显示全部楼层
可以写成M文件的,之前做过,有的还行,有的拟合效果很一般,和预期的很不一致……
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2014-6-12 20:27:15 | 显示全部楼层
可以写成M文件的,之前做过,有的还行,有的拟合效果很一般,和预期的很不一致……
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-13 22:07:14 | 显示全部楼层
真心不错,感谢楼主分享,正在努力学习中
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-14 09:37:58 | 显示全部楼层
二爷名声在外 发表于 2014-6-12 20:27
可以写成M文件的,之前做过,有的还行,有的拟合效果很一般,和预期的很不一致……

其实我不是很会matlab的,所以不会写m文件
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2014-6-14 09:50:40 | 显示全部楼层
shrily521 发表于 2014-6-14 09:37
其实我不是很会matlab的,所以不会写m文件

好吧…不过没关系,很好学的!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-14 10:19:22 | 显示全部楼层
拟合优度是?R2是多少呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-16 10:18:49 | 显示全部楼层
kongfeng0824 发表于 2014-6-14 10:19
拟合优度是?R2是多少呢?

这是二次多项式拟合。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-17 12:53:11 | 显示全部楼层
shrily521 发表于 2014-6-16 10:18
这是二次多项式拟合。

应该有个拟合优度才对。cftool。这样对拟合的程度才有判断
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-17 15:48:16 | 显示全部楼层
kongfeng0824 发表于 2014-6-17 12:53
应该有个拟合优度才对。cftool。这样对拟合的程度才有判断

这个我也不清楚,不好意思,有人知道的话麻烦回答下吧
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表