How can I substitute a variable given by a struct?

5 ビュー (過去 30 日間)
Abner Ojeda
Abner Ojeda 2020 年 5 月 2 日
コメント済み: Abner Ojeda 2020 年 5 月 3 日
Hey, everyone. I'm trying to solve an optimization problem without the Optimization toolbox. I already have the roots from the equation (by "solve"), it gives to me an struct array, with fields x1, x2... xn, my problem is I already have an x1, x2... xn. I've been trying by so many ways but simply I couldn't solve it, I'm looking to asign the value of the fields to the syms variables and then apply "subs". Code below.
Thanks in advance for help. :)
clear; close all; clc;
var= input ('Dimention of domain');
syms x [1 var];
func= input ('Enter the function, i.e [x1+x2+...+xn]: \n');
y=[];
for s=x(1:end)
dp= diff (func, [s]);
y=[y, dp];
end
pc=solve (y);
pucr= fieldnames (pc);
for i= 1:numel(pucr)
vals(i, 1)= getfield(pc,pucr{i});
end
Hf=hessian (func);
Hff= subs (Hf);
  2 件のコメント
darova
darova 2020 年 5 月 2 日
Try this
for i= 1:numel(pucr)
v = getfield(pc,pucr{i});
v1 = subs(v,oldv, newv);
vals(i, 1) = doulbe(v1);
end
Abner Ojeda
Abner Ojeda 2020 年 5 月 3 日
Thanks for taking your time in answer my question, it doesn't work, but I appreciate your answer.

サインインしてコメントする。

採用された回答

Walter Roberson
Walter Roberson 2020 年 5 月 2 日
subs(ExpressionToSubstituteInto, pc)
When pc is a struct with field names, then the corresponding variable names will have the appropriate values substituted.
Note: this will not work if any of the field names contain an empty value, even if the corresponding variable is not used in the expression. In particular, if you happened to do
pc = solve(y, 'returnconditions', true)
then pc.parameters would be created even if no parameters were generated, and with pc.parameters being empty, subs() with pc would fail. You can get it working by removing the empty field from the struct.
  1 件のコメント
Abner Ojeda
Abner Ojeda 2020 年 5 月 3 日
It works, Walter, it works pretty well. Thank you.

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeMathematics についてさらに検索

製品


リリース

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by