Unrecognized function or variable in a code

Question

Mathew 2024-6-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2130181-unrecognized-function-or-variable-in-a-code

编辑： Mathew 2024-6-19

% param value
clear all; close all;
% p is constant and q varies
k0 = 0.0000169; u=0.6; p=536.2; q=0.0000376;
% your function
 k = @(t,u) (k0 + (p*k0.^u-q*k0).*t.^a/gamma(a+1) +(p*k0.^u-q*k0).*(p*u*k0.^u-q).*t.^(2*a)/gamma(2*a+1)...
  +(p*k0.^u-q*k0).*((p*u*k0.^u-q).^2-p*u.*(u-1)*k0.^(u-1)./2).*t.^(3*a)/gamma(3*a+1));
% grid
%t = linspace(0,0.1);
%u = linspace(0.3,0.9);
%[T,U] = meshgrid(t,q);
t = linspace(0,0.1,50);
a = linspace(0.3,0.9,50);
[T,A] = meshgrid(t,a);
% evaluate function
Z = k(T,A);
Unrecognized function or variable 'a'.

Error in solution>@(t,u)(k0+(p*k0.^u-q*k0).*t.^a/gamma(a+1)+(p*k0.^u-q*k0).*(p*u*k0.^u-q).*t.^(2*a)/gamma(2*a+1)+(p*k0.^u-q*k0).*((p*u*k0.^u-q).^2-p*u.*(u-1)*k0.^(u-1)./2).*t.^(3*a)/gamma(3*a+1)) (line 6)
 k = @(t,u) (k0 + (p*k0.^u-q*k0).*t.^a/gamma(a+1) +(p*k0.^u-q*k0).*(p*u*k0.^u-q).*t.^(2*a)/gamma(2*a+1)...
% plot
figure
surf(T,A,Z)
%surf(T,Q,Z,'facecolor','none')
xlabel('t');
ylabel('\mu');
zlabel('k(t)')

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Cris LaPierre 2024-6-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2130181-unrecognized-function-or-variable-in-a-code#answer_1474461

You create an anonymous function, k, that uses a variable a in the function. However, this variable was not defined at the time the anonymous function was defined.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Mathew 2024-6-19

在 MATLAB Online 中打开

Here is exactly the code that I wanted to run.

The code does not produce an image with MATLAB 2023b.

Warning: Matrix is singular to working precision.

% Warning: Matrix is singular to working precision.

% That is the response that I get.

% param value

clear all; close all;

% p is constant and q varies

k0 = 0.0000169; u=0.6; p=536.2; q=0.0000376;

k = @(t,a) (k0 + (p*k0.^u-q*k0).*t.^a/gamma(a+1) +(p*k0.^u-q*k0).*(p*u*k0.^u-q).*t.^(2*a)/gamma(2*a+1)...

+(p*k0.^u-q*k0).*((p*u*k0.^u-q).^2-p*u.*(u-1)*k0.^(u-1)./2).*t.^(3*a)/gamma(3*a+1));

% grid

t = linspace(0,1,50);

a = linspace(0.25,0.85,50);

[T,A] = meshgrid(t,a);

% evaluate function

Z = k(T,A);

Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 4.818625e-53.

Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.825871e-53.

Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 6.413452e-54.

% plot

figure

surf(T,A,Z)

xlabel('t');

ylabel('\alpha');

zlabel('k(t)')

Steven Lord 2024-6-19

在 MATLAB Online 中打开

You likely want to use array division (element-wise) rather than matrix division.

% p is constant and q varies

k0 = 0.0000169; u=0.6; p=536.2; q=0.0000376;

k = @(t,a) (k0 + (p*k0.^u-q*k0).*t.^a./gamma(a+1) +(p*k0.^u-q*k0).*(p*u*k0.^u-q).*t.^(2*a)./gamma(2*a+1)...

+(p*k0.^u-q*k0).*((p*u*k0.^u-q).^2-p*u.*(u-1)*k0.^(u-1)./2).*t.^(3*a)./gamma(3*a+1));

% grid

t = linspace(0,1,50);

a = linspace(0.25,0.85,50);

[T,A] = meshgrid(t,a);

% evaluate function

Z = k(T,A);

% plot

figure

surf(T,A,Z)

xlabel('t');

ylabel('\alpha');

zlabel('k(t)')

Mathew 2024-6-19

编辑：Mathew 2024-6-19

Thanks @ Steven Lord and @ Cris LaPierre

请先登录，再进行评论。

Unrecognized function or variable in a code

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Unrecognized function or variable in a code

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

4 个评论 显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

WeChat

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

4 个评论
显示 2更早的评论隐藏 2更早的评论