01方程的几何解释

一、n个未知数,n个方程

1.1 两个未知数,两个方程

{ 2 x − y = 0 − x + 2 y = 3 (1) \left\{ \begin{aligned} &2x-y=0\\ &-x+2y=3\\ \end{aligned} \right.\tag{1} {2xy=0x+2y=3(1)
矩阵形式为:

[ 2 − 1 − 1 2 ] [ x y ] = [ 0 3 ] (2) \begin{bmatrix}2& -1 \\-1& 2 \\ \end{bmatrix} \begin{bmatrix}x\\y\\ \end{bmatrix} =\begin{bmatrix}0\\3\\ \end{bmatrix}\tag{2} [2112][xy]=[03](2)
记作: A x = b Ax=b Ax=b

  • 行图像
    行图像就是分别作出两条直线的图像,根据交点求出方程的解。如图, B ( 1 , 2 ) B(1,2) B(1,2)为方程组(1)对应图像的交点。
    在这里插入图片描述

  • 列图像
    即将两个方程组看成是列向量的线性组合:
    x [ 2 − 1 ] + y [ − 1 2 ] = [ 0 3 ] (3) x\begin{bmatrix}2\\-1\\ \end{bmatrix} +y\begin{bmatrix}-1\\2 \end{bmatrix} =\begin{bmatrix}0\\3\\ \end{bmatrix}\tag{3} x[21]+y[12]=[03](3)
    在这里插入图片描述

如上图,找到一组系数 ( x , y ) (x,y) (x,y)使得图像其组合恰好等于右向量 b = [ 0 3 ] b=\begin{bmatrix}0\\3\\ \end{bmatrix} b=[03]。此例中,对于任意的系数组合,我们将会得到很多个右向量组合,这将会遍布整个平面。

1.2 三个未知数,三个方程

{ 2 x − y = 0 − x + 2 y − z = − 1 − 3 y + 4 z = 4 (4) \left\{ \begin{aligned} &2x-y=0\\ &-x+2y-z=-1\\ &-3y+4z=4 \end{aligned} \right.\tag{4} 2xy=0x+2yz=13y+4z=4(4)
矩阵形式为:
[ 2 − 1 0 − 1 2 − 1 0 − 3 4 ] [ x y z ] = [ 0 − 1 4 ] (5) \begin{bmatrix}2& -1&0 \\-1& 2&-1 \\0&-3& 4\end{bmatrix} \begin{bmatrix}x\\y\\z \end{bmatrix} =\begin{bmatrix}0\\-1\\4 \end{bmatrix}\tag{5} 210123014 xyz = 014 (5)
同样记作 A x = b Ax=b Ax=b

  • 行图像

依次画出三个平面的图像,交点即为方程的解。

  • 列图像

列向量的线性组合系数即为方程的解。


【20220728】原来的图太大了,改小了些

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是基于MATLAB的R-D方程几何纠正代码: ```matlab % 读取SAR图像 img = imread('sar_image.tif'); % 定义SAR参数 fc = 5.3e9; % 雷达频率 lambda = 3e8/fc; % 波长 R0 = 3000; % 雷达距离 PRF = 5000; % 脉冲重复频率 V = 700; % 雷达平台速度 h = 1000; % 雷达平台高度 theta_s = 20; % 俯仰角 theta_i = 30; % 方位角 % 定义地球参数 Re = 6371000; % 地球平均半径 ecc = 0.08181919; % 地球椭球体偏心率 Rp = Re*(1-ecc^2); % 地球极半径 % 计算SAR图像像素坐标系 [N, M] = size(img); x = linspace(-M/2, M/2, M); y = linspace(-N/2, N/2, N); [X, Y] = meshgrid(x, y); % 计算地球坐标系 theta = atan2(Y, R0+X); R = sqrt((R0+X).^2 + Y.^2); phi = asin(cos(theta_s)*sin(theta_i) + sin(theta_s)*cos(theta_i).*cos(theta)); lambda_s = atan2(sin(theta_i).*sin(theta), cos(theta_i).*cos(theta_s)-sin(theta_i).*sin(theta_s).*cos(theta)); phi_s = asin(sin(theta_s).*sin(theta_i) - cos(theta_s).*cos(theta_i).*cos(theta)); R_s = h./cos(phi_s); % 计算R-D方程参数 a = V^2/(2*fc*Rp*cos(phi_s)); b = V^2/(2*fc*Rp*cos(phi)); c = V^2/(2*fc*Rp); % 计算几何校正后的R-D图像 Rc = sqrt(R.^2 + (a-b)^2 - 2*R.*(a-b).*cos(phi)); D = sqrt((R_s-Rc).^2 + c^2 - 2*(R_s-Rc).*c.*cos(phi_s-phi)); img_corr = img.*(R./Rc).^2.*cos(phi)./D; % 显示几何校正后的R-D图像 figure; imshow(abs(img_corr), []); ``` 需要注意的是,这个代码只是一个简单的示例,实际应用中需要根据具体的SAR参数和地球参数进行调整。此外,由于R-D方程几何校正需要用到DEM数据,因此在实际应用中还需要考虑如何获取和处理DEM数据。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值