|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
使用fsolve求解方程,显示no solution found .但此时仍然会赋于未知数一个值;请问如何丢弃这个结果(或者说让其不赋值),程序在一个大循环中,无法一个一个丢弃;例如:% [& ~ ~- F, x" I5 `0 s1 G% T9 e' _
clc;, P3 D/ Z P3 R6 O9 D% c
clear;! ]. y$ F9 { }- s) o+ P+ b2 t9 V
x1=[];
! `# }: ]+ X: \f1=@(x)([x(1)-x(2);x(1)+2-x(2)]);4 {& ]* o ^ t; H: b. ]
options = optimoptions('fminunc', 'TolX', 1e-100, 'MaxIter', 1000);
+ H$ O6 w% S4 S8 J; a! M5 n3 ^x=fsolve(f1,[1,1],options);! t& w* _0 s! G
if isempty(x)8 W" B/ q: P* G8 q5 k! T
throw(x);
* m3 J+ i/ u$ I: Xelse/ F3 v3 _2 k; m
x1=[x1,x];& P! S% Y" w4 Y. x' N$ p1 Y# N
end, s, W* l& E; j4 M2 F
得到
+ v$ ~/ M& G% f5 G, t>> x
! D0 _' G. H1 l' b( N# Fx =
4 }* W) D$ |; T5 J( X2 x( ^ 0.5000 1.5000
, x: s: M/ y" r& l/ K7 O: k>> x14 a _! i; e% Q: |- q+ O" y
x1 =6 y" o1 W9 C7 Y9 ^% c) R- {$ j) {6 p
0.5000 1.5000。
0 ]' w) ]+ m. }, l; s- g
7 \3 [9 |) W1 s. ^9 r5000;" @8 A/ M1 g" [+ ` e/ r
请教大神,非常感谢!. p3 \! k* P; x+ D1 e
|
|