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

busahuang

function [TD] = menxian2(L)
# e( W/ K9 ^; E" u1 j%UNTITLED4 此处显示有关此函数的摘要
%   此处显示详细说明
8 C# r0 l1 {* D/ m' esyms 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);
TD=f1(1-L);
$ k% [6 J' L; d+ B$ C! V) Cend
. h( ^1 l# D" o# P
代码如上，解不了图片中要求得值，请教大神，非常感谢！
8 Q  F2 Y4 `4 h# R6 g6 ?9 I( D7 o
关于这个函数得解，Pfa已知，求Td

6 T. X3 d" o4 h# }: ?

4 h5 @/ C; j; ~5 e6 w
* Z: C, q; `/ h! J2 I" @

Ferrya

$ Q9 M& h; m% _% [! H4 k% I# E' X直接用TD = norminv( 1- L );就行了
非要解方程的话用vpasolve求数值解
0 p# a, L/ B, p5 g' [# Ysyms x t real
s = int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
  F) y$ W+ n1 N: q: Q+ Y- jy = matlabFunction(s);
TD = vpasolve( y(x)==1-L );

busahuang

Ferrya 发表于 2020-5-28 15:03
4 B: O0 T) |  c直接用TD = norminv( 1- L );就行了
0 g! t+ l8 \: |; q3 J非要解方程的话用vpasolve求数值解
syms x t real

. k2 M2 w. [, Q2 P1 F9 F大佬看看求这个Td用的代码跟前面这个一样吗，n已知。


+ L( z& t0 ]  O6 d/ @/ X9 |: G# Z