matlab中二维插值函数interp2的使用详解

作者:天糊土 时间:2023-08-11 00:28:45 

下面是一段产生log-normal分布的代码,以此进行说明。


clear all;
clc;
for t=1:100
Traffic(t) =curve(t);
end
MaxTraffic = max(Traffic);
w = 0.2;
Wmax = 2*pi*w/3000;
x=[0:10:300];
y=[0:10:300];
Nx=length(x);
Ny=length(y);
Sigma = 0.53;
t = 0;
M = 10*curve(t)/MaxTraffic;
sum = 0;
for i=1:Nx
forj=1:Ny
 Mu = log(M)-0.5*Sigma^2;
 Rho(i,j) = RhoFromCoordination(x(i),y(j),Wmax,Sigma,Mu);
 Lognrnd(i,j) = round(exp(Sigma*Rho(i,j)+Mu));
 sum = Lognrnd(i,j)+sum;
end
end
sum
[xi,yi]=meshgrid(0:2:300,0:2:300);
z1=interp2(x,y,Lognrnd,xi,yi,'spline');%三次样条插值
surf(xi,yi,z1)

(1)首先理解meshgrid的原理和用法。简单地说,就是产生Oxy平面的网格坐标。

在进行3-D绘图操作时,涉及到x、y、z三组数据,而x、y这两组数据可以看做是在Oxy平面内对坐标进行采样得到的坐标对(x,y)。例如,要在“3<=x<=5,6<=y<=9,z不限制区间”这个区域内绘制一个3-D图形,如果只需要整数坐标为采样点的话。

我们可能需要下面这样一个坐标构成的矩阵:

(3,9),(4,9),(5,9);

(3,8),(4,8),(5,8);

(3,7),(4,7),(5,7);

(3,6),(4,6),(5,6);

在matlab中我们可以这样描述这个坐标矩阵

把各个点的x坐标独立出来,得:

3,4,5;

3,4,5;

3,4,5;

3,4,5;

再把各个点的y坐标也独立出来:

9,9,9;

8,8,8;

7,7,7;

6,6,6;

这样对应的x、y结合,便表示了上面的坐标矩阵。meshgrid就是产生这样两个矩阵,来简化我们的操作。然后根据(x,y)计算获得z,并绘制出三维图形。

(2)理解interp2的参数含义和用法,如ZI= interp2(X,Y,Z,XI,YI,'spline')

A、返回矩阵ZI,ZI的元素包含对应于参量XI与YI(可以是向量、或同型矩阵)的元素, 即ZI(i,j)←(XI(i),YI(j))

B、用户可以输入行向量和列向量XI与YI。

C、若XI与YI中有在X与Y范围之外的点,则相应地返回nan(Not a Number)。

D、用指定的算法method计算二维插值:

'linear' :双线性插值算法(缺省算法);
'nearest' :最临近插值;
'spline' :三次样条插值;
'cubic' :双三次插值。

E、如以下的运用:

[xi,yi] = meshgrid(0:2:300,0:2:300);
ZI = interp2(x,y,Lognrnd,xi,yi,'spline');%三次样条插值
surf(xi,yi,ZI)%这里已经不再是(x,y),而是(xi,yi)。

(3)上述的代码效果

插值前:

matlab中二维插值函数interp2的使用详解

插值后:

matlab中二维插值函数interp2的使用详解

补充知识:Matlab 二维插值,求面积

matlab中二维插值函数interp2的使用详解

matlab中二维插值函数interp2的使用详解

先将表中数据复制到EXCEL中,再导入到MATLAB中

这里只做了前两问,第三位实在不会,等学会了再补

第二问本来想着用差分求出来导数,再用面积公式,结果发现连z=f(x,y)我都不会表示。。。。。。

直接用的海伦公式,每一个方块内分成两个三角形,分开求面积


x=0:100:1200;
y=0:100:1000;
[x y]=meshgrid(x,y);
z=mydata1;
x1=0:10:1200;
y1=0:10:1000;
[x1 y1]=meshgrid(x1,y1);
z1=interp2(x,y,z,x1,y1)
surf(x1,y1,z1)
shading flat
square=0;
[r c]=size(z1)
for n=1:c-1
for m=1:r-1
 a=x1(m,n+1)-x1(m,n);
 b=y1(m+1,n)-y1(m,n);
 temp=z1(m+1,n+1)-z1(m,n);
 c=sqrt(a*a+b*b+temp*temp);
 temp=z1(m,n+1)-z1(m,n);
 a1=sqrt(temp*temp+a*a);
 temp=z1(m+1,n+1)-z1(m,n+1);
 b1=sqrt(temp*temp+b*b);
 p=(a1+b1+c)/2;
 square=square+sqrt(p*(p-a1)*(p-b1)*(p-c));
 temp=z1(m+1,n)-z1(m,n);
 b2=sqrt(b*b+temp*temp);
 temp=z1(m+1,n+1)-z1(m+1,n);
 a2=sqrt(temp*temp+a*a);
 p=(a2+b2+c)/2;
 square=square+sqrt(p*(p-a2)*(p-b2)*(p-c));
end
end
square

square =

1.2210e+06

matlab中二维插值函数interp2的使用详解

来源:https://blog.csdn.net/oqqHuTu12345678/article/details/52664433

标签:matlab,二维插值,interp2
0
投稿

猜你喜欢

  • 部署.Net6项目到docker

    2024-06-05 15:43:46
  • perl比较两个文件字符串的实例代码

    2023-02-15 15:40:18
  • web前端vue之CSS过渡效果示例

    2024-04-10 10:33:20
  • 一文带你上手Vue新的状态管理Pinia

    2024-05-09 15:11:33
  • python利用元类和描述器实现ORM模型的详细步骤

    2023-11-13 14:54:12
  • Python Opencv实现最强美颜滤镜效果

    2021-11-09 11:23:15
  • 基于python的mysql复制工具详解

    2023-06-24 01:08:28
  • Python异常学习笔记

    2021-03-24 01:33:33
  • 一文教你用python编写Dijkstra算法进行机器人路径规划

    2021-09-16 21:09:28
  • Pytorch学习笔记DCGAN极简入门教程

    2022-05-28 17:29:02
  • DD_belatedPNG,解决IE6不支持PNG绝佳方案

    2009-04-27 16:30:00
  • oracle索引总结

    2024-01-17 05:08:37
  • python中几种自动微分库解析

    2022-12-20 18:24:02
  • Python 中random 库的详细使用

    2022-01-19 05:35:14
  • Python:Numpy 求平均向量的实例

    2021-09-25 08:01:15
  • 信息分类是为了更好的索引

    2010-02-04 17:51:00
  • MYSQL初学者使用指南[适用自己安装mysql者]

    2007-08-06 14:53:00
  • 分享整理的12条sql语句连同数据

    2012-07-11 16:14:59
  • vue之webpack -v报错解决方案总结

    2024-05-09 09:51:23
  • JavaScript高级程序设计 阅读笔记(十五) 浏览器中的JavaScript

    2024-04-10 10:57:11
  • asp之家 网络编程 m.aspxhome.com