Issue in plotting of output after function

2021 9 8
1 个回答
回答已采纳
更新时间：2021 9 8
4 次查看（30 天）

0 个投票

Unable to plot output after function.
I gave my script below
syms theta
A = cos(theta).*[1 5 6; 2 9 3; 5 1 0];
B = tan(2*theta).*[8 5 1; 0 1 6; 2 4 7];
myfun = @(t,y)scriptname(t,y,A,B);
% dummy values for tspan and y0
tspan = [0 1];
y0 = [0; 0; 0];
% ode solver
sol = ode45(myfun,tspan,y0);
function dydt = scriptname(t,y,A,B)
V = [1; 5; 3];
% evaluation of A and B (numerical) with theta = y(3)
An = double(subs(A,y(3)));
Bn = double(subs(B,y(3)));
dydt = V-((An+Bn)*y);
end
for i=1.:3.
stator=y(:,i);
figure(1.);
subplot(2.,2.,i);
plot(t,stator);
xlabel('time');
ylabel('Current');
title(['stat current # ',int2str(i)]);
end
iam getting error like this
Function definitions in a script must appear at the end of the file.
Move all statements after the "scriptname" function definition to before the first local function definition.
How to resolve this??
Thanks in advance

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

 采纳的回答

KSSV
KSSV 2021-9-8
编辑：KSSV 2021-9-8

0 个投票

function myfun()
syms theta
A = cos(theta).*[1 5 6; 2 9 3; 5 1 0];
B = tan(2*theta).*[8 5 1; 0 1 6; 2 4 7];
myfun = @(t,y)scriptname(t,y,A,B);
% dummy values for tspan and y0
tspan = [0 1];
y0 = [0; 0; 0];
% ode solver
sol = ode45(myfun,tspan,y0);
%% Plot here
plot(sol.x,sol.y)
function dydt = scriptname(t,y,A,B)
V = [1; 5; 3];
% evaluation of A and B (numerical) with theta = y(3)
An = double(subs(A,y(3)));
Bn = double(subs(B,y(3)));
dydt = V-((An+Bn)*y);
Save the aboove code in myfun.m and run it.

10 个评论
显示 8更早的评论 隐藏 8更早的评论

KSSV
KSSV 2021-9-8
Check Sol, it will be a structure with required values. See it's fields and plot what you want.
Bathala Teja
Bathala Teja 2021-9-8
Iam new to matlab.
i didnt get what your are suggesting.
I want the output waveforms after compiling script.
Can you explian clearly.
KSSV
KSSV 2021-9-8
I have shown it in code already....check it.
Bathala Teja
Bathala Teja 2021-9-8
After compiling your code(which you gave), it is giving this error
Unable to resolve the name Sol.x.
Error in Untitled5 (line 151)
plot(Sol.x,sol.y)
KSSV
KSSV 2021-9-8
The variable is sol not sol (No capitlal S). Typo error. Edited the code.
Bathala Teja
Bathala Teja 2021-9-8
Is there any way to plot without making entire script as function(means without using function myfun() at starting)???
KSSV
KSSV 2021-9-8
syms theta
A = cos(theta).*[1 5 6; 2 9 3; 5 1 0];
B = tan(2*theta).*[8 5 1; 0 1 6; 2 4 7];
myfun = @(t,y)scriptname(t,y,A,B);
% dummy values for tspan and y0
tspan = [0 1];
y0 = [0; 0; 0];
% ode solver
sol = ode45(myfun,tspan,y0);
%% Plot here
plot(sol.x,sol.y)
Write the above in a min file. Save the below into a function.
function dydt = scriptname(t,y,A,B)
V = [1; 5; 3];
% evaluation of A and B (numerical) with theta = y(3)
An = double(subs(A,y(3)));
Bn = double(subs(B,y(3)));
dydt = V-((An+Bn)*y);
Now run the main file.
Bathala Teja
Bathala Teja 2021-9-8
sorry to bothering you again, one last doubt.
i followed your 1st method(keeping whole script as function), i am getting all waveforms in one graph. How to plot each solution separatley.
KSSV
KSSV 2021-9-8
figure
plot(sol.x,sol.y(1,:))
figure
plot(sol.x,sol.y(2,:))
figure
plot(sol.x,sol.y(3,:))

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Mathematics and Optimization 的更多信息

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by