Trying to Plot a Function

2022 1 31
1 个回答
更新时间：2022 1 31
1 次查看（30 天）

0 个投票

Working on a Nuclear Engineering problem, we have to plot the relationship between Keff and the diameter of a cylinder. Here is where I'm stuck:
>> %%Keff%%
nu=2.65;
Density=19.1;
MW=235;
N=(Density/MW)*6.022*10^23;
MicF=1.235*10^-24;
MicA=1.335*10^-24;
MicS=4.566*10^-24;
MacF=MicF*N;
MacA=MicA*N;
MacS=MicS*N;
A=235;
MuBar=2/(3*A);
D=1/(3*MacS*(1-MuBar));
H=30;
Diam=(1:180);
Keff=(nu*MacF)/(MacA+(D*((pi/H)^2+((2*2.4050)/Diam)^2)));
plot(Keff,Diam);
xlabel("Keff");
ylabel("Diameter");
title("Relationship Between Keff and Diameter");
Error using /
Matrix dimensions must agree.
Any ideas on how to fix this code. I just want a plot of the relationship.

5 个评论
显示 3更早的评论 隐藏 3更早的评论

Arif Hoq
Arif Hoq 2022-1-31
编辑：Arif Hoq 2022-1-31
Ran in:
Ran in:
Try with this
clear;
clc;
nu=2.65;
Density=19.1;
MW=235;
N=(Density/MW)*6.022*10^23;
MicF=1.235*10^-24;
MicA=1.335*10^-24;
MicS=4.566*10^-24;
MacF=MicF*N;
MacA=MicA*N;
MacS=MicS*N;
A=235;
MuBar=2/(3*A);
D=1/(3*MacS*(1-MuBar));
H=30;
Diam=(1:180);
Keff=(nu*MacF)./(MacA+(D*((pi/H)^2+((2*2.4050)./Diam).^2))); % element by element
plot(Keff,Diam);
xlabel("Keff");
ylabel("Diameter");
title("Relationship Between Keff and Diameter");
I think you have to create your matrix element by element, like ./ and .^
For example
A = [1 1 0 0];
B = [1; 2; 3; 4];
C = A*B
C = 3
D= A.*B % element by element
D = 4×4
1 1 0 0 2 2 0 0 3 3 0 0 4 4 0 0
Kevin Smith
Kevin Smith 2022-1-31
Arif: you are a steely-eyed missle man! Thank you!
Arif Hoq
Arif Hoq 2022-1-31
@Kevin Smith Thank you very much for your compliment. If you find your answer, then please consider marking it as accepted answer. That will make it easier for others with similar issues to find an answer.
KSSV
KSSV 2022-1-31
@Arif Hoq you have commented the code. It cannot be accepted unless you copy it in answer section.
Arif Hoq
Arif Hoq 2022-1-31
@KSSV thanks a lot. I was not such aware of it.

请先登录，再进行评论。

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

回答（1 个）

Arif Hoq
Arif Hoq 2022-1-31
Ran in:

0 个投票

Ran in:
Try with this:
clear;
clc;
nu=2.65;
Density=19.1;
MW=235;
N=(Density/MW)*6.022*10^23;
MicF=1.235*10^-24;
MicA=1.335*10^-24;
MicS=4.566*10^-24;
MacF=MicF*N;
MacA=MicA*N;
MacS=MicS*N;
A=235;
MuBar=2/(3*A);
D=1/(3*MacS*(1-MuBar));
H=30;
Diam=(1:180);
Keff=(nu*MacF)./(MacA+(D*((pi/H)^2+((2*2.4050)./Diam).^2))); % element by element
plot(Keff,Diam);
xlabel("Keff");
ylabel("Diameter");
title("Relationship Between Keff and Diameter");
I think you have to create your matrix element by element, like ./ and .^
For example
A = [1 1 0 0];
B = [1; 2; 3; 4];
C = A*B
C = 3
D= A.*B % element by element
D = 4×4
1 1 0 0 2 2 0 0 3 3 0 0 4 4 0 0

类别

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

标签

编辑：

KSSV
KSSV
2022-1-31

Community Treasure Hunt

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

Start Hunting!

Translated by