matlab用solve解方程组，解出来有很多组解，如何编程只保留全部未知数都大于0的一组解（该方

一：V1 = double(V1); V2 = double(V2); V3 = double(V3); Vges = double(Vges);
二：ind = find(V1>0&V2>0&V3>0&Vges>0);
三：V1_end = V1(ind); V2_end = V2(ind); V3_end = V3(ind); Vges_end = Vges(ind);

看看fsolve的源代码：>>typefsolvefunction[x,FVAL,EXITFLAG,OUTPUT,JACOB]=fsolve(FUN,x,options,varargin)%FSOLVEsolvessystemsofnonlinearequationsofseveralvariables.%%FSOLVEattemptstosolveequationsoftheform:%%F(X)=0whereFandXmaybevectorsormatrices.%%X=FSOLVE(FUN,X0)startsatthematrixX0andtriestosolvethe%equationsinFUN.FUNacceptsinputXandreturnsavector(matrix)of%equationvaluesFevaluatedatX.%%X=FSOLVE(FUN,X0,OPTIONS)solvestheequationswiththedefaultoptimization%parametersreplacedbyvaluesinthestructureOPTIONS,anargument%createdwiththeOPTIMSETfunction.SeeOPTIMSETfordetails.Used%optionsareDisplay,TolX,TolFun,DerivativeCheck,Diagnostics,%FunValCheck,Jacobian,JacobMult,JacobPattern,LineSearchType,%NonlEqnAlgorithm,MaxFunEvals,MaxIter,PlotFcns,OutputFcn,%DiffMinChangeandDiffMaxChange,LargeScale,MaxPCGIter,%PrecondBandWidth,TolPCG,andTypicalX.UsetheJacobianoptionto%specifythatFUNalsoreturnsasecondoutputargumentJthatisthe%JacobianmatrixatthepointX.IfFUNreturnsavectorFofm%componentswhenXhaslengthn,thenJisanm-by-nmatrixwhereJ(i,j)%isthepartialderivativeofF(i)withrespecttox(j).(Notethatthe%JacobianJisthetransposeofthegradientofF.)%%X=FSOLVE(PROBLEM)solvessystemdefinedinPROBLEM.PROBLEMisa%structurewiththefunctionFUNinPROBLEM.objective,thestartpoint%inPROBLEM.x0,theoptionsstructureinPROBLEM.options,andsolver%name'fsolve'inPROBLEM.solver.Usethissyntaxtosolveatthe%commandlineaproblemexportedfromOPTIMTOOL.ThestructurePROBLEM%musthaveallthefields.%%

我用其他软件的，应该会有一个选择条件的，可以筛选结果，
然后显示出来