解不了图片中要求得值，请教大神，非常感谢！

busahuang

function [TD] = menxian2(L)
) c* l4 j- M1 V9 C( G- j  ]" ]%UNTITLED4 此处显示有关此函数的摘要
+ L; a" u* s) y2 U( \4 u# p: @%   此处显示详细说明
# l1 R# R$ b" O- d! O+ g: gsyms s y x t
s=int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
y=matlabFunction(s);
f1=finverse(s,t);
5 [4 B2 w$ R6 D  df2=matlabFunction(f1);
TD=f1(1-L);
. C/ e9 Y7 X' d6 G$ a6 b& K) Xend

2 p: c' T, E( Q# g, G代码如上，解不了图片中要求得值，请教大神，非常感谢！

7 [8 z- Q/ R: R关于这个函数得解，Pfa已知，求Td
9 M' ?. n# G: X9 x+ S8 A1 h% K# z


- g1 Q: ]6 H7 B

Ferrya

: b7 B$ z% ]5 g% c% j直接用TD = norminv( 1- L );就行了
6 s& D4 s# w9 f1 i, y非要解方程的话用vpasolve求数值解
syms x t real
s = int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
9 l6 K) D8 E9 }# k$ x' j- e- ^y = matlabFunction(s);
6 e- ~6 W( `2 z% ?; f4 BTD = vpasolve( y(x)==1-L );

busahuang

Ferrya 发表于 2020-5-28 15:03
# `& X# s, S( `+ b+ N( h" F7 o直接用TD = norminv( 1- L );就行了
非要解方程的话用vpasolve求数值解
4 N5 [: E- e! jsyms x t real

大佬看看求这个Td用的代码跟前面这个一样吗，n已知。
- b4 X4 c* K0 v1 g( C6 b' Y! c

# B8 H2 e9 D. G; K. M