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

busahuang

function [TD] = menxian2(L)
%UNTITLED4 此处显示有关此函数的摘要
2 T  i5 t" s5 J8 [8 `: l%   此处显示详细说明
syms s y x t
s=int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
y=matlabFunction(s);
f1=finverse(s,t);
f2=matlabFunction(f1);
1 \) w) k, ^$ C: _) e; j/ QTD=f1(1-L);
. h2 Y! A# P+ R( L% M$ U) Rend
) w& ?5 g7 X' A: Z" S& A
代码如上，解不了图片中要求得值，请教大神，非常感谢！
% M. S* [0 h5 S5 f4 y; H! Y5 ]
关于这个函数得解，Pfa已知，求Td

Ferrya

直接用TD = norminv( 1- L );就行了
: s3 G/ X' U: }: B$ W非要解方程的话用vpasolve求数值解
& u; W' s! J4 c! T, ^5 A: Qsyms x t real
s = int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
y = matlabFunction(s);
1 Z, J2 e; ]/ k% K8 j! pTD = vpasolve( y(x)==1-L );

busahuang

Ferrya 发表于 2020-5-28 15:03
1 k7 [8 {6 O: w$ N直接用TD = norminv( 1- L );就行了
非要解方程的话用vpasolve求数值解
+ q: z# i7 x; X9 ?% z% Rsyms x t real

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

5 J8 j# F/ d( e1 U
& L( H0 X8 m. G! \/ B  I9 s4 M