 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 1859 | 回復(fù): 9 | ||
xhw2211金蟲 (正式寫手)
|
[求助]
mathematica 編程求助
|
|
頭一次用mathtematica編的程序，總是沒有畫出圖形，不知道錯(cuò)在什么地方了，大�？匆幌�。 z1 = LaplaceTransform[Exp[-10*t - 5*t^2], t, s]; z2 = LaplaceTransform[(10 + 10*t)*Exp[-10*t - 5*t^2], t, s]; z3 = Simplify[z2]; equ1 = z1 + z3*y - z; y1 = LaplaceTransform[7*Exp[-7*t], t, s]; y2 = LaplaceTransform[110*Exp[-110*t], t, s]; equ2 = y1*x + y2*z - y; x1 = LaplaceTransform[120*Exp[-120*t], t, s]; equ3 = x1*y - x; Solve[{equ1 == 0, equ2 == 0, equ3 == 0}, {x, y, z}] m = InverseLaplaceTransform[z, s, t] Plot[z, {t, 0, 1}, PlotStyle -> {Red, Thickness[0.008]}] |
回復(fù)此樓» 猜你喜歡
291求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有13人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有5人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有17人回復(fù)
-
考研調(diào)劑
已經(jīng)有7人回復(fù)
-
343求調(diào)劑085601
已經(jīng)有3人回復(fù)
-
08工科求調(diào)劑286
已經(jīng)有3人回復(fù)
-
337求調(diào)劑
已經(jīng)有3人回復(fù)
-
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求mathematica編程
已經(jīng)有6人回復(fù)
- 【求助】開平方問題【已完結(jié)】
已經(jīng)有3人回復(fù)
- 【求助】求解氫原子電子密度的時(shí)候遇到的4個(gè)問題
已經(jīng)有4人回復(fù)
- 【求助】如何利用matlab計(jì)算這個(gè)無窮級(jí)數(shù)的和
已經(jīng)有12人回復(fù)
- 【求助】100金幣求助高人一個(gè)程序解方程
已經(jīng)有15人回復(fù)
- 【求助】有關(guān)fortran進(jìn)行數(shù)值模擬
已經(jīng)有10人回復(fù)
jjdg
版主 (知名作家)
- 程序強(qiáng)帖: 9
- 應(yīng)助: 223 (大學(xué)生)
- 貴賓: 3.926
- 金幣: 51005.2
- 紅花: 238
- 帖子: 7227
- 在線: 474.2小時(shí)
- 蟲號(hào): 27000
- 注冊(cè): 2003-10-21
- 專業(yè): 藥物代謝與藥物動(dòng)力學(xué)
- 管轄: 程序語言
2樓2011-07-25 17:51:26
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時(shí)
- 蟲號(hào): 570836
- 注冊(cè): 2008-06-08
- 專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
lot[m, {t, 0, 1}, PlotStyle -> {Red, Thickness[0.008]}] 3樓2011-07-25 18:33:04
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時(shí)
- 蟲號(hào): 570836
- 注冊(cè): 2008-06-08
- 專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
4樓2011-07-25 18:35:01
mshwangg
至尊木蟲 (正式寫手)
- 程序強(qiáng)帖: 5
- 應(yīng)助: 206 (大學(xué)生)
- 金幣: 10702.8
- 散金: 100
- 紅花: 19
- 帖子: 597
- 在線: 195.4小時(shí)
- 蟲號(hào): 576702
- 注冊(cè): 2008-06-21
- 專業(yè): 物理學(xué)I
5樓2011-07-25 21:47:56
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時(shí)
- 蟲號(hào): 570836
- 注冊(cè): 2008-06-08
- 專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
|
我按你的方法改了程序，但是還是沒有結(jié)果，感覺是不是拉普拉斯變換后解出的Z表達(dá)式太復(fù)雜，mathematica無法進(jìn)行拉普拉斯反變換？？？下面是程序運(yùn)行后得到的m m={1/2 Sqrt[\[Pi]/5] InverseLaplaceTransform[( E^(5 + s + s^2/20) (1 - 840/((7 + s) (120 + s))) Erfc[(10 + s)/( 2 Sqrt[5])])/( 1 - 840/((7 + s) (120 + s)) - ( 110 (1 - 1/2 E^(1/20 (10 + s)^2) Sqrt[\[Pi]/5] s Erfc[(10 + s)/(2 Sqrt[5])]))/(110 + s)), s, t]} |
6樓2011-07-26 09:25:25
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時(shí)
- 蟲號(hào): 570836
- 注冊(cè): 2008-06-08
- 專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
7樓2011-07-26 09:26:19
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時(shí)
- 蟲號(hào): 570836
- 注冊(cè): 2008-06-08
- 專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
|
把修改后的程序貼出來： z1 = LaplaceTransform[Exp[-10*t - 5*t^2], t, s]; z2 = LaplaceTransform[(10 + 10*t)*Exp[-10*t - 5*t^2], t, s]; z3 = Simplify[z2]; equ1 = z1 + z3*y - z; y1 = LaplaceTransform[7*Exp[-7*t], t, s]; y2 = LaplaceTransform[110*Exp[-110*t], t, s]; equ2 = y1*x + y2*z - y; x1 = LaplaceTransform[120*Exp[-120*t], t, s]; equ3 = x1*y - x; rr = Solve[{equ1 == 0, equ2 == 0, equ3 == 0}, {x, y, z}]; z = z /. rr; m = InverseLaplaceTransform[z, s, t] Plot[z, {t, 0, 1}, PlotStyle -> {Red, Thickness[0.008]}] |
8樓2011-07-26 09:27:33
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時(shí)
- 蟲號(hào): 570836
- 注冊(cè): 2008-06-08
- 專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
★
jjdg(金幣+1): 感謝參與 2011-07-30 14:50:08
jjdg(金幣+1): 感謝參與 2011-07-30 14:50:08
|
下面這個(gè)是對(duì)的： z1 = LaplaceTransform[Exp[-10*t - 5*t^2], t, s]; z2 = LaplaceTransform[(10 + 10*t)*Exp[-10*t - 5*t^2], t, s]; z3 = Simplify[z2]; equ1 = z1 + z3*y - z; y1 = LaplaceTransform[7*Exp[-7*t], t, s]; y2 = LaplaceTransform[110*Exp[-110*t], t, s]; equ2 = y1*x + y2*z - y; x1 = LaplaceTransform[120*Exp[-120*t], t, s]; equ3 = x1*y - x; rr = Solve[{equ1 == 0, equ2 == 0, equ3 == 0}, {x, y, z}]; z = z /. rr; m = InverseLaplaceTransform[z, s, t] Plot[m, {t, 0, 1}, PlotStyle -> {Red, Thickness[0.008]}] |
9樓2011-07-26 09:33:09
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時(shí)
- 蟲號(hào): 570836
- 注冊(cè): 2008-06-08
- 專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
| 另外在程序中，我的理解是反拉普拉斯變換“InverseLaplaceTransform”這個(gè)函數(shù)是對(duì)字符變量進(jìn)行處理，得到的是一個(gè)m關(guān)于t的表達(dá)式。其實(shí)我也不需要明晰的表達(dá)式，因此要是有反拉普拉斯變換的數(shù)值解法也是好的。我知道反拉普拉斯變換應(yīng)該是一個(gè)積分運(yùn)算，但是具體什么公式就不曉得了，如果有哪位達(dá)人知道反拉普拉斯變換的公式也可以應(yīng)助阿，這樣我可以嘗試用matlab或mathematica的積分?jǐn)?shù)值解法進(jìn)行求解 |
10樓2011-07-26 09:46:47
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xhw2211 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277跪求調(diào)劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 309求調(diào)劑 +7 | 誰不是少年 2026-03-29 | 7/350 |
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +8 | Koxui 2026-03-28 | 8/400 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5 | SYA！ 2026-03-23 | 5/250 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|