Data must be numeric, datetime, duration or an array convertible to double.

1 回表示 (過去 30 日間)
I wrote a following code. It shows error:"Data must be numeric, datetime, duration or an array convertible to double."
clear all
clc
syms x
h=4;
q=0;
delta=26;
phi=39;
kh=0;
kv=0;
psi=atan(kh/(1-kv));
da1=delta*(pi/180); pha1=phi*(pi/180);
m=pha1+da1;
b=pha1-psi;
c=psi+da1;
gma=18.4; nq=2*q/(gma*(x));
delta=26;
phi=39;
lam=0;
A=lam*nq/(1+nq);
alphac=atan((sin(m)*sin(b)+(sin(m)^2+sin(b)^2+sin(m)*cos(m)*sin(b)*cos(b)+A*cos(c)*cos(m)*sin(b))^0.5)/(A*cos(c)+sin(m)*cos(b)));
w=0.5*gma*(x)^2*(1/tan(alphac));
B=(x)*((1/tan(alphac))-lam);
ka1=((1+nq)*(1-A*tan(alphac))*(cos(b)-(sin(b)/tan(alphac))))/(cos(psi)*(cos(m)+tan(alphac)*sin(m)))
pa1=0.5*gma*x^2*ka1
x=0:0.1:5;
plot(pa1,x)

採用された回答

madhan ravi
madhan ravi 2019 年 3 月 15 日
syms x
h=4;
q=0;
delta=26;
phi=39;
kh = 0;
kv=0;
s=...; % didn't define s
psi = atan (kh / (1-kv));
da1 = delta * (pi / 180); pha1 = phi * (pi / 180);
m=pha1+da1;
b=pha1-psi;
c=psi+da1;
gma = 18.4;
nq = 2 * q / (gma * (x));
delta=26;
phi=39;
lam = 0;
leg=...; % didn't define leg
Null=...;
A = lam * s / (1 + Null);
% ^---only one
alphac = atan ((sin (m) * sin (b) +...
(sin (m) ^ 2 + sin (b) ^ 2 +...
sin (m) * cos (m) * sin (b) * cos (b) +...
A * cos (c) * cos (m) * sin (b) ^ 0.5) / (A * cos (c) +...
sin (m) * cos (b))));
w=0.5*gma*(x)^2*(1/tan(alphac));
B = (x) * ((1 / tan (alphac)) - lam);
ka1 = ((1 + nq) * (1-A * tan (alphac)) *...
(cos (b) - (sin (b) / tan (alphac))))...
* (cos (psi) * (leg*(m ) + tan (alphac) * sin (m)))
% ^---- missed it
pa1 = matlabFunction(0.5 * gma * x ^ 2 * ka1);
x=0:0.1:5;
plot(pa1(x),x)

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeCalculus についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by