Matlab联立隐函数H(x,y1)和M(x,y2)求解新函数C(y1,y2)和因变量x的关系

目录

1. 问题描述

2. 函数图像

2.1 隐函数H的图像

2.2 隐函数M的图像

3. 解决思路

3.1 初步思路

3.2 实际情况

3.3 解决方法

4. 完整代码

5. 交流讨论


1. 问题描述

已知条件如下:

①隐函数H(x,y1) = exp(2*x-y1)*log(x-3*y1)-4 = 0

②隐函数M(x,y2) = log(2*x-y2)/exp(x-2*y2)+3 = 0

现有函数y =  C(y1,y2) = exp(1-y1*y2),求y和x在坐标轴范围x∈[2,10],y∈[0,25]的函数图像


2. 函数图像

我们借助matlab的官方函数fimplicit绘制隐函数H和M的图像。

fimplicit函数接收两个参数,一个是匿名函数,一个是坐标轴范围。

2.1 隐函数H的图像

figure
H = @(x,y) exp(2*x-y)*log(x-3*y)-4;
bound = [2,10,0,25];
h = fimplicit(H, bound);
title('隐函数H')

Matlab联立隐函数H(x,y1)和M(x,y2)求解新函数C(y1,y2)和因变量x的关系

2.2 隐函数M的图像

%% 绘制隐函数M
figure
M = @(x,y) log(2*x-y)/exp(x-2*y)+3;
bound = [2,10,0,25];
m = fimplicit(M, bound);
title('隐函数M')

 Matlab联立隐函数H(x,y1)和M(x,y2)求解新函数C(y1,y2)和因变量x的关系

3. 解决思路

3.1 初步思路

  1. 对于隐函数H(x1,y1)的图像,可以得到一组一一对应的向量x1和y1,同理对于隐函数M(x2,y2),也可以得到一组一一对应的向量x2和y2;
  2. 如果这里的x1和x2完美契合,即x1向量和x2向量完全一致,不妨用x来表示x1和x2,那么只需要将y1向量和y2向量代入函数表达式y =  C(y1,y2) = exp(1-y1*y2),计算得到y向量,利用plot(x, y)就可以得到目标函数图像。

3.2 实际情况

对于隐函数H,函数fimplicit不光绘制隐函数H在bound上的图像,还会返回该图像窗口,这里用h接收其返回的图像窗口

h = fimplicit(H, bound);

在工作区内查看该h变量,可以发现其有子属性XData和YData,分别对应图像上各点的横纵坐标。

Matlab联立隐函数H(x,y1)和M(x,y2)求解新函数C(y1,y2)和因变量x的关系

 于是,利用fimplicit的返回值,可以得到隐函数H和M图像的x、y向量

%% 获取隐函数H的x1、y1向量
x1 = h.XData;
y1 = h.YData;
%% 获取隐函数M的x2、y2向量
x2 = m.XData;
y2 = m.YData;

遗憾的是,x1和x2并没有完全一致(完全一致的要求是x1和x2的长度首先需要一样,其次两个向量相减为0向量)。

Matlab联立隐函数H(x,y1)和M(x,y2)求解新函数C(y1,y2)和因变量x的关系

3.3 解决方法

对于这种情况,我们可以考虑对H的图像进行插值,人为地选取向量x = [x1, x2, ..., xn],并得到与其对应的yp1;同理,同样使用向量x = [x1, x2, ..., xn]对M的图像进行插值,得到与其对应的yp2,这样就建立了x与yp1、x与yp2的映射关系。

%% 设置x的插值点,对y1和y2进行插值
x = 2:0.5:10;             %设置x的插值点
yp1 = interp1(x1, y1, x); %用x对y1重新插值得到yp1
yp2 = interp1(x2, y2, x); %用x对y2重新插值得到yp2

进一步地,通过y = C(y1,y2) = exp(1-y1*y2),可以得到x与y的映射关系,从而绘制目标图像。

%% 计算及绘图
y = exp(1-yp1.*yp2);
figure
plot(x, y), title('x和y的图像')

4. 完整代码

clear, close all
%% 绘制隐函数H
figure
H = @(x,y) exp(2*x-y)*log(x-3*y)-4;
bound = [2,10,0,25];
h = fimplicit(H, bound);
title('隐函数H')
%% 绘制隐函数M
figure
M = @(x,y) log(2*x-y)/exp(x-2*y)+3;
bound = [2,10,0,25];
m = fimplicit(M, bound);
title('隐函数M')
%% 获取隐函数H的x1、y1向量
x1 = h.XData;
y1 = h.YData;
%% 获取隐函数M的x2、y2向量
x2 = m.XData;
y2 = m.YData;
%% 设置x的插值点,对y1和y2进行插值
x = 2:0.5:10;             %设置x的插值点
yp1 = interp1(x1, y1, x); %用x对y1重新插值得到yp1
yp2 = interp1(x2, y2, x); %用x对y2重新插值得到yp2
%% 计算及绘图
y = exp(1-yp1.*yp2);
figure
plot(x, y), title('x和y的图像')

5. 交流讨论

Matlab联立隐函数H(x,y1)和M(x,y2)求解新函数C(y1,y2)和因变量x的关系

上一篇:appium页面滑动封装思路


下一篇:时区概念