版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請盡快對帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

fabaoyi

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 28
帖子: 11
在線: 8.6小時(shí)
蟲號(hào): 2066628
注冊: 2012-10-17
性別: GG
專業(yè): 應(yīng)用數(shù)學(xué)方法

【答案】應(yīng)助回帖

引用回帖:

8樓: Originally posted by wangjunyahit at 2012-11-25 17:04:02
看不懂啊，能不能給下面的語句加點(diǎn)注釋
=disha(B);             %??
x=(B(:,1:n))\B(:,n+1)    %??
s1=A(i,;                      %??
            A(i,=A(row+i-1,;
            A(row+i ...

表情不是我加的,是一些字符自動(dòng)轉(zhuǎn)成了表情

贊一下

回復(fù)此樓

活著為了實(shí)現(xiàn)價(jià)值

11樓2012-11-25 21:57:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zsy312

銀蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: 495.5
紅花: 1
帖子: 145
在線: 20.7小時(shí)
蟲號(hào): 756330
注冊: 2009-04-24
性別: GG
專業(yè): 數(shù)學(xué)

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與！ 2013-04-27 23:43:07

這種線性方程組直接有符號(hào)解, 還可以比較和數(shù)值解之間的誤差

CODE:

A=[-1.2e-21,1.4e6,1.4e6,-1.4e6,1.3e6;

    8.9e-22,-1.4e6,4.9e3,3.8e3,0;

    1.8e-22,1.4e3,-1.4e6,4.3e3,0;

    4.4e-23,3.2e2,1.2e3,-5.1e6,4.3e3;

    5.3e-23,0,0,1.3e3,-2.3e6];

b=[0;-1.8e-4;-1.3e-2;-1.5e-2;-4.0e-3];

x=A\b;

A1=sym(A);

b1=sym(b);

x1=A1\b;

disp(['數(shù)值解: '])

disp(num2str(x))

disp(['符號(hào)解: '])

disp(x1)

disp(['符號(hào)解近似值: '])

disp(num2str(double(x1)))

disp(' ')

disp(['二者的誤差為: '])

disp(num2str(x-double(x1)))

贊一下

回復(fù)此樓

12樓2012-11-25 22:08:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zsy312

銀蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: 495.5
紅花: 1
帖子: 145
在線: 20.7小時(shí)
蟲號(hào): 756330
注冊: 2009-04-24
性別: GG
專業(yè): 數(shù)學(xué)

輸出:

CODE:

數(shù)值解: 

1.029393177754837e+020

6.565808888911405e-008

 2.25982285154494e-008

3.842185487752932e-009

4.113382036189007e-009

符號(hào)解: 

 388474900292779720441528931100854424320492441698304/3773824314049418181639511136225

              1776102037346606295811692690717/27050772683106229525992015824460800000

             3820622141326942446909828077309/169067329269413934537450098902880000000

                  7423863302538155608170509459/1932198048793302108999429701747200000

                  6409571567956304127562009059/1558224232897824281451152985280000000

符號(hào)解近似值: 

1.029393177754837e+020

6.565808888911403e-008

 2.25982285154494e-008

3.842185487752931e-009

4.113382036189006e-009

二者的誤差為: 

          32768

2.64697796e-023

3.30872245e-024

8.27180613e-025

1.65436123e-024

贊一下

回復(fù)此樓

13樓2012-11-25 22:11:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

barryyee

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 90.8
散金: 180
紅花: 3
帖子: 293
在線: 100.5小時(shí)
蟲號(hào): 1612694
注冊: 2012-02-12
專業(yè): 敏感電子學(xué)與傳感器

引用回帖:

3樓: Originally posted by fabaoyi at 2012-11-25 01:28:13
你學(xué)過數(shù)值分析沒有,這個(gè)可以用高斯迭代法和雅可比迭代法求解,不過雖然用迭代法,跑程序時(shí)還是有警告,可能求解還是不是非常可靠,我求解了一個(gè),下面我把高斯迭代法那個(gè)我求解用的程序和運(yùn)行結(jié)果貼在下面:
function g ...

matlab 解線性方程組默認(rèn)是高斯法，所以結(jié)果是一樣的。

贊一下

回復(fù)此樓

14樓2013-04-27 14:31:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

change0618

鐵桿木蟲 (著名寫手)

方丈大師

應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
紅花: 17
帖子: 2413
在線: 546.7小時(shí)
蟲號(hào): 496517
注冊: 2008-01-19
專業(yè): 化學(xué)反應(yīng)工程

【答案】應(yīng)助回帖

★
xiegangmai: 金幣+1, 謝謝參與！ 2013-04-27 23:43:15

對于剛性的矩陣，采用x=pinv(A)*b 就可以了

贊一下(1人)

回復(fù)此樓

15樓2013-04-27 15:35:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

change0618

鐵桿木蟲 (著名寫手)

方丈大師

應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
紅花: 17
帖子: 2413
在線: 546.7小時(shí)
蟲號(hào): 496517
注冊: 2008-01-19
專業(yè): 化學(xué)反應(yīng)工程

【答案】應(yīng)助回帖

★ ★
xiegangmai: 金幣+2, 謝謝參與！ 2013-04-27 23:43:24

CODE:

>> A=[-1.2e-21,1.4e6,1.4e6,-1.4e6,1.3e6;

     8.9e-22,-1.4e6,4.9e3,3.8e3,0;

     1.8e-22,1.4e3,-1.4e6,4.3e3,0;

     4.4e-23,3.2e2,1.2e3,-5.1e6,4.3e3;

     5.3e-23,0,0,1.3e3,-2.3e6];

 b=[0;-1.8e-4;-1.3e-2;-1.5e-2;-4.0e-3];

>> x=pinv(A)*b

x =

   1.0e-08 *

   0.000000000000000

  -0.223269102876553

   0.689363755668183

   0.312251529518066

   0.091881357445947

>> vpa(x,30)

ans =

       3.70460811145685954277266468961e-37

  -0.0000000022326910287655349973263795146

   0.0000000068936375566818299856437653433

  0.00000000312251529518065781430357617809

 0.000000000918813574459469995008142319883

贊一下

回復(fù)此樓

16樓2013-04-27 15:45:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

change0618

鐵桿木蟲 (著名寫手)

方丈大師

應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
紅花: 17
帖子: 2413
在線: 546.7小時(shí)
蟲號(hào): 496517
注冊: 2008-01-19
專業(yè): 化學(xué)反應(yīng)工程

【答案】應(yīng)助回帖

★ ★ ★ ★
xiegangmai: 金幣+2, 謝謝參與！ 2013-04-27 23:43:35
wangjunyahit: 金幣+2, ★★★很有幫助 2013-04-30 16:19:30

CODE:

>> A=[-1.070906744447462E-017   244564571.598802        303866255.019473      303856255.019473        303756255.019473  ;  

    9.735515836938967E-017   -437040200.438244        546168396.621813        415770404.263350      0.000000000000000E+000 ; 

    1.947103215071310E-017   157778358.747585     -986775896.684545        473345946.492110       0.000000000000000E+000; 

    4.867758037678274E-018   34697270.0918581        136741245.043259       -1333413850.81819        473345946.492110   ;   

     5.841309406796348E-018  0.000000000000000E+000  0.000000000000000E+000   136741245.043259      -778102201.511583 ];

 b=[ 104348331231.704       -2072871.66710357      -154296552.906917       -178403143.022063       -48501206.8009638 ]';

>> x=pinv(A)*b

x =

   1.0e+02 *

  -0.000000000000000

   1.682414884915961

   0.668383858892869

   0.388902999746503

   0.336174005405975

>> vpa(x,30)

ans =

 -0.0000000000000000000000285841068286876334276912171331

                         168.241488491596101084724068642

                         66.8383858892869398005132097751

                         38.8902999746503041933465283364

                         33.6174005405974725135820335709

贊一下(1人)

回復(fù)此樓

17樓2013-04-27 15:48:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

diqiuyazhou

木蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: 1828.8
散金: 1045
紅花: 1
帖子: 106
在線: 96.4小時(shí)
蟲號(hào): 1089840
注冊: 2010-09-05
性別: GG
專業(yè): 聲學(xué)

想問一下樓主，問題最后解決了嗎，我最近也遇到這類問題，不知怎么辦

贊一下

回復(fù)此樓

18樓2015-04-22 23:37:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 wangjunyahit 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4	share_joy 2026-03-16	4/200	2026-03-17 21:32 by hmn_wj
[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 312求調(diào)劑 +4	陌宸希 2026-03-16	5/250	2026-03-17 17:09 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7	Ncdx123456 2026-03-13	8/400	2026-03-16 12:15 by karry wen
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 材料工程調(diào)劑 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 工科調(diào)劑 +4	Jiang191123！ 2026-03-11	4/200	2026-03-13 15:15 by Miko19
[考研] 085600材料與化工 309分請求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant