Can you store symbolic functions in a matrix/vector?

21 次查看（过去 30 天）
KieranSQ
KieranSQ 2019-2-4
评论： Walter Roberson 2020-9-6
I am trying to write a finite element code. I want to consider the basis functions .
So I need to be able to store them in a matrix/vector such as: . In matlab I tried the following ,
phi=@(x,y)1/h.*[1-x-y;x;y]
but I either get an error for variables not being defined (which is the point at this stage). From here I want to do something like the following, but I can't see how. Any help would be appreciated!
phi(1)=1/h[1-x-y]

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2019-2-4
Symbolic functions cannot be stored in a regular vector or array. If you concatenate symbolic functions, then what you get back is not a vector of symbolic functions but rather a single symbolic function that returns a vector. This is because the () indexing operation for symbolic functions is always interpreted as invocation of the function, never as indexing.
This is also part of the same reason that function handles cannot be stored in a regular vector or array, because () for function handles always means invocation, never indexing.
You can store symbolic functions as elements of a cell array.
Symbolic expressions can be stored in a vector or array. Most of the time it is easier to work with symbolic expressions rather than symbolic functions. However if you need to express differentiation of a symbolic function whose definition is not fully known, such as a differential equation, then you need a symbolic function.

更多回答（1 个）

Sebastian Pakula
Sebastian Pakula 2020-9-6
There is a way of picking one of the function from matrix (or vector) of functions. Try this:
syms x y
h=10;
phi=(1/h)*[1-x-y;x;y];
% if you need second function write:
f2=[0 1 0]*phi;
% if you need third one:
f3=[0 0 1]*phi;
  1 个评论
Walter Roberson
Walter Roberson 2020-9-6
You do not show any functions, only expressions.
For expressions you can just index directly, phi(2) for example.
You also need to be careful about 0 times infinity, which would normally be nan.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Code Performance 的更多信息

标签

产品


版本

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by