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

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

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

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

加油大楠楠

銅蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 236.4
散金: 5
帖子: 99
在線: 26.9小時(shí)
蟲(chóng)號(hào): 1806539
注冊(cè): 2012-05-09
專(zhuān)業(yè): 機(jī)構(gòu)學(xué)與機(jī)器人

[求助] 請(qǐng)教matlab用數(shù)組給參數(shù)矩陣賦值的問(wèn)題

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為我需要的比矩陣
到這一步矩陣的計(jì)算都不成問(wèn)題
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行對(duì)應(yīng)的列來(lái)賦值的
e = norm(d,'fro'); %求范數(shù)
f(n,1)=e;%將范數(shù)記錄
n=n+1;
if n==101;%B數(shù)組有100行所以到一百終止
end

最后結(jié)果f數(shù)列中只有一個(gè)數(shù) ，請(qǐng)幫忙

回復(fù)此樓

» 猜你喜歡

一志愿西電085401數(shù)一英一299求調(diào)劑六級(jí)521 已經(jīng)有3人回復(fù)
生物考研337分求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑：085600材料與化工，考材科基，總分319 已經(jīng)有14人回復(fù)
285求調(diào)劑已經(jīng)有10人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有29人回復(fù)
理學(xué)07化學(xué) 303求調(diào)劑已經(jīng)有5人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有3人回復(fù)
土木304求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國(guó)科學(xué)院大學(xué)265求調(diào)劑已經(jīng)有7人回復(fù)
286求調(diào)劑已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

matlab稀疏矩陣問(wèn)題求教已經(jīng)有3人回復(fù)
如何用matlab求解矩陣系數(shù)的二階微分方程已經(jīng)有13人回復(fù)
matlab函數(shù)調(diào)用與數(shù)組的問(wèn)題已經(jīng)有19人回復(fù)
用matlab來(lái)實(shí)現(xiàn)分類(lèi)匯總已經(jīng)有3人回復(fù)
如何用Matlab畫(huà)矩陣的散點(diǎn)圖？已經(jīng)有3人回復(fù)
菜鳥(niǎo)求助~matlab怎么樣讓一個(gè)1xn的向量中的每個(gè)元素是個(gè)矩陣？已經(jīng)有11人回復(fù)
請(qǐng)教向量Q滿足高斯分布，用matlab如何表示？已經(jīng)有11人回復(fù)
請(qǐng)教關(guān)于matlab矩陣轉(zhuǎn)換問(wèn)題已經(jīng)有6人回復(fù)
matlab線性矩陣不等式求解做優(yōu)化問(wèn)題已經(jīng)有10人回復(fù)
如何用c++來(lái)求解矩陣的相關(guān)問(wèn)題。已經(jīng)有12人回復(fù)
【分享】Matlab經(jīng)典教程——從入門(mén)到精通【已搜索無(wú)重復(fù)】已經(jīng)有797人回復(fù)
VS2005下matlab調(diào)用已經(jīng)有8人回復(fù)
【求助】請(qǐng)教在在matlab里面怎么輸入含參數(shù)矩陣指數(shù)形式-----------【已完結(jié)】已經(jīng)有5人回復(fù)
【求助】關(guān)于matlab的高維數(shù)組已經(jīng)有6人回復(fù)
【求助】初學(xué)matlab 已經(jīng)有11人回復(fù)
【求助】柱形統(tǒng)計(jì)圖matlab怎么畫(huà)？已經(jīng)有4人回復(fù)
【求助】MATLAB的矩陣擴(kuò)展問(wèn)題已經(jīng)有2人回復(fù)

1樓 2013-04-15 09:14:09

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

材料廖

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 15 (小學(xué)生)
金幣: 4889.5
散金: 200
紅花: 11
帖子: 891
在線: 295.1小時(shí)
蟲(chóng)號(hào): 1794189
注冊(cè): 2012-05-04
性別: GG
專(zhuān)業(yè): 金屬材料的合金相、相變及

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

★ ★
jjdg: 金幣+2, 辛苦了 2013-04-22 10:54:37

引用回帖:

6樓: Originally posted by 加油大楠楠 at 2013-04-21 10:17:29
程序中的嗎，m,n就是xy,為了減少計(jì)算量我把z固定了是20，其實(shí)xyz是空間坐標(biāo)...

開(kāi)始那個(gè)由于你在rad中用了0，導(dǎo)致了出現(xiàn)了NAN，我照你的程序修改了下，可以出部分結(jié)果，程序如下：
j=1;z=20;r=32;
for l=32:34 %設(shè)定范圍
i=1;  tezheng=0;
for rad=1:l; %這兩個(gè)for的步長(zhǎng)可以根據(jù)運(yùn)算效率適當(dāng)放大，此處我改成了從1開(kāi)始
      for thea=0:pi/10:2*pi;

         x=rad*cos(thea);  %x
         y=rad*sin(thea);  %y
         d11=1/3*(3*r+3*x+sqrt(9*l^2-9*y^2-9*z^2));
         d12=1/3*(6*r-3*x+3*y*sqrt(3)+sqrt(36*l^2-27*x^2-18*x*y*sqrt(3)-9*y^2-36*z^2));
         d13=1/3*(6*r-3*x-3*y*sqrt(3)+sqrt(36*l^2-27*x^2+18*x*y*sqrt(3)-9*y^2-36*z^2));
         B(i,:)=[x,y,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/(eps+e);%把這個(gè)范數(shù)記錄在hh中
         tezheng=hh(i,1)+tezheng;%相加對(duì)于相同l值的對(duì)應(yīng)的范數(shù)
         i=i+1;
      end
end
tezheng1(j,:)=tezheng/(i-1+eps);%把不同l值對(duì)應(yīng)的范數(shù)求和再除以相應(yīng)的l值  存在數(shù)組tezheng1中
j=j+1;
end

根據(jù)你說(shuō)的和我理解的，我自己寫(xiě)了一個(gè)，運(yùn)行速度理論上會(huì)快點(diǎn)，也可以把你計(jì)算過(guò)程中的大部分?jǐn)?shù)據(jù)記錄下來(lái)，程序如下：

%對(duì)于l我感覺(jué)你要是后面是求和的話，沒(méi)有必要那樣，你的求和在多重循環(huán)里面，重復(fù)了很多次，我就將l設(shè)成了你程序中的最大值，這個(gè)你自己再考慮下，看著修改吧。

j=1;r=32;l = 34;  %注意我把l設(shè)成了一個(gè)常數(shù)
for rad = 1:34
%theta = 0:pi/10:2*pi
x(rad,:) = rad*cos(0:pi/10:2*pi);  %x
y(rad,:) = rad*sin(0:pi/10:2*pi);  %y
z(rad,:) = 20*ones(size(0:pi/10:2*pi));
d11(rad,:)=1/3*(3*r+3*x(rad,:)+sqrt(9*l^2-9*y(rad,:).^2-9*z(rad,:).^2));
d12(rad,:)=1/3*(6*r-3*x(rad,:)+3*y(rad,:)*sqrt(3)+sqrt(36*l^2-...
      27*x(rad,:).^2-18*x(rad,:).*y(rad,:)*sqrt(3)-9*y(rad,:).^2-36*z(rad,:).^2));
d13(rad,:)=1/3*(6*r-3*x(rad,:)-3*y(rad,:)*sqrt(3)+sqrt(36*l^2-...
      27*x(rad,:).^2+18*x(rad,:).*y(rad,:)*sqrt(3)-9*y(rad,:).^2-36*z(rad,:).^2));
B(rad,1:21,1:6) = [(x(rad,:))',(y(rad,:))',(z(rad,:))',(d11(rad,:))',(d12(rad,:))',(d13(rad,:))'];
n = length(B(rad,:,1));
for k = 1:n
      a(:,:,rad,k) = (1/B(rad,k,1)).*[B(rad,k,3)+r-B(rad,k,4) 0 0;0 -0.866.*B(rad,k,2)-0.5.*B(rad,k,3)+...
         r-B(rad,k,5) 0;0 0 0.866.*B(rad,k,2)-0.5.*B(rad,k,3)+r-B(rad,k,6)];
      b(:,:,rad,k) = [B(rad,k,1) B(rad,k,2) B(rad,k,3)+r-B(rad,k,4);B(rad,k,1) B(rad,k,2)-0.866.*r+...
         0.866.*B(rad,k,5) B(rad,k,3)-0.5.*r+0.5.*B(rad,k,5);B(rad,k,1) B(rad,k,2)+...
         0.866.*r-0.866.*B(rad,k,6) B(rad,k,3)-0.5.*r+0.5.*B(rad,k,6)];
      d(:,:,rad,k)=a(:,:,rad,k)\b(:,:,rad,k); %雅克比矩陣
      e(rad,k) = norm(d(:,:,rad,k),'fro');  %二范數(shù)
end
hh = 1./e;  %同你程序中的hh，只是此處為一個(gè)矩陣，即所有你的hh都在這里
end
tezheng = sum((sum(hh)));%所有的hh求和
%到這里，求和都做了，所有的e，hh都有數(shù)據(jù)在里面，后面我不是很明白，你自己看著修改吧。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

加油大楠楠

7樓2013-04-21 18:10:57

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

查看全部 12 個(gè)回答

材料廖

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 15 (小學(xué)生)
金幣: 4889.5
散金: 200
紅花: 11
帖子: 891
在線: 295.1小時(shí)
蟲(chóng)號(hào): 1794189
注冊(cè): 2012-05-04
性別: GG
專(zhuān)業(yè): 金屬材料的合金相、相變及

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

感謝參與，應(yīng)助指數(shù) +1

你這個(gè)肯定只有n=1的值呀，都沒(méi)有循環(huán)，你可以將你的這段程序
f(n,1)=e;%將范數(shù)記錄
n=n+1;
if n==101;%B數(shù)組有100行所以到一百終止
end
改成：
for n = 1:100
f(n,1) = e;
end
這樣應(yīng)該就可以了！

贊一下

回復(fù)此樓

2樓2013-04-15 12:12:41

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

加油大楠楠

銅蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 236.4
散金: 5
帖子: 99
在線: 26.9小時(shí)
蟲(chóng)號(hào): 1806539
注冊(cè): 2012-05-09
專(zhuān)業(yè): 機(jī)構(gòu)學(xué)與機(jī)器人

引用回帖:

2樓: Originally posted by 材料廖 at 2013-04-15 12:12:41
你這個(gè)肯定只有n=1的值呀，都沒(méi)有循環(huán)，你可以將你的這段程序
f(n,1)=e;%將范數(shù)記錄
n=n+1;
if n==101;%B數(shù)組有100行所以到一百終止
end
改成：
for n = 1:100
f(n,1) = e;
end
這樣應(yīng)該就可以了！

不好意思回復(fù)晚了  那你再幫我看看這整段程序我把上次的程序稍微改了下，就是對(duì)于每個(gè)l，我要求他xyx為不同值的時(shí)候的范數(shù)，然后相加，最后把每個(gè)l的范數(shù)和提取出來(lái)，十分感謝了  我運(yùn)行總是報(bào)錯(cuò)
j=1;z=20;r=32;
for l=32:34 %設(shè)定范圍
i=1;  tezheng=0;
for rad=0:l; %這兩個(gè)for的步長(zhǎng)可以根據(jù)運(yùn)算效率適當(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;%把這個(gè)范數(shù)記錄在hh中
tezheng=hh(i,1)+tezheng;%相加對(duì)于相同l值的對(duì)應(yīng)的范數(shù)
i=i+1;
end
end
tezheng1(j,

=tezheng/(i-1);%把不同l值對(duì)應(yīng)的范數(shù)求和再除以相應(yīng)的l值存在數(shù)組tezheng1中
j=j+1;
end

錯(cuò)誤信息是
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];

>>

贊一下

回復(fù)此樓

3樓2013-04-18 09:38:21

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

材料廖

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 15 (小學(xué)生)
金幣: 4889.5
散金: 200
紅花: 11
帖子: 891
在線: 295.1小時(shí)
蟲(chóng)號(hào): 1794189
注冊(cè): 2012-05-04
性別: GG
專(zhuān)業(yè): 金屬材料的合金相、相變及

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

我還是不是很明白你要做什么，你這個(gè)程序?qū)懙奈矣悬c(diǎn)看不懂。但有點(diǎn)建議，一：你這個(gè)應(yīng)該沒(méi)有必要用符號(hào)計(jì)算，符號(hào)計(jì)算很慢。二：你這個(gè)程序編寫(xiě)的不怎樣，基礎(chǔ)不是很好，寫(xiě)的有C語(yǔ)言的感覺(jué)，可能你的想法你沒(méi)用程序表達(dá)出來(lái)，要不你說(shuō)說(shuō)你要干啥吧，我看能不能幫你。話說(shuō)直了點(diǎn)，望不怪。
下面是我沒(méi)用符號(hào)計(jì)算程序，速度快了很多，但是還是存在錯(cuò)誤，出不了結(jié)果，因?yàn)槲疫€是沒(méi)有理解你要干啥。
j=1;z=20;r=32;
for l=32:34 %設(shè)定范圍
i=1;  tezheng=0;
for rad=0:l; %這兩個(gè)for的步長(zhǎng)可以根據(jù)運(yùn)算效率適當(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;%把這個(gè)范數(shù)記錄在hh中
         tezheng=hh(i,1)+tezheng;%相加對(duì)于相同l值的對(duì)應(yīng)的范數(shù)
         i=i+1;
      end
end
tezheng1(j,

=tezheng/(i-1+eps);%把不同l值對(duì)應(yīng)的范數(shù)求和再除以相應(yīng)的l值存在數(shù)組tezheng1中
j=j+1;
end

贊一下

回復(fù)此樓

4樓2013-04-19 18:36:16

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

查看全部 12 個(gè)回答

普通表情龍兔虎貓

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

[考研] 化學(xué)工程085602 305分求調(diào)劑 +27	RichLi_ 2026-03-25	29/1450	2026-03-31 11:33 by winsuccess
[考研] 287求調(diào)劑 +17	land xuxu 2026-03-26	17/850	2026-03-31 11:16 by Zzxxxs
[考研] 0856調(diào)劑 +6	曲聽(tīng)筠 2026-03-30	6/300	2026-03-31 09:44 by Jaylen.
[考研] 一志愿食品科學(xué)與工程083200求調(diào)劑 +4	XQTJZ 2026-03-30	4/200	2026-03-31 04:10 by fmesaito
[考研] 福建理工大學(xué)材料學(xué)院先進(jìn)合金團(tuán)隊(duì)招收考研調(diào)劑學(xué)生 +3	大華金商都 2026-03-30	4/200	2026-03-31 01:04 by 方英俊602
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 材料工程專(zhuān)碩求調(diào)劑 +8	hyl3153942 2026-03-29	8/400	2026-03-30 20:45 by dophin1985
[考研] 332求調(diào)劑 +14	032500 2026-03-25	14/700	2026-03-30 20:03 by 無(wú)際的草原
[考研] 一志愿廈門(mén)大學(xué)材料工程專(zhuān)碩354找調(diào)劑�。�！ +5	貝唄鋇鋇 2026-03-30	5/250	2026-03-30 18:16 by 無(wú)際的草原
[考研] 085600，材料與化工321分求調(diào)劑 +10	大饞小子 2026-03-28	10/500	2026-03-29 23:35 by 飛行日記西
[考研] 275求調(diào)劑 +15	Micky11223 2026-03-25	20/1000	2026-03-29 20:44 by 唐沐兒
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8	柳溪邊 2026-03-26	8/400	2026-03-29 20:42 by 唐沐兒
[考研] 086000生物與醫(yī)藥調(diào)劑 +5	Feisty。 2026-03-28	9/450	2026-03-29 12:02 by longlotian
[考研] 085600，專(zhuān)業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 275求調(diào)劑 +10	jjjjjjjjjjl 2026-03-27	10/500	2026-03-27 23:47 by barnett0632
[考研] 0856調(diào)劑 +5	求求讓我有書(shū)讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 環(huán)境專(zhuān)碩324分求調(diào)劑推薦 +5	軒小寧—— 2026-03-26	5/250	2026-03-26 12:05 by i_cooler