Cannot figure out what went wrong

2018 10 10
2 个回答
回答已采纳
更新时间：2018 10 10
5 次查看（30 天）

0 个投票

Hi all, I am actually using Octave, but it is similar to Matlab so I figure I should get it a shoot here. I have a problem while trying to do numerical solution. I got the following error while trying to run the script.
When I tried to remove the ";" I will get the following error.
I also noticed that there is nothing in the workspace. The weird part is that the code works just fine on my friend's Octave. I just copy and paste the code into my computer, but it does not run the script. Can anyone please help?

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

KSSV
KSSV 2018-10-10
Copy the code here...not image.....image snippet is not going to hlep any manner for checking.
Daniel Yip
Daniel Yip 2018-10-10
z = 25;
l = 2500;
L0 = (z^2+l^2)^.5;
ks = 1.35;
EA = 50000000;
F_ext = 200;
v = [0:0.1:70]';
f_ext = EA*(L0-((z-v).^2+l^2).^0.5)/L0.*(z-v)/(((z-vI(i))^2+l^2)^0.5)+ks*v;
plot(v,f_ext); xlabel('v (mm)'); ylabel('f_e_x_t (N)');
vI = 0;
fI_ext = 0;
n = 50;
DF = F_ext/n;
Ki = [];
Dvi = [];
for i = 1 : n
Li = ((z-vI(i))^2+l^2)^0.5;
K = EA/L0 * ( (Li-L0)/Li+((z-vI(i))^2/(Li^2)) + ks;
Dv = K^-1*DF;
vI(i+1) = vI(i)+Dv;
fI_ext(i+1) = fI_ext(i)+DF;
Ki = [Ki;K];
Dvi = [Dvi;Dv];
end
hold on; plot(vI,fI_ext,'o--')
madhan ravi
madhan ravi 2018-10-10
@Daniel check Steven's answer

请先登录，再进行评论。

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

 采纳的回答

Steven Lord
Steven Lord 2018-10-10

1 个投票

Count your parentheses.

2 个评论
显示 无 隐藏 无

madhan ravi
madhan ravi 2018-10-10
+1 @steven scrutinised precisely.
madhan ravi
madhan ravi 2018-10-10
z = 25;
l = 2500;
L0 = (z^2+l^2)^.5;
ks = 1.35;
EA = 50000000;
F_ext = 200;
v = [0:0.1:70]';
% f_ext = EA*(L0-((z-v).^2+l^2).^0.5)/L0.*(z-v)/(((z-vI(i))^2+l^2)^0.5)+ks*v; error was here also because you can't subtract scalar with vector
%
% plot(v,f_ext); xlabel('v (mm)'); ylabel('f_e_x_t (N)');
vI = 0;
fI_ext = 0;
n = 50;
DF = F_ext/n;
Ki = [];
Dvi = [];
for i = 1 : n
Li = ((z-vI(i))^2+l^2)^0.5;
K = EA/L0 * ( (Li-L0)/Li+((z-vI(i))^2/(Li^2))) + ks; %error here as steven said one closing parentheses was missing here
Dv = K^-1*DF;
vI(i+1) = vI(i)+Dv;
fI_ext(i+1) = fI_ext(i)+DF;
Ki = [Ki;K];
Dvi = [Dvi;Dv];
end
hold on; plot(vI,fI_ext,'o--')

请先登录，再进行评论。

更多回答（1 个）

类别

在 帮助中心 和 File Exchange 中查找有关 C Shared Library Integration 的更多信息

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by