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

busahuang

function [TD] = menxian2(L)
+ W( L+ O7 f4 T. f" g+ k%UNTITLED4 此处显示有关此函数的摘要
%   此处显示详细说明
, o0 j( _+ S0 @. [7 Msyms s y x t
* y( Y' x! J) G8 |: {% Vs=int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
+ J) G' \& \8 o) U) iy=matlabFunction(s);
f1=finverse(s,t);
f2=matlabFunction(f1);
. d7 j1 _% `0 z# b# h+ l  ?( _TD=f1(1-L);
end
$ l  W9 ^  F" f3 g; ]% j5 D5 n2 r
, L/ y8 _. Z* M6 K* b& X$ @代码如上，解不了图片中要求得值，请教大神，非常感谢！
4 X! q5 T' x' X% W5 J
关于这个函数得解，Pfa已知，求Td



* M: c. @. G. l; i* ~
2 p" C9 l' ^, p' M$ Q! B

Ferrya

直接用TD = norminv( 1- L );就行了
# m6 ]% M' |0 O( e3 A) }! f5 U* ~# q非要解方程的话用vpasolve求数值解
syms x t real
' Y. q% N9 X3 b# b* i. r' u  zs = int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
8 b. i: b( i* G  z& R1 ^5 Z5 ?y = matlabFunction(s);
4 b0 K) R0 o- kTD = vpasolve( y(x)==1-L );

busahuang

Ferrya 发表于 2020-5-28 15:03
# A# E7 V9 {7 X+ t( U3 L直接用TD = norminv( 1- L );就行了
5 e% y! n& a7 X) ?2 B1 k; {非要解方程的话用vpasolve求数值解
7 p% d) n& Z8 U# N" o  ysyms x t real

! }2 m$ a* {  a1 v大佬看看求这个Td用的代码跟前面这个一样吗，n已知。

! `4 M4 d" T4 J7 O: T; ^1 h3 C
) C2 J# f- }. {9 \& K, O8 J