 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 1861 | 回復(fù): 9 | ||
xhw2211金蟲 (正式寫手)
|
[求助]
mathematica 編程求助
|
|
頭一次用mathtematica編的程序，總是沒有畫出圖形，不知道錯在什么地方了，大�？匆幌�。 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ù)此樓» 猜你喜歡
一志愿南開大學0710生物學359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
327求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600 286分 材料求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿河北工業(yè)大學0817化工278分求調(diào)劑
已經(jīng)有15人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有15人回復(fù)
-
南京大學化學調(diào)劑
已經(jīng)有14人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求mathematica編程
已經(jīng)有6人回復(fù)
- 【求助】開平方問題【已完結(jié)】
已經(jīng)有3人回復(fù)
- 【求助】求解氫原子電子密度的時候遇到的4個問題
已經(jīng)有4人回復(fù)
- 【求助】如何利用matlab計算這個無窮級數(shù)的和
已經(jīng)有12人回復(fù)
- 【求助】100金幣求助高人一個程序解方程
已經(jīng)有15人回復(fù)
- 【求助】有關(guān)fortran進行數(shù)值模擬
已經(jīng)有10人回復(fù)
jjdg
版主 (知名作家)
- 程序強帖: 9
- 應(yīng)助: 223 (大學生)
- 貴賓: 3.926
- 金幣: 51005.2
- 紅花: 238
- 帖子: 7227
- 在線: 474.2小時
- 蟲號: 27000
- 注冊: 2003-10-21
- 專業(yè): 藥物代謝與藥物動力學
- 管轄: 程序語言
2樓2011-07-25 17:51:26
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時
- 蟲號: 570836
- 注冊: 2008-06-08
- 專業(yè): 機械結(jié)構(gòu)強度學
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小時
- 蟲號: 570836
- 注冊: 2008-06-08
- 專業(yè): 機械結(jié)構(gòu)強度學
4樓2011-07-25 18:35:01
mshwangg
至尊木蟲 (正式寫手)
- 程序強帖: 5
- 應(yīng)助: 206 (大學生)
- 金幣: 10702.8
- 散金: 100
- 紅花: 19
- 帖子: 597
- 在線: 195.4小時
- 蟲號: 576702
- 注冊: 2008-06-21
- 專業(yè): 物理學I
5樓2011-07-25 21:47:56
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時
- 蟲號: 570836
- 注冊: 2008-06-08
- 專業(yè): 機械結(jié)構(gòu)強度學
|
我按你的方法改了程序，但是還是沒有結(jié)果，感覺是不是拉普拉斯變換后解出的Z表達式太復(fù)雜，mathematica無法進行拉普拉斯反變換？？？下面是程序運行后得到的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小時
- 蟲號: 570836
- 注冊: 2008-06-08
- 專業(yè): 機械結(jié)構(gòu)強度學
7樓2011-07-26 09:26:19
xhw2211
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9240.1
- 紅花: 1
- 帖子: 301
- 在線: 56.3小時
- 蟲號: 570836
- 注冊: 2008-06-08
- 專業(yè): 機械結(jié)構(gòu)強度學
|
把修改后的程序貼出來： 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小時
- 蟲號: 570836
- 注冊: 2008-06-08
- 專業(yè): 機械結(jié)構(gòu)強度學
★
jjdg(金幣+1): 感謝參與 2011-07-30 14:50:08
jjdg(金幣+1): 感謝參與 2011-07-30 14:50:08
|
下面這個是對的： 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小時
- 蟲號: 570836
- 注冊: 2008-06-08
- 專業(yè): 機械結(jié)構(gòu)強度學
10樓2011-07-26 09:46:47
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xhw2211 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 南京大學化學調(diào)劑 +9 | 景隨風 2026-03-29 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 材料科學與工程求調(diào)劑 +6 | 深V宿舍吧 2026-03-29 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
|
[考研] 0703化學 +11 | 妮妮ninicgb 2026-03-27 | 11/550 |
|
|
|
[考研] 071000生物學求調(diào)劑，初試成績343 +7 | 小小甜面團 2026-03-25 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
|
[碩博家園] 北京林業(yè)大學碩導招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|