MATLAB的答案

我怎么找到这两个功能的拦截i1和i2

1视图(30天)
亚历山大•辛格
亚历山大•辛格 2021年10月10日
评论道: 明星黾 2021年10月12日
% hemuppgift 3
一个= 0.5;
ih = 10 ^ 7;
工业区= 10 ^ 5;
F = 96485;
R = 8.134;
T = 298.15;
h = 0.0001;
x = 1: h: 1;
z1 = -1.87: h: 0.13;
z2 = -2.32: h: -0.32;
%
% E1 = @ (z) (x - 0.87);
% E2 = @ (z) (x - 1.37);
i1 = @ (x) (ih * (exp ((* F * x) / (R * T))——exp (F (——* * x) / (R * T))));
i2 = @ (x)(工业区* (exp ((* F * x) / (R * T))——exp (F (——* * x) / (R * T))));
logh = log10 (i1 (x));
logz = log10 (i2 (x));
情节(logh z1,“b”)
持有在
情节(logz z2,“r”)

在回答这个问题。

接受的答案

明星黾
明星黾 2021年10月10日
跑:
这是一个挑战!
之前我试过这几次有适当的洞察力计算交点。(我做了一些小的改变代码,具体变化 “logh” 和 “logz” 为清晰的匿名函数,但是除了我的添加,以便于计算交点,这都是commente-documented,另有不变。)
只有一个真正的交集, (-4.06,-1.21) 。
% hemuppgift 3
一个= 0.5;
ih = 10 ^ 7;
工业区= 10 ^ 5;
F = 96485;
R = 8.134;
T = 298.15;
h = 0.0001;
x = 1: h: 1;
z1 = -1.87: h: 0.13;
z2 = -2.32: h: -0.32;
%
% E1 = @ (z) (x - 0.87);
% E2 = @ (z) (x - 1.37);
i1 = @ (x) (ih * (exp ((* F * x) / (R * T))——exp (F (——* * x) / (R * T))));
i2 = @ (x)(工业区* (exp ((* F * x) / (R * T))——exp (F (——* * x) / (R * T))));
logh = @ (x) log10 (i1 (x));
logz = @ (x) log10 (i2 (x));
common_z = linspace (min ([z1, z2],[], 2)、马克斯([z1, z2], [], 2), 1 e + 3);%常见的“z”向量(由“z1”&“z2”)
common_logh = interp1 (z1,实际(logh (x)), common_z);%的logh插入“common_z”
common_logz = interp1 (z2,实际(logz (x)), common_z);%的logz插入“common_z”
common_dif = common_logh - common_logz;%的差异
Lidx = ~ isnan (common_dif);%去除“NaN”值
z_intx = interp1 (common_dif (Lidx) common_z (Lidx), 0)%计算“z”的十字路口
z_intx = -1.2104
logh_intx = interp1 (z1, logh (x) z_intx)%计算logh的十字路口(也可以使用“logz”)
我logh_intx = -4.0590 + 1.3644
图
情节(真实(logh (x)), z1,“b”)
持有在
情节(真实(logz (x)), z2,“r”)
yline (z_intx“g”)
参照线(真实(logh_intx),“g”)
传奇(“logh”,“logz”,交叉线的,“位置”,“最佳”)
图
情节(x,真实(logh (x)),“b”)
持有在
情节(x,图像放大(logh (x)),“b”)
情节(x,真实(logz (x)),“r”)
情节(x,图像放大(logz (x)),“——r”)
网格
标题(绘制功能的x”)
实验,得到不同的结果。
。

更多的答案(0)

类别

找到更多的在字符和字符串在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

翻译的