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

busahuang

function [TD] = menxian2(L)
%UNTITLED4 此处显示有关此函数的摘要
%   此处显示详细说明
syms s y x t
, T- [# y8 {5 b% Y6 ts=int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
: m7 P6 R0 [! A' Ly=matlabFunction(s);
7 P6 y" |; }8 g8 Z/ ~+ x; w  Nf1=finverse(s,t);
f2=matlabFunction(f1);
TD=f1(1-L);
( ?6 i% q# v( |* \( Zend
; t2 w7 O0 [! y) v! w( ~
代码如上，解不了图片中要求得值，请教大神，非常感谢！

关于这个函数得解，Pfa已知，求Td
2 _- }$ P: B0 o
: q; u* d9 }% g* X/ m2 `/ R$ L8 h
+ Q' q5 Z. n8 n' Y/ R

busahuang

Ferrya 发表于 2020-5-28 15:03
直接用TD = norminv( 1- L );就行了
5 M* _7 x: C8 E& e) [非要解方程的话用vpasolve求数值解
syms x t real

大佬看看求这个Td用的代码跟前面这个一样吗，n已知。
" @" K# j3 Q( E
! @, H6 u8 m" `

Ferrya

直接用TD = norminv( 1- L );就行了
非要解方程的话用vpasolve求数值解
syms x t real
s = int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
y = matlabFunction(s);
TD = vpasolve( y(x)==1-L );