 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2558 | 回復(fù): 11 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
加油大楠楠銅蟲 (小有名氣)
|
[求助]
請教matlab用數(shù)組給參數(shù)矩陣賦值的問題
|
|
|
syms x y z d1 d2 d3 r l;n=1;%定義矩陣參數(shù) a=1/l*[z+r-d1 0 0;0 -0.866*y-0.5*z+r-d2 0;0 0 0.866*y-0.5*z+r-d3]; b=[x y z+r-d1;x y-0.866*r+0.866*d2 z-0.5*r+0.5*d2;x y+0.866*r-0.866*d3 z-0.5*r+0.5*d3]; c=inv(a); d=c*b;%d為我需要的比矩陣 到這一步矩陣的計算都不成問題 d=subs(d,{'r','l','x','y','z','d1','d2','d3'},{31.8,40.2,B(n,1),B(n,2),B(n,3),B(n,4),B(n,5),B(n,6)}); %給參數(shù)賦值，其中x y z d1 d2 d3是用B矩陣的第n行對應(yīng)的列來賦值的 e = norm(d,'fro'); %求范數(shù) f(n,1)=e;%將范數(shù)記錄 n=n+1; if n==101;%B數(shù)組有100行所以到一百終止 end 最后結(jié)果f數(shù)列中只有一個數(shù) ，請幫忙 |
回復(fù)此樓» 猜你喜歡
305分求調(diào)劑（食品工程）
已經(jīng)有7人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有5人回復(fù)
-
299求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料292調(diào)劑
已經(jīng)有5人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有4人回復(fù)
-
收08調(diào)劑生
已經(jīng)有3人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有9人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
086003食品工程求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856材料專碩353求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab稀疏矩陣問題求教
已經(jīng)有3人回復(fù)
- 如何用matlab求解矩陣系數(shù)的二階微分方程
已經(jīng)有13人回復(fù)
- matlab函數(shù)調(diào)用與數(shù)組的問題
已經(jīng)有19人回復(fù)
- 用matlab來實現(xiàn)分類匯總
已經(jīng)有3人回復(fù)
- 如何用Matlab畫矩陣的散點圖？
已經(jīng)有3人回復(fù)
- 菜鳥求助~matlab怎么樣讓一個1xn的向量中的每個元素是個矩陣？
已經(jīng)有11人回復(fù)
- 請教向量Q滿足高斯分布，用matlab如何表示？
已經(jīng)有11人回復(fù)
- 請教關(guān)于matlab矩陣轉(zhuǎn)換問題
已經(jīng)有6人回復(fù)
- matlab線性矩陣不等式求解做優(yōu)化問題
已經(jīng)有10人回復(fù)
- 如何用c++來求解矩陣的相關(guān)問題。
已經(jīng)有12人回復(fù)
- 【分享】Matlab經(jīng)典教程——從入門到精通【已搜索無重復(fù)】
已經(jīng)有797人回復(fù)
- VS2005下matlab調(diào)用
已經(jīng)有8人回復(fù)
- 【求助】請教在在matlab里面怎么輸入含參數(shù)矩陣指數(shù)形式-----------【已完結(jié)】
已經(jīng)有5人回復(fù)
- 【求助】關(guān)于matlab的高維數(shù)組
已經(jīng)有6人回復(fù)
- 【求助】初學(xué)matlab
已經(jīng)有11人回復(fù)
- 【求助】柱形統(tǒng)計圖matlab怎么畫？
已經(jīng)有4人回復(fù)
- 【求助】MATLAB的矩陣擴(kuò)展問題
已經(jīng)有2人回復(fù)
材料廖
木蟲 (正式寫手)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 4889.5
- 散金: 200
- 紅花: 11
- 帖子: 891
- 在線: 295.1小時
- 蟲號: 1794189
- 注冊: 2012-05-04
- 性別: GG
- 專業(yè): 金屬材料的合金相、相變及
9樓2013-04-22 22:01:14
材料廖
木蟲 (正式寫手)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 4889.5
- 散金: 200
- 紅花: 11
- 帖子: 891
- 在線: 295.1小時
- 蟲號: 1794189
- 注冊: 2012-05-04
- 性別: GG
- 專業(yè): 金屬材料的合金相、相變及
2樓2013-04-15 12:12:41
加油大楠楠
銅蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 236.4
- 散金: 5
- 帖子: 99
- 在線: 26.9小時
- 蟲號: 1806539
- 注冊: 2012-05-09
- 專業(yè): 機(jī)構(gòu)學(xué)與機(jī)器人
|
不好意思回復(fù)晚了 那你再幫我看看這整段程序 我把上次的程序稍微改了下，就是對于每個l，我要求他xyx為不同值的時候的范數(shù)，然后相加，最后把每個l的范數(shù)和提取出來，十分感謝了 我運行總是報錯 j=1;z=20;r=32; for l=32:34 %設(shè)定范圍 i=1; tezheng=0; for rad=0:l; %這兩個for的步長可以根據(jù)運算效率適當(dāng)放大 for thea=0:pi/10:2*pi; m=rad*cos(thea); n=rad*sin(thea); d11=1/3*(3*r+3*m+sqrt(9*l^2-9*n^2-9*z^2)); d12=1/3*(6*r-3*m+3*n*sqrt(3)+sqrt(36*l^2-27*m^2-18*m*n*sqrt(3)-9*n^2-36*z^2)); d13=1/3*(6*r-3*m-3*n*sqrt(3)+sqrt(36*l^2-27*m^2+18*m*n*sqrt(3)-9*n^2-36*z^2)); B(i,  =[m,n,z,d11,d12,d13];syms x1 y1 z1 d1 d2 d3 l1; a=(1/l1)*[z1+r-d1 0 0;0 -0.866*y1-0.5*z1+r-d2 0;0 0 0.866*y1-0.5*z1+r-d3]; b=[x1 y1 z1+r-d1;x1 y1-0.866*r+0.866*d2 z1-0.5*r+0.5*d2;x1 y1+0.866*r-0.866*d3 z1-0.5*r+0.5*d3]; c=inv(a); d=c*b;%d為雅克比矩陣 d=subs(d,{'l1','x1','y1','z1','d1','d2','d3'},{l,B(i,1),B(i,2),B(i,3),B(i,4),B(i,5),B(i,6)}); e = norm(d,'fro'); %2范數(shù) hh(i,1)=1/e;%把這個范數(shù)記錄在hh中 tezheng=hh(i,1)+tezheng;%相加對于相同l值的對應(yīng)的范數(shù) i=i+1; end end tezheng1(j, =tezheng/(i-1);%把不同l值對應(yīng)的范數(shù)求和再除以相應(yīng)的l值 存在數(shù)組tezheng1中j=j+1; end 錯誤信息是 clear ??? Error using ==> mupadmex Error in MuPAD command: Out of memory Error in ==> sym.horzcat at 25 y = mupadmex('mllib::horzcat',strs{:}); Error in ==> jacobi at 16 a=(1/l1)*[z1+r-d1 0 0;0 -0.866*y1-0.5*z1+r-d2 0;0 0 0.866*y1-0.5*z1+r-d3]; >> |
3樓2013-04-18 09:38:21
材料廖
木蟲 (正式寫手)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 4889.5
- 散金: 200
- 紅花: 11
- 帖子: 891
- 在線: 295.1小時
- 蟲號: 1794189
- 注冊: 2012-05-04
- 性別: GG
- 專業(yè): 金屬材料的合金相、相變及
【答案】應(yīng)助回帖
|
我還是不是很明白你要做什么，你這個程序?qū)懙奈矣悬c看不懂。但有點建議，一：你這個應(yīng)該沒有必要用符號計算，符號計算很慢。二：你這個程序編寫的不怎樣，基礎(chǔ)不是很好，寫的有C語言的感覺，可能你的想法你沒用程序表達(dá)出來，要不你說說你要干啥吧，我看能不能幫你。話說直了點，望不怪。 下面是我沒用符號計算程序，速度快了很多，但是還是存在錯誤，出不了結(jié)果，因為我還是沒有理解你要干啥。 j=1;z=20;r=32; for l=32:34 %設(shè)定范圍 i=1; tezheng=0; for rad=0:l; %這兩個for的步長可以根據(jù)運算效率適當(dāng)放大 for thea=0:pi/10:2*pi; m=rad*cos(thea); n=rad*sin(thea); d11=1/3*(3*r+3*m+sqrt(9*l^2-9*n^2-9*z^2)); d12=1/3*(6*r-3*m+3*n*sqrt(3)+sqrt(36*l^2-27*m^2-18*m*n*sqrt(3)-9*n^2-36*z^2)); d13=1/3*(6*r-3*m-3*n*sqrt(3)+sqrt(36*l^2-27*m^2+18*m*n*sqrt(3)-9*n^2-36*z^2)); B(i, =[m,n,z,d11,d12,d13];a=(1/B(i,1))*[B(i,3)+r-B(i,4) 0 0;0 -0.866*B(i,2)-0.5*B(i,3)+r-B(i,5) 0;0 0 0.866*B(i,2)-0.5*B(i,3)+r-B(i,6)]; b=[B(i,1) B(i,2) B(i,3)+r-B(i,4);B(i,1) B(i,2)-0.866*r+0.866*B(i,5) B(i,3)-0.5*r+0.5*B(i,5);B(i,1) B(i,2)+0.866*r-0.866*B(i,6) B(i,3)-0.5*r+0.5*B(i,6)]; % c=inv(a); d=a\b;%d為雅克比矩陣 e = norm(d,'fro'); %2范數(shù) hh(i,1)=1/e;%把這個范數(shù)記錄在hh中 tezheng=hh(i,1)+tezheng;%相加對于相同l值的對應(yīng)的范數(shù) i=i+1; end end tezheng1(j, =tezheng/(i-1+eps);%把不同l值對應(yīng)的范數(shù)求和再除以相應(yīng)的l值 存在數(shù)組tezheng1中j=j+1; end |
4樓2013-04-19 18:36:16
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +11 | 手機(jī)用戶 2026-03-17 | 12/600 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +11 | hanamiko 2026-03-18 | 11/550 |
|
|
|
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 3/150 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|