【代码】矩阵的线性映射(平面转换)

相同x、y坐标下,改变z坐标

如下图所示,希望将原矩阵由上面的平面转换到下面的平面

已知上面平面所有点的三维坐标和下面平面的角点(四个顶点)的三维坐标

代码:

matlab 复制代码
function main()
%矩阵线性变换,相同x、y,仅z坐标不同的两个平面相互转换
clc
clear
close all
%% 0.原始数据、目标数据
x=[1,2,3,4,5];
y=[1,2,3]*10;
%原始数据
res_ori=[1,2,3,4,5;
    2,3,4,5,6;
    3,4,5,6,7];
%目标数据
% res_tar=[5,4,3,2,1;
%     6,5,4,3,2;
%     7,6,5,4,3]*-1;

%% 1.目标数据的角点
r11=-5;
r12=-1;
r21=-7;
r22=-3;
mat_tar_cor=[r11,r12;r21,r22];%目标矩阵(平面)的四个角点
%% 2.目标数据
results=trans_func(res_ori,mat_tar_cor);

%% 3.画图
[X,Y]=meshgrid(x,y); 
mesh(X,Y,res_ori)
hold on
mesh(X,Y,results)

trans_func.m

matlab 复制代码
function results=trans_func(mat_ori,mat_tar)
%线性映射
rrr=mat_ori;
qwe=mat_tar;
q1=length(rrr(1,:));%列数
q2=length(rrr(:,1));%行数
for i=1:q1
    qweqwe(1,i)=(qwe(1,2)-qwe(1,1))/(rrr(1,q1)-rrr(1,1))*(rrr(1,i)-rrr(1,1))+qwe(1,1);
end
for i=1:q1
    qweqwe(q2,i)=(qwe(2,2)-qwe(2,1))/(rrr(1,q1)-rrr(1,1))*(rrr(1,i)-rrr(1,1))+qwe(2,1);
end

for j=1:q1%列
    for i=1:q2%行
        qweqwe(i,j)=(qweqwe(q2,j)-qweqwe(1,j))/(rrr(q2,j)-rrr(1,j))*(rrr(i,j)-rrr(1,j))+qweqwe(1,j);
    end
end
results=qweqwe;
相关推荐
豆沙沙包?1 小时前
2025年--Lc201- 378. 有序矩阵中第 K 小的元素(排序)--Java版
java·线性代数·矩阵
CLubiy3 小时前
【研究生随笔】Pytorch中的线性代数
pytorch·python·深度学习·线性代数·机器学习
_码力全开_6 小时前
P1005 [NOIP 2007 提高组] 矩阵取数游戏
java·c语言·c++·python·算法·矩阵·go
张晓~183399481219 小时前
碰一碰发视频 系统源码 /PHP 语言开发方案
开发语言·线性代数·矩阵·aigc·php·音视频·文心一言
dxnb221 天前
Datawhale25年10月组队学习:math for AI+Task2线性代数
人工智能·学习·线性代数
豆沙沙包?1 天前
2025年--Lc187--120. 三角形最小路径和(多维动态规划,矩阵)--Java版
java·矩阵·动态规划
豆沙沙包?2 天前
2025年--Lc186--64. 最小路径和(多维动态规划,矩阵)--Java版
java·矩阵·动态规划
lingchen19062 天前
矩阵的除法
人工智能·算法·矩阵
鼓掌MVP2 天前
图生3D技术解析:从二维平面到立体世界的智能飞跃
平面·3d
老歌老听老掉牙2 天前
使用 OpenCASCADE 提取布尔运算后平面图形的外轮廓
c++·平面·opencascade