ncl绘制直方图和概率密度图

礼乐王

这种图谁会呀  求大神帮忙

尽头的尽头

http://www.ncl.ucar.edu/Applications/pdf.shtml
NCL有概率函数

礼乐王

尽头的尽头 发表于 2017-3-31 08:09
http://www.ncl.ucar.edu/Applications/pdf.shtml
NCL有概率函数

我想两个放一起的，难不成是要叠加吗

几米的猜想

origin可以实现吧

尽头的尽头

礼乐王 发表于 2017-3-31 11:32
我想两个放一起的，难不成是要叠加吗

叠加也可以，xy也可以

礼乐王

几米的猜想 发表于 2017-3-31 11:40
origin可以实现吧

origin??一个软件吗 还是ncl的一个函数，兄弟，能说具体点不

礼乐王

礼乐王 发表于 2017-3-31 14:41
origin??一个软件吗 还是ncl的一个函数，兄弟，能说具体点不

begin
dev   = asciiread("./dev.txt",-1,"float")
opt  = True
opt@bin_min = -20.
opt@bin_max =  20.
a    =  pdfx(dev,0,opt)

wks  = gsn_open_wks ("X11","PDF")  
res  = True
res@xyLineThicknessF       = 2.
res@xyLineColors           = "blue"
res@xyDashPatterns            = "dash"
plot = gsn_csm_xy (wks,a@bin_center,a,res)

res1                          = True
res1@gsnHistogramComputePercentages  = True
res1@gsnHistogramPercentSign        = True
res1@tiYAxisString          = "number"
res1@tiXAxisString          = "error   (m/s)"
res1@gsnHistogramBarWidthPercent = 90.
res1@tiMainString                = "Bar width is 80% of bin width"
res1@tmXBLabelStride          = 10
res1@gsnHistogramBinIntervals = ispan(-20,20,1)
mplot = gsn_histogram(wks,dev,res1)

overlay(mplot,plot)
draw(mplot)
frame(wks)
end

兄弟这是脚本，你看一下哪的问题，叠加没效果呀，出不来那种图

礼乐王

begin
dev   = asciiread("./dev.txt",-1,"float")
opt  = True
opt@bin_min = -20.
opt@bin_max =  20.
a    =  pdfx(dev,0,opt)

wks  = gsn_open_wks ("X11","PDF")  
res  = True
res@xyLineThicknessF       = 2.
res@xyLineColors           = "blue"
res@xyDashPatterns            = "dash"
plot = gsn_csm_xy (wks,a@bin_center,a,res)

res1                          = True
res1@gsnHistogramComputePercentages  = True
res1@gsnHistogramPercentSign        = True
res1@tiYAxisString          = "number"
res1@tiXAxisString          = "error   (m/s)"
res1@gsnHistogramBarWidthPercent = 90.
res1@tiMainString                = "Bar width is 80% of bin width"
res1@tmXBLabelStride          = 10
res1@gsnHistogramBinIntervals = ispan(-20,20,1)
mplot = gsn_histogram(wks,dev,res1)

overlay(mplot,plot)
draw(mplot)
frame(wks)
end
这是脚本  ，大神帮忙看一下问题所在

月满西楼

礼乐王 发表于 2017-3-31 14:41
origin??一个软件吗 还是ncl的一个函数，兄弟，能说具体点不

origin是一个软件。。。

礼乐王

月满西楼 发表于 2017-3-31 21:58
origin是一个软件。。。

想用ncl呀，想精通一种软件