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

busahuang

function [TD] = menxian2(L)
%UNTITLED4 此处显示有关此函数的摘要
%   此处显示详细说明
4 `7 O- o% f/ xsyms s y x t
2 k) s/ L3 D6 M( Ks=int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
. v8 W" ?" d* w- ~. jy=matlabFunction(s);
. o1 w. n1 P& o* G0 Mf1=finverse(s,t);
: ^# r" k; u* K  z+ \f2=matlabFunction(f1);
! R; }- `& W6 d; A1 ?TD=f1(1-L);
1 w9 V! Y- P- [1 t8 D1 A, vend

( }. n5 q* n( h, U代码如上，解不了图片中要求得值，请教大神，非常感谢！
8 W+ c7 b4 N; U
关于这个函数得解，Pfa已知，求Td
+ V4 S4 [3 W2 Y& i7 [; ]; ]


, u. G9 s2 d5 J: f
- o; z+ }0 S2 s% _( t- L+ b- e% @2 d

Ferrya

' h% |! e1 o7 @) O直接用TD = norminv( 1- L );就行了
非要解方程的话用vpasolve求数值解
  h7 `. h6 o2 Ssyms x t real
s = int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
% |4 W7 l! e' jy = matlabFunction(s);
. w" Z& f1 ^9 f( @9 `TD = vpasolve( y(x)==1-L );

busahuang

Ferrya 发表于 2020-5-28 15:03
直接用TD = norminv( 1- L );就行了
非要解方程的话用vpasolve求数值解
. t+ ?" b1 q; K: S. qsyms x t real

大佬看看求这个Td用的代码跟前面这个一样吗，n已知。

) m1 |1 L! V: ]& c5 w