MATLAB Answers

How to find the inverse of a function numerically

92 ビュー (過去 30 日間)
BeeTiaw
BeeTiaw 2019 年 1 月 28 日
コメント済み: BeeTiaw 2019 年 1 月 28 日
Hi expert,
Can someone tell me how is it possible to find the inverse of this function
in which z is a complex number and cannot be zero. , a and b are constant.
How to solve this function for ?
Data (Example)
The following is obtained using the following input:
for
theta=[
0
0.6981
1.3963
2.0944
2.7925
3.4907
4.1888
4.8869
5.5851
6.2832];
The results is:
z = [
2.9400 + 0.0000i
3.2277 + 2.1618i
2.2730 + 1.2986i
0.4557 - 2.2605i
-1.4094 - 6.6606i
-2.3950 - 9.7020i
-1.9857 -10.1025i
-0.4083 - 7.8642i
1.5321 - 3.9596i
2.9400 - 0.0000i
];
Example:
I used Matlab function "roots" to solve the following inversion problem
by rearrange the function into:
and use to function "roots" to find the solution.
  9 件のコメント
BeeTiaw
BeeTiaw 2019 年 1 月 28 日
Oh probably I can do it by multiplying them with

サインインしてコメントする。

採用された回答

John D'Errico
John D'Errico 2019 年 1 月 28 日
Multiply by zeta^2, and collect terms. As long as zeta is not zero, that is not a problem. Your equation reduces to
b*m2 + (a + b*m1)*zeta - z*zeta^2 + (a*m1 + b)*zeta^3 + (a*m2)*zeta^4 == 0
We only need to worry about zeta==0 if either of b or m2 was zero. In that case, zeta==0 would be one of the roots of the above equation.
Of the coefficients of the above equation, all are apparently known, and have fixed values. So there are 4 roots.
Then the "inverse" is given as any of the 4 roots of that equation, thus:
zetaroots = solve(b*m2 + (a + b*m1)*zeta - z*zeta^2 + (a*m1 + b)*zeta^3 + (a*m2)*zeta^4,zeta,'maxdegree',4);
You don't want me to write the entire expression in here, as it is a massive mess of terms.
The problem is, the "inverse" is a rather nasty mess of a function of z. There are 4 solutions. Even if I show only 5 digit numbers in that expression for all coefficients, it is still a nasty mess.
vpa(expand(subs(zetaroots,{a,b,m1,m2},[-2.0800,4.0800,0.5,-0.03])),5)
ans =
- (0.16667*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2))/(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/6) - (0.16667*(10680.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 70.15*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 256.82*z^2*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) + 2868.6*z*(5.1962*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 6767.6*z - 8231.4*z^3 - 125699.0)^(1/2) + 106211.0*(5.1962*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 6767.6*z - 8231.4*z^3 - 125699.0)^(1/2) + 192.31*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2))^(1/2))/((0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/6)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/4)) - 12.179
(0.16667*(10680.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 70.15*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 256.82*z^2*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) + 2868.6*z*(5.1962*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 6767.6*z - 8231.4*z^3 - 125699.0)^(1/2) + 106211.0*(5.1962*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 6767.6*z - 8231.4*z^3 - 125699.0)^(1/2) + 192.31*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2))^(1/2))/((0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/6)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/4)) - (0.16667*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2))/(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/6) - 12.179
(0.16667*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2))/(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/6) - (0.16667*(10680.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 70.15*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 256.82*z^2*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 2868.6*z*(5.1962*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 6767.6*z - 8231.4*z^3 - 125699.0)^(1/2) - 106211.0*(5.1962*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 6767.6*z - 8231.4*z^3 - 125699.0)^(1/2) + 192.31*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2))^(1/2))/((0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/6)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/4)) - 12.179
(0.16667*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2))/(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/6) + (0.16667*(10680.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 70.15*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 256.82*z^2*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2) - 2868.6*z*(5.1962*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 6767.6*z - 8231.4*z^3 - 125699.0)^(1/2) - 106211.0*(5.1962*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 6767.6*z - 8231.4*z^3 - 125699.0)^(1/2) + 192.31*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/2))^(1/2))/((0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/6)*(96.154*z*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 5340.2*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(1/3) + 9.0*(0.096225*(2.07e6*z^4 + 7.6638e7*z^3 + 1.1346e6*z^2 + 6.3008e7*z + 5.8506e8)^(1/2) - 125.33*z - 152.43*z^3 - 2327.6)^(2/3) + 256.82*z^2 + 70.15)^(1/4)) - 12.179
I'm not at all sure what you expected the inverse of your function would look like. But it is not pretty.
  2 件のコメント
BeeTiaw
BeeTiaw 2019 年 1 月 28 日
An answer for a much more generalised form of function is available here https://uk.mathworks.com/matlabcentral/answers/441867-tthe-inverse-of-a-function-numerically-with-n-terms

サインインしてコメントする。

その他の回答 (0 件)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by