Производное длинное уравнение в Matlab, затем вычисление результата

Question

У меня есть пример, подобный этому: у меня есть скорость одного элемента === производная это ====> ускорение элемента. Скорость рассчитывается по многим уравнениям, например:

    syms m w k E mI c w t L z;
E=2.7*10^10;
L=21;
mI=0.001332;
z=5;
m=500;
k=1*10^7;
c=30000;
w=300;
t=1;
a=((((m*w^2-k)/(E*mI))^2+(c*w/(E*mI))^2)^(1/8));
b=(atan(-c*w/(m*w^2-k))/4+pi/4);
syms cosi(w);
cosi=a*exp(t*b);
delta=2*(1-cosh(cosi*L)*cos(cosi*L));
syms v(w);
v=((sin(cosi*L)-sinh(cosi*L))*(cosh(cosi*z)+cos(cosi*z))+(cosh(cosi*L)...
    -cos(cosi*L))*(sinh(cosi*z)+sin(cosi*z)))...
    /(E*mI*cosi^3*delta);

f=jacobian(v,w);

С m,k,E,mI,c постоянны, w является независимой переменной. 1. Теперь я хочу вычислить производную от v 2. Теперь я хочу скомпилировать этот код в сборку Net, чтобы вызвать его с помощью программы на C#. Как я могу это исправить.

function [resultv,resultf]=Untitled4(m1, w1, k1, E1, mI1, c1, t1, L1, z1)
syms m w k E mI c w t L z cosi(w)
a=((((m*w^2-k)/(E*mI))^2+(c*w/(E*mI))^2)^(1/8));
b=(atan(-c*w/(m*w^2-k))/4+pi/4);
cosi=a*exp(t*b);
delta=2*(1-cosh(cosi*L)*cos(cosi*L));
v=((sin(cosi*L)-sinh(cosi*L))*(cosh(cosi*z)+cos(cosi*z))+(cosh(cosi*L)...
    -cos(cosi*L))*(sinh(cosi*z)+sin(cosi*z)))...
    /(E*mI*cosi^3*delta);

f=jacobian(v,w);

E=E1;
L=L1;
mI=mI1;
z=z1;
m=m1;
k=k1;
c=c1;
w=w1;
t=t1;
resultv=double(subs(v));
resultf = double(subs(f));
end

Большое спасибо.

Answer 1

Использование jacobian, Допустим, вы хотите получить WRT w

syms m w k E mI c w t L z 

a=((((m*w^2-k)/(E*mI))^2+(c*w/(E*mI))^2)^(1/8));
b=(atan(-c*w/(m*w^2-k))/4+pi/4);
cosi=(a*exp(t*b));
delta=2*(1-cosh(cosi*L)*cos(cosi*L));
v=((sin(cosi*L)-sinh(cosi*L))*(cosh(cosi*z)+cos(cosi*z))+(cosh(cosi*L)...
    -cos(cosi*L))*(sinh(cosi*z)+sin(cosi*z)))...
    /(E*mI*cosi^3*delta);

jacobian(v,w)

ОБНОВЛЕНИЕ ПОСЛЕ ВОПРОСА ОБНОВЛЕНО

syms m w k E mI c w t L z cosi(w)
a=((((m*w^2-k)/(E*mI))^2+(c*w/(E*mI))^2)^(1/8));
b=(atan(-c*w/(m*w^2-k))/4+pi/4);
cosi=a*exp(t*b);
delta=2*(1-cosh(cosi*L)*cos(cosi*L));
v=((sin(cosi*L)-sinh(cosi*L))*(cosh(cosi*z)+cos(cosi*z))+(cosh(cosi*L)...
    -cos(cosi*L))*(sinh(cosi*z)+sin(cosi*z)))...
    /(E*mI*cosi^3*delta);

f=jacobian(v,w);

E=2.7*10^10;
L=21;
mI=0.001332;
z=5;
m=500;
k=1*10^7;
c=30000;
w=300;
t=1;

results = double(subs(f))