matlab如何用网格截取计算三维图?已经有了一个截断中间段的三维图,想用100*100的网格(按坐标系)将三维图分开,分别计算网格中nan出现的概率(有三维图的矩阵),显示出nan出现概率最小的

来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/19 06:50:49
matlab如何用网格截取计算三维图?已经有了一个截断中间段的三维图,想用100*100的网格(按坐标系)将三维图分开,分别计算网格中nan出现的概率(有三维图的矩阵),显示出nan出现概率最小的ma

matlab如何用网格截取计算三维图?已经有了一个截断中间段的三维图,想用100*100的网格(按坐标系)将三维图分开,分别计算网格中nan出现的概率(有三维图的矩阵),显示出nan出现概率最小的
matlab如何用网格截取计算三维图?
已经有了一个截断中间段的三维图,想用100*100的网格(按坐标系)将三维图分开,分别计算网格中nan出现的概率(有三维图的矩阵),显示出nan出现概率最小的网格.

matlab如何用网格截取计算三维图?已经有了一个截断中间段的三维图,想用100*100的网格(按坐标系)将三维图分开,分别计算网格中nan出现的概率(有三维图的矩阵),显示出nan出现概率最小的
clear all
[X,Y] = meshgrid(linspace(-3,3));
Z = peaks(X,Y);%生成数据z
 
zR=[1,6];%设置z范围,你自己可以改
z1=Z.*(Z>=zR(1)&Z<=zR(2));%为取出设定范围数据而得到逻辑数组
z2=Z;%为不覆盖原始Z,新建变量
z2(z1(:)==0)=NaN;%设定范围外的值改为NAN,得到初始数据

s=size(z2);
s1=[5,10];%设置分块大小(按数据点),你自己可以改,但必须根据原始数据大小能够整除,如不能整除可以剪裁或添加数据
%如果是按坐标系,你需要计算100×100的坐标中有多少个数据点,再来更改上面的设置
s2=[s(1)/s1(1),s(2)/s1(2)];%计算分多少块
for i=1:s2(1);
    for j=1:s2(2);
        b=z2((i-1)*s1(1)+1:i*s1(1),(j-1)*s1(2)+1:j*s1(2));%分块数据
        a(i,j)=length(find(isnan(b)));%计算每个分块中NaN的个数
    end
end

for i=1:s2(1);
    for j=1:s2(2);
        c((i-1)*s1(1)+1:i*s1(1),(j-1)*s1(2)+1:j*s1(2))=a(i,j);%让每一分块赋值为该块的个数       
    end
end

z3=(c<=min(c(:)));%为取出概率最小的数据重新赋颜色而得到逻辑数组
z4=Z;%为不覆盖原始Z,新建变量
z4(z3(:)==1)=max(Z(:));%设定概率最小的数据颜色暗红色

figure(2)
mesh(X,Y,z2,Z);
% surf(X,Y,z2,a2);
axis([-3 3 -3 3 -10 10])
% shading interp
title('裁剪的原始数据');
figure(3)
mesh(X,Y,z2,c);
% surf(X,Y,z2,a2);
axis([-3 3 -3 3 -10 10])
% shading interp
title('裁剪的原始数据分块NaN个数的分布');
figure(4)
mesh(X,Y,z2,z4);
% surf(X,Y,z2,a2);
axis([-3 3 -3 3 -10 10])
% shading interp
title('裁剪的原始数据NaN最少的分块区域');