Проверка на началното условие:
>> exp(0)
ans =
1
Втори начин на решаване - числено
1. Записваме уравнението като функция във файл oneurawn.m в m-редактора:
function damp=oneurawn(t,x)
damp=-3*x
2. Извикваме odе45:
>>[t,х]=ode45(@oneurawn,[0 4],1)
3. Резултатите се получават в workspace в променливата t – стойности на времето t и в променливата x - стойност на функцията x(t). Изчертаваме графиката на x(t).
>> plot(t,x)
Допълнително могат да се поставят етикети по осите х и y и заглавие на графиката с командите:
>>xlabel('t')
>> ylabel('x(t)')
>> title('Solution of the equation dx/dt=-3x(t)')
Графиката изглежда така:
14.В Command Window пишем:
>>diary off
за да се прекрати дневника на работата.
15.От менюто File избираме Save Workspace as и записваме Workspace във файл.
16.Маркираме всички редове в Command History. Натискаме десен бутон на мишката върху маркираната област. От контекстното меню избираме Create M-File и се отваря М-редактора. Оттам избираме File/Save as и записваме файла.
Всички файлове се съхраняват в директорията, която Matlab подразбира за записване на файловете си. Тя е работната директория на продукта.
Сподели с приятели: |