网站建设 视频教程,太原网站优化多少钱,做一个个人网站的步骤,邢台网警1. 有关 Matlab 获取代码关注WZZHHH回复关键词#xff0c;或者咸鱼关注#xff1a;WZZHHH123 怀俄明探空站数据解算PWV和Tm#xff1a;怀俄明探空站数据解算PWV和Tm 怀俄明多线程下载探空站数据#xff08;包括检查和下载遗漏数据的代码#xff09;#xff1a;怀俄明多线…1. 有关 Matlab 获取代码关注WZZHHH回复关键词或者咸鱼关注WZZHHH123 怀俄明探空站数据解算PWV和Tm怀俄明探空站数据解算PWV和Tm 怀俄明多线程下载探空站数据包括检查和下载遗漏数据的代码怀俄明多线程下载 对IGRAv2进行质量控制得到PWV和Tm的 matlab 代码IGRAv2进行质量控制得到PWV和Tm 算 IGRAv2 探空站的 Tm 和 PWVIGRAv2计算Tm和PWV 提取探空站 IGRAv2 全部数据 ERA5 解算合集温度、气压、PWV、水汽压和 Tm代码获取ERA5合集
MATLAB画散点密度图 代码、测试数据的压缩包在文末可自行下载测试出图效果
散点误差图出图效果 代码如下
具体使用请自行修改
function scatter_plot(x,y,derror)
% 点状密度图
% 根据测试数据画的实际应用自己修改
% --------------------改变密度大小------------------
numbins 50;
% --------------------------------------------------X x;
Y y;
% 找到最大值最小值
MaxSPmax(x);MaxVmax(y);Maxi1.1*max(MaxSP,MaxV);[values, centers] hist3([X Y], [numbins numbins]);
centers_X centers{1,1};
centers_Y centers{1,2};
binsize_X abs(centers_X(2) - centers_X(1)) / 2;
binsize_Y abs(centers_Y(2) - centers_Y(1)) / 2;
bins_X zeros(numbins, 2);
bins_Y zeros(numbins, 2);
for i 1:numbinsbins_X(i, 1) centers_X(i) - binsize_X;bins_X(i, 2) centers_X(i) binsize_X;bins_Y(i, 1) centers_Y(i) - binsize_Y;bins_Y(i, 2) centers_Y(i) binsize_Y;
end
scatter_COL zeros(length(X), 1);
onepercent round(length(X) / 100);
for i 1:length(X)if (mod(i,onepercent) 0)fprintf(.);endlast_higher_X NaN;id_X NaN;c_X X(i);last_lower_X find(c_X bins_X(:,1));if (~isempty(last_lower_X))last_lower_X last_lower_X(end);elselast_higher_X find(c_X bins_X(:,2));if (~isempty(last_higher_X))last_higher_X last_higher_X(1);endendif (~isnan(last_lower_X))id_X last_lower_X;elseif (~isnan(last_higher_X))id_X last_higher_X;endendlast_higher_Y NaN;id_Y NaN;c_Y Y(i);last_lower_Y find(c_Y bins_Y(:,1));if (~isempty(last_lower_Y))last_lower_Y last_lower_Y(end);elselast_higher_Y find(c_Y bins_Y(:,2));if (~isempty(last_higher_Y))last_higher_Y last_higher_Y(1);endendif (~isnan(last_lower_Y))id_Y last_lower_Y;elseif (~isnan(last_higher_Y))id_Y last_higher_Y;endendscatter_COL(i) values(id_X, id_Y);
end% 拟合曲线
apolyfit(x,y,1);
axlinspace(0,Maxi10,2000);
aya(1)*axa(2);Str_yx[Y ,num2str(round(a(1),2)),*X ,num2str(round(a(2),2))];scatter(x, y, 20, scatter_COL, . );
hold oncolormap(jet);
colorbar;xlabel(GNSS-PWV/mm,FontSize,12,FontName,Times New Roman);
ylabel(ERA5-PWV/mm,FontSize,12,FontName,Times New Roman);bx linspace(0,120,2000);% 设置xy轴的界限
% xlim([0 120]);ylim([0 120]);% 把误差表示在图上位置根据自己需要改我这是根据测试数据填的
text(5,115,Str_yx,FontSize,12,FontName,Times New Roman)
text(5,107,[Bias ,num2str(round(100*derror(1))/100,%.2f)],FontSize,12,FontName,Times New Roman)
text(5,99,[RMSE ,num2str(round(100*derror(2))/100,%.2f)],FontSize,12,FontName,Times New Roman)
text(5,91,[Cor ,num2str(derror(3),%.2f)],FontSize,12,FontName,Times New Roman)% 画拟合曲线
plot(ax,ay,r)
hold on% 画对角线
plot(bx,bx,k);
hold ontitle(测试数据,FontSize,12)
box on
grid onset(gca,LineWidth,1.2)
end压缩包地址如下 压缩包里面的内容 压缩包地址 关注: WZZHHH回复 MATLAB画散点密度图
. 有问题请你问你的老师师兄师姐。我也挺忙的。
