爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 9922|回复: 8

[作图] ncl拟合二次曲线问题

[复制链接]
发表于 2017-4-9 16:25:32 | 显示全部楼层 |阅读模式

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

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

x
我想用txt文本数据做散点图,并拟合绘制成一条二次曲线,为啥会出现这么多线??
代码如下:
  fn = "11.txt"
  ncols = numAsciiCol(fn)
  nrows = numAsciiRow(fn)
  data1 = asciiread(fn ,(/nrows,ncols/),"float")
; print(data1)

  x=data1(:,0)
  ;print(x)
  ts=data1(:,1)
  ;ts=runave(ts,20,0)
  rc = lspoly(x,ts,1,3)
  print(rc)
  data=new((/2,dimsizes(ts)/),"float")
  data(0,:)=ts
  data(1,:)=rc(0)+rc(1)*x+rc(2)*x^2
  ;print(data)

  wks   = gsn_open_wks ("png","scatter2")          ; send graphics to PNG file
res                   = True                     ; plot mods desired
res@gsnMaximize         = True                   ; maximize plot in frame
res@xyMarkLineModes     = (/"Markers","Lines"/)  ; choose which have markers
res@xyMarkers           = 16                     ; choose type of marker
res@xyMarkerColor       = "red"                  ; Marker color
res@xyMarkerSizeF       = 0.005                ; Marker size (default 0.01)
res@xyDashPatterns      = 1                      ; solid line
res@xyLineThicknesses   = (/1,2/)                ; set second line to 2
res@xyComputeXMin       =True
res@trXMinF             =210
plot  = gsn_csm_xy(wks,x,data,res)                    ; create plot


画出来的图

画出来的图

所用文本数据

所用文本数据
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2018-1-12 08:53:50 | 显示全部楼层
小小蛋 发表于 2018-1-11 19:41
楼主好,ncl拟合二次曲线的问题解决了,我遇到了同样的情况,出现多条曲线,恳请赐教

http://www.ncl.ucar.edu/Applications/regress.shtml
你要的都在这
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

发表于 2017-4-9 18:15:08 | 显示全部楼层
关注一下
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

发表于 2018-1-11 19:41:34 | 显示全部楼层
楼主好,ncl拟合二次曲线的问题解决了,我遇到了同样的情况,出现多条曲线,恳请赐教
密码修改失败请联系微信:mofangbao
发表于 2018-6-13 09:45:20 | 显示全部楼层
Soaring 发表于 2018-1-12 08:53
http://www.ncl.ucar.edu/Applications/regress.shtml
你要的都在这

非常感谢~越来越爱这个平台~一起加油  
密码修改失败请联系微信:mofangbao
发表于 2018-6-13 09:45:23 | 显示全部楼层
Soaring 发表于 2018-1-12 08:53
http://www.ncl.ucar.edu/Applications/regress.shtml
你要的都在这

非常感谢~越来越爱这个平台~一起加油  
密码修改失败请联系微信:mofangbao
发表于 2022-11-4 20:08:49 | 显示全部楼层
你好,我想问一下这个怎么进行显著性检验呢?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-11-5 09:26:08 | 显示全部楼层
MARZOLINI 发表于 2022-11-4 20:08
你好,我想问一下这个怎么进行显著性检验呢?

可以用rtest函数
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2022-11-5 10:24:45 | 显示全部楼层

好的谢谢,我试试
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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