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

busahuang

function [TD] = menxian2(L)
%UNTITLED4 此处显示有关此函数的摘要
1 P, X7 W" d1 E7 C4 e%   此处显示详细说明
8 w! r1 j6 T, w  M4 Ysyms s y x t
s=int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
y=matlabFunction(s);
" @: E3 o' O# O4 m3 p/ af1=finverse(s,t);
f2=matlabFunction(f1);
# x9 Q) A9 b) H5 N; E! H' g& q- `TD=f1(1-L);
! C0 R& O5 y  E& O% ?9 C  ?end

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

关于这个函数得解，Pfa已知，求Td

9 Y# c# ?4 R8 Z- o1 Q" a8 ]( q+ _

Ferrya

直接用TD = norminv( 1- L );就行了
) A6 S1 {$ p" X' S! f; E非要解方程的话用vpasolve求数值解
. `: l6 x2 f5 D$ s' d0 P- P# Ksyms x t real
s = int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
& g' p+ d1 C: [" f9 d7 {y = matlabFunction(s);
0 ?0 S+ {7 D! L! ^0 Q6 h5 F  YTD = vpasolve( y(x)==1-L );

busahuang

Ferrya 发表于 2020-5-28 15:03
; M2 q1 t9 K$ B2 ~4 o直接用TD = norminv( 1- L );就行了
; Y. q: t5 v* O- H8 [非要解方程的话用vpasolve求数值解
! y: ~# w7 s8 s4 N+ Vsyms x t real

) f: u' W' J) _: v% g9 z大佬看看求这个Td用的代码跟前面这个一样吗，n已知。


. f8 _# z  w0 a& w0 S