How to code for an equation for numbers of input?

3 次查看（过去 30 天）
Khairun Nisa Khamil
评论： Khairun Nisa Khamil 2018-2-15
Hi,
I am new in Matlab.
Would like to know how to code for this problem.
I have this equation
Qc = 2*N*((alpha*I0*Tc)-(0.5*I0*I0*(rho/G))-(k*G*(Th-Tc)))
And i want to compute this equation for I0 from 1 to 10 and Tc from -60 to 240.
I try coding it this way.
//
N= 96;
G = 0.072;
alpha = 0.000203;
rho = 0.001095;
k = 0.01438;
Th = 323;
for I0 = 0.0
for Tc = -60.0
Qc = 2*N*((alpha*I0*Tc)-(0.5*I0*I0*(rho/G))-(k*G*(Th-Tc)))
I0 = I0+1
Tc = Qc+30
if I0>=11 && Tc >=270
end
end
end
But, it only give me 1 answer instead of I need the Qc to give me the output in terms of I(1:10) and Tc (-60:240)
Please help

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

采纳的回答

Walter Roberson
Walter Roberson 2018-2-13
N = 96;
G = 0.072;
alpha = 0.000203;
rho = 0.001095;
k = 0.01438;
Th = 323;
[I0, Tc] = ndgrid(1:10, -60:240);
Qc = 2*N*((alpha .* I0 .* Tc) - (0.5 .* I0 .* I0 .* (rho./G)) - (k .* G .* (Th-Tc)));
surf(I0, Tc, Qc, 'edgecolor', 'none')
  4 个评论
显示 2更早的评论
Walter Roberson
Walter Roberson 2018-2-14
Yes, it is possible. See https://blogs.mathworks.com/graphics/2014/10/21/double_pendulum/ for a more complicated example. I do not think you need to create a class for this, but you should have the general structure of an initialization phase and an update phase.
Khairun Nisa Khamil
Khairun Nisa Khamil 2018-2-15
The example is very complicated for me to understand.
i tried doing it this way using matlab function in simulink.
I used this code
function [y1,y2,y3] = fcn(u)
h = u(1);
deltatmax = u(2);
Imax = u(3);
Umax = u(4);
z0 = (2.*u(2))/((u(1)-u(2))^2);
alphaM = u(4)/u(1);
rhoM = ((u(1)-u(2))*u(4)*u(3))/(2*u(1)*u(2));
kM = ((u(1)-u(2))*u(4))/(u(1)*u(3));
N = 127;
[I0, Tc] = ndgrid(0:10, 300:-20:240);
y1 = ((alphaM .* I0 .* Tc) - (0.5 .* I0 .* I0 .*rhoM) - (kM .*(u(1)-Tc)));
y2 = alphaM.*(u(1)-Tc)+(I0.*rhoM);
y3 = y1./(y2.*I0);
figure (1);
plot (I0,y1);
figure (2);
plot (y2,y1);
But it doesn't give the output i want in workspace.
The results i want is the graph on left but the scope shows result on right. Not sure what i did wrong. Or should i use S Function instead?

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Creating, Deleting, and Querying Graphics Objects 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by