2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Дифференциальное уравнение бесселя
Сообщение06.04.2010, 22:13 
Добрый день, уважаемые коллеги,

передо мной стоит задача решения дифференциального уравнения (бесселя) с последующим интегрированием полученных результатов. Задачка решалась в матлабе. Использовал для решения диф. ур. dsolve - вроде без проблем получилось, но при интегрировании полученные результаты не становятся числовыми, как я понимаю, решение остаётся символьным. В матлабе новичок, использую его как калькулятор.
Текст программы такой:

clear all
close all
clc
format
syms E0 unreal
syms r a b c unreal

f=0.01;
Edef=450.4;
fef=f/0.34;
hdef=0.9;
hd=0.34*hdef;

pd=1;



hm=(hd/f)-hd;
Em = 3.8;
num=0.4;
Gm = Em/(2+2*num);

pm=1;
p=1;
y1=1;
x1=1;


Ed=(Em*(1-fef)-Em*(1-f)+Edef*fef)/f;
nud=0.45;
Gd = Ed/(2+2*nud);

ld=150 ;


[R1m, R1d]=dsolve('(r^2)*D2R1d+r*DR1d-(r^2)*((12*Ed*R1d/(Gd*(ld^2)))+12*pd/(Gd*(ld^2)))' , '(r^2)*D2R1m+r*DR1m-(r^2)*((12*Em*R1m/(Gm*(ld^2)))+12*pm/(Gm*(ld^2)))' , 'R1d(hd)=R1m(hd)' , 'Gm*DR1m(hd)=Gd*DR1d(hd)')



Ad=sym('R1d*r');
Am=sym('R1m*r');

U1=pi*ld*p*(int(Ad, r, a, b)+int(Am, r, b, c));

U=subs(U1, {a,b,c}, {(0), (hd), (hd+hm)})

U0=(1/2)*ld*(hd+hm)*(p^2)/E0;
Eq5=U-U0;
Eeqv1=solve(Eq5,E0)


Вопрос в следующем:как в решении вместо символа получить число, в чём у меня ошибка?

 
 
 [ 1 сообщение ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group