色哟哟视频在线观看-色哟哟视频在线-色哟哟欧美15最新在线-色哟哟免费在线观看-国产l精品国产亚洲区在线观看-国产l精品国产亚洲区久久

0
  • 聊天消息
  • 系統消息
  • 評論與回復
登錄后你可以
  • 下載海量資料
  • 學習在線課程
  • 觀看技術視頻
  • 寫文章/發帖/加入社區
會員中心
創作中心

完善資料讓更多小伙伴認識你,還能領取20積分哦,立即完善>

3天內不再提示

關于MATLAB連續模型求解方法介紹和分析

MATLAB ? 來源:djl ? 作者:卓金武 ? 2019-09-16 11:41 ? 次閱讀

連續模型是指模型是連續函數的一類模型總稱,具體建模方法主要是微分方程建模。微分方程建模是數學建模的重要方法,因為許多實際問題的數學描述將導致求解微分方程的定解問題。把形形色色的實際問題化成微分方程的定解問題,大體上可以按以下幾步:

1.根據實際要求確定要研究的量(自變量、未知函數、必要的參數等)并確定坐標系。

2.找出這些量所滿足的基本規律(物理的、幾何的、化學的或生物學的等等)。

3.運用這些規律列出方程和定解條件。

MATLAB 在微分模型建模過程中的主要作用是求解微分方程的解析解, 將微分方程轉化為一般的函數形式。 另外, 微分方程建模, 一定要做數值模擬, 即根據方程的表達形式, 給出變量間關系的圖形, 做數值模擬也需要用 MATLAB 來實現。

微分方程的形式多樣,微分方程的求解也是根據不同的形式采用不同的方法, 在建模比賽中, 常用的方法有三種:

1.用 dsolve 求解常見的微分方程解析解

2.用 ODE 家族的求解器求解數值解

3.使用專用的求解器求解

1.常規微分方程的求解

微分方程在 MATLAB 中固定的表達方式, 這些基本的表達方式如下表所示:

關于MATLAB連續模型求解方法介紹和分析

對于通常的微分方程, 一般需要先求解析解, 那么 dsolve 是首先考慮的求解器,因為dsolve 能夠求解解析解,其具體的用法如下:

[實例]求微分方程 xy'+y-e^x=0 在初始條件 y(1)=2e 下的特解,并畫出解函數的圖形.

求解本問題的 Matlab 程序為:

syms x y

y=dsolve('x*Dy+y-exp(x)=0','y(1)=2*exp(1)','x')

ezplot(y)

微分方程的特解為:y=1/x*exp(x)+1/x*exp (1) (Matlab格式),即 y=(e+e^x)/x ,此函數的圖形如圖 1:

關于MATLAB連續模型求解方法介紹和分析

圖1y關于x的函數圖象

2.ODE 家族求解器

如果微分方程的解析形式求解不出來, 那么退而求其次的辦法是求解數值解,那么這個時候就需要用 ODE 家族的求解器求解微分方程的數值解啦。

因為沒有一種算法可以有效地解決所有的 ODE 問題,為此,MATLAB 提供了多種求解器,對于不同的 ODE 問題,采用不同的 Solver。MATLAB 中常用的微分方程數值解的求解器及特點如下表所示。

關于MATLAB連續模型求解方法介紹和分析

要特別提醒的是:ode23、ode45 是極其常用的用來求解非剛性標準形式一階常微分方程(組)初值問題解的 Matlab 的常用程序,其中:

ode23 采用龍格-庫塔2 階算法,用3 階公式作誤差估計來調節步長,具有低等的精度.

ode45 則采用龍格-庫塔4 階算法,用5 階公式作誤差估計來調節步長,具有中等的精度.

[實例]導彈追蹤問題

設位于坐標原點的甲艦向位于 x 軸上點 A(1, 0) 處的乙艦發射導彈,導彈頭始終對準乙艦。如果乙艦以最大的速度v0(是常數)沿平行于 y 軸的直線行駛,導彈的速度是 5*v0,求導彈運行的曲線方程,以及乙艦行駛多遠時,導彈將它擊中?

記導彈的速度為 w,乙艦的速率恒為v0。設時刻 t乙艦的坐標為 (X(t),Y(t)),導彈的坐標為 (x(t),y(t))。當零時刻,(X(0),Y(0))=(1,0),(x(0),y(0))=(0,0),建立微分方程模型:

關于MATLAB連續模型求解方法介紹和分析

因乙艦以速度v0沿直線 x=1 運動,設v0=1,w=5,X=1,Y=t,因此導彈運動軌跡的參數方程為:

關于MATLAB連續模型求解方法介紹和分析

MATLAB 求解數值解程序如下:

(1)定義方程的函數形式:

function dy=eq2(t,y)

dy=zeros(2,1);

dy(1)=5*(1-y(1))/sqrt((1-y(1))^2+(t-y(2))^2);

dy(2)=5*(t-y(2))/sqrt((1-y(1))^2+(t-y(2))^2);

(2)求解微分方程的數值解

t0=0,tf=0.21;

[t,y]=ode45('eq2',[t0 tf],[0 0]);

X=1;Y=00.21;plot(X,Y,'-')

plot(y(:,1),y(:,2),'*'),hold on

x=01; y=-5*(1-x).^(4/5)/8+5*(1-x).^(6/5)/12+5/24;

plot(x,y,'r')

關于MATLAB連續模型求解方法介紹和分析

3.專用求解器

對于復雜的微分方程模型的求解, 可以借助 MATLAB 偏微分方程工具箱中的專用求解器。以下將以一個實例來看看如何借助偏微分方程工具箱來實現一個微分方程的求解與數值仿真

所研究的對象是一個二階波的方程:

這個時候要查看一下 MALTAB 中哪個函數能求解相類似的方程, solvepde 可以求解的方程形式為:

可以發現只要通過參數設定就可以將所要求解的方程轉化成這種標準形式。

具體求解步驟如下:

(1) 設置參數

c = 1;

a = 0;

f = 0;

m = 1;

(2) 定義波的空間位置

numberOfPDE = 1;

model = createpde(numberOfPDE);

geometryFromEdges(model,@squareg);

pdegplot(model,'EdgeLabels','on');

ylim([-1.1 1.1]);

axis equal

title'Geometry With Edge Labels Displayed';

xlabel x

ylabel y

關于MATLAB連續模型求解方法介紹和分析

(3) 定義微分方程模型的系數和邊界條件

specifyCoefficients(model,'m',m,'d',0,'c',c,'a',a,'f',f);

applyBoundaryCondition(model,'dirichlet','Edge',[2,4],'u',0);

applyBoundaryCondition(model,'neumann','Edge',([1 3]),'g',0);

(4)定義該問題的有限元網格

generateMesh(model);

figure

pdemesh(model);

ylim([-1.1 1.1]);

axis equal

xlabel x

ylabel y

關于MATLAB連續模型求解方法介紹和分析

(5)定義初始條件

u0 = @(location) atan(cos(pi/2*location.x));

ut0 = @(location) 3*sin(pi*location.x).*exp(sin(pi/2*location.y));

setInitialConditions(model,u0,ut0);

(6)方程的求解

n = 31;% 求解次數

tlist = linspace(0,5,n);

model.SolverOptions.ReportStatistics ='on';

result = solvepde(model,tlist);

u = result.NodalSolution;

(7)模型的數值仿真

figure

umax = max(max(u));

umin = min(min(u));

for i = 1:n

pdeplot(model,'XYData',u(:,i),'ZData',u(:,i),'ZStyle','continuous',...

'Mesh','off','XYGrid','on','ColorBar','off');

axis([-1 1 -1 1 umin umax]);

caxis([umin umax]);

xlabel x

ylabel y

zlabel u

M(i) = getframe;

end

關于MATLAB連續模型求解方法介紹和分析

聲明:本文內容及配圖由入駐作者撰寫或者入駐合作網站授權轉載。文章觀點僅代表作者本人,不代表電子發燒友網立場。文章及其配圖僅供工程師學習之用,如有內容侵權或者其他違規問題,請聯系本站處理。 舉報投訴
  • 函數
    +關注

    關注

    3

    文章

    4345

    瀏覽量

    62867
  • 生物學
    +關注

    關注

    0

    文章

    26

    瀏覽量

    12991
  • 數值模擬
    +關注

    關注

    0

    文章

    29

    瀏覽量

    9605
收藏 人收藏

    評論

    相關推薦

    Simulink與 MATLAB 的結合使用 Simulink中的信號處理方法

    在工程和科學研究中,信號處理是一個重要的領域,涉及到信號的采集、分析、處理和生成。MATLAB 提供了豐富的信號處理工具箱,而 Simulink 提供了一個直觀的圖形界面,使得復雜的信號處理系統可以
    的頭像 發表于 12-12 09:25 ?451次閱讀

    介紹FIR濾波模型的建立,分4個步驟

    本帖介紹FIR濾波模型的建立,分以下幾個步驟: 選定濾波結構:低通、高通、帶通、帶阻; 選定合適的窗函數,常見的有hamming、hanning、blackman、ExactBlackman
    發表于 09-04 09:08

    matlab纜索動力學?

    關于這個我想獲得細纜在水中的張力和構型的matlab求解,編程時候這個向量我怎么辦
    發表于 08-29 15:57

    支路電流法是以什么為求解對象

    的定義 支路電流法是一種基于基爾霍夫電流定律(KCL)的電路分析方法。它通過將電路分解為若干個支路,并以這些支路的電流為變量,建立一組線性方程組來求解電路中的電流分布。 1.2 支路電流法的適用范圍 支路電流法適用于各種線性和非
    的頭像 發表于 08-08 17:00 ?1255次閱讀

    對放大電路的分析方法介紹

    放大電路是電子技術中非常重要的組成部分,它廣泛應用于信號處理、通信、測量等領域。對于放大電路的分析,通常可以分為直流分析和交流分析兩種情況。下面我們將介紹這兩種
    的頭像 發表于 08-07 10:08 ?770次閱讀

    【《大語言模型應用指南》閱讀體驗】+ 俯瞰全書

    的大語言模型設計技術人員閱讀,主要包括大語言模型的優化方法、Agent系統調優以及模型的安全技術。 展望篇分析了大語言
    發表于 07-21 13:35

    matlab預測模型怎么用

    MATLAB預測模型是一種基于統計和數學方法的預測工具,廣泛應用于各種領域,如金融、氣象、生物醫學等。本文將介紹MATLAB預測
    的頭像 發表于 07-11 14:33 ?668次閱讀

    MATLAB預測模型哪個好

    MATLAB中,預測模型的選擇取決于數據類型、問題復雜度和預測目標。以下是一些常見的預測模型及其適用場景的介紹: 線性回歸(Linear Regression): 線性回歸是最基本的
    的頭像 發表于 07-11 14:31 ?457次閱讀

    如何使用MATLAB創建預測模型

    在這篇文章中,我們將討論如何使用MATLAB創建預測模型MATLAB(矩陣實驗室)是一種用于數值計算、可視化和編程的高級編程語言和交互式環境。它廣泛應用于工程、科學研究和金融領域。 1.
    的頭像 發表于 07-11 14:29 ?726次閱讀

    matlab預測模型有哪些

    環境,使其成為預測模型開發和實現的理想平臺。本文將詳細介紹MATLAB中常用的預測模型及其應用。 線性回歸模型 線性回歸是一種簡單的預測
    的頭像 發表于 07-11 14:27 ?817次閱讀

    基于MATLAB 的質量守恒空化模型(JFO 模型

    可有大佬會基于MATLAB 的質量守恒空化模型(JFO 模型
    發表于 07-05 23:32

    matlab神經網絡工具箱結果分析

    神經網絡是一種強大的機器學習技術,廣泛應用于各種領域,如圖像識別、語音識別、自然語言處理等。MATLAB提供了一個功能強大的神經網絡工具箱,可以幫助用戶快速構建和訓練神經網絡模型。本文將介紹
    的頭像 發表于 07-03 10:32 ?674次閱讀

    基于MATLAB的信號處理系統與分析

    在信息技術日新月異的今天,信號處理與分析技術在通信、控制、生物醫學、圖像處理等領域的應用日益廣泛。MATLAB作為一款功能強大的數學軟件,為信號處理與分析提供了強大的工具和支持。本文旨在探討
    的頭像 發表于 05-17 14:24 ?1192次閱讀

    求助,關于AD采集到的數據分析問題

    MATLAB中畫圖后能看到幾個周期的圖像 數據特征:在matlab中能看到圖像是由兩部分構成,一部分是基波及其n次諧波,即存在上升沿,過沖,另一部分是隨機噪聲 求各位大神科普數據分析方法
    發表于 05-09 07:40

    gis中常用的空間分析方法

    將詳細介紹GIS中常用的空間分析方法,包括空間插值、緩沖區分析、空間統計、領域分析、網絡分析和多
    的頭像 發表于 02-25 13:44 ?5984次閱讀
    主站蜘蛛池模板: 久久热在线视频精品1| 92精品国产成人观看免费| 中文字幕a有搜索网站| 岛国大片在线播放免费| 久久yy99re66| 色99蜜臀AV无码| 嘴巴舔着她的私处插| 国产蜜臀AV在线一区视频| 欧美大片免费| 一二三四在线视频社区8| 国产av在线播放| 蜜芽手机在线观看| 亚洲精品一线二线三线无人区| 办公室激情在线观看| 久久无码人妻AV精品一区| 特黄特色大片免费播放器试看| 91麻豆久久| 精品国产自在天天线2019| 手机在线观看你懂的| aaaaaaa一级毛片| 九九99热久久999精品| 牲高潮99爽久久久久777| 9988电影网| 久99久热只有精品国产99| 乌克兰10一12x video| bl(高h)文| 久久这里只有热精品18| 亚洲精品成人A8198A片漫画| 大乳牛奶女magnet| 女厕所边摸边吃奶边做爽视频| 亚洲视频国产| 国产美女裸身网站免费观看视频 | 亚洲成AV人电影在线观看| XXX老姥群交| 麻豆沈芯语| 伊人大香线蕉精品在线播放| 国产在线观看免费观看不卡| 偷拍精品视频一区二区三区| 啊灬啊别停灬用力啊老师| 免费观看久久| 中文文字幕文字幕亚洲色|