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

登錄注冊

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

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

返回列表

521438

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 82.9
帖子: 19
在線: 13.2小時(shí)
蟲號(hào): 1801287
注冊: 2012-05-07
性別: MM
專業(yè): 電磁場與電路

引用回帖:

10樓: Originally posted by xiegangmai at 2012-05-09 12:31:19:
你的函數(shù)里包含符號(hào)I1,I2,肯定解不出數(shù)值解的。
不用簡化，把數(shù)學(xué)模型貼出來看看。

我的函數(shù)表達(dá)式如下：B=I*r^2/2*(r^2+z^2)^(-3/2)；
方程組分別是:
I1*d3Bdz3(z=a1)+I2*d3Bdz3(z=a2)=0
I1*d5Bdz5(z=a1)+I2*d5Bdz5(z=a2)=0
I1*d7Bdz7(z=a)+I2*d7Bdz7(z=a2)=0備注：這里的d3Bdz3，d5Bdz5，d7Bdz7是B對(duì)z求3，5，7階導(dǎo)數(shù)。然后我程序里面的式子就是分別求出3，5，7階導(dǎo)數(shù)后的表達(dá)式！版主幫我看看怎么修改程序��？謝謝！

贊一下

回復(fù)此樓

11樓2012-05-09 16:33:15

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

wjx0162

銀蟲 (著名寫手)

應(yīng)助: 3 (幼兒園)
金幣: 468.9
紅花: 1
帖子: 1152
在線: 68.4小時(shí)
蟲號(hào): 1187390
注冊: 2011-01-11
專業(yè): 機(jī)械動(dòng)力學(xué)

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

★
感謝參與，應(yīng)助指數(shù) +1
臭水溝: 金幣+1, 謝謝交流~~ 2012-05-10 12:05:02

MATLAB求解參數(shù)方程基本是忽悠人，建議你把參數(shù)代入方程，直接用他的非線性方程組求解方法來解。因?yàn)榇蠖喾蔷€性方程組沒有解析解。

贊一下

回復(fù)此樓

12樓2012-05-10 11:30:37

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

521438

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 82.9
帖子: 19
在線: 13.2小時(shí)
蟲號(hào): 1801287
注冊: 2012-05-07
性別: MM
專業(yè): 電磁場與電路

引用回帖:

12樓: Originally posted by wjx0162 at 2012-05-10 11:30:37:
MATLAB求解參數(shù)方程基本是忽悠人，建議你把參數(shù)代入方程，直接用他的非線性方程組求解方法來解。因?yàn)榇蠖喾蔷€性方程組沒有解析解。

我是把參數(shù)帶入求解了，是不是我的初值設(shè)置的有問題��！fsolve函數(shù)需要給定初值，大概是多少我也不知道，有沒有什么更好的辦法啊？

贊一下

回復(fù)此樓

13樓2012-05-10 11:36:04

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

xiegangmai

版主 (職業(yè)作家)

我沒頭銜

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
散金: 2394
紅花: 63
沙發(fā): 4
帖子: 4046
在線: 1330.2小時(shí)
蟲號(hào): 653607
注冊: 2008-11-13
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 仿真模擬

引用回帖:

13樓: Originally posted by 521438 at 2012-05-10 11:36:04:
我是把參數(shù)帶入求解了，是不是我的初值設(shè)置的有問題��！fsolve函數(shù)需要給定初值，大概是多少我也不知道，有沒有什么更好的辦法��？

多數(shù)情況下是不知道初值的，只有多試幾次，解出數(shù)值解了，再把這組解作為初值，代入后求解，多迭代幾次，再算算誤差。

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

明德厚學(xué)、求是創(chuàng)新

14樓2012-05-10 15:23:19

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

521438

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 82.9
帖子: 19
在線: 13.2小時(shí)
蟲號(hào): 1801287
注冊: 2012-05-07
性別: MM
專業(yè): 電磁場與電路

這里我想給表達(dá)式里的符號(hào)變量賦予初值，出現(xiàn)如下錯(cuò)誤：??? Conversion to double from sym is not possible.

Error in ==> sfdnls at 95
J(:,k) = (v-valx)/alpha(k);

Error in ==> optim\private\trustnleqn at 108
[JACfindiff,numFDfevals] = sfdnls(x,Fvec,JACfindiff,group,[], ...

Error in ==> fsolve at 295
[x,FVAL,JACOB,EXITFLAG,OUTPUT,msg]=...

Error in ==> fsolveDemon1 at 4
x=fsolve(@f,x0,opt);
就是不能把符號(hào)型變量轉(zhuǎn)換成double型的，大家知道我應(yīng)該怎么賦值么？需要如何修改�。�

贊一下

回復(fù)此樓

15樓2012-05-10 16:19:14

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

521438

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 82.9
帖子: 19
在線: 13.2小時(shí)
蟲號(hào): 1801287
注冊: 2012-05-07
性別: MM
專業(yè): 電磁場與電路

我的程序如下：
function x=fsolveDemon1
x0=[1.0;0.25;0.4];
opt=optimset('Display','off')
x=fsolve(@f,x0,opt);
function F=f(x)
syms I1 I2 r1 r2 a1 a2
r1=x(1);
r2=x(2);
a1=x(3);
a2=x(3);
F=[I1*r1^2*(3*a1*(r1^2+a1^2)^(-7/2)-7*a1^3*(r1^2+a1^2)^(-9/2))+I2*r2^2*(3*a2*(r2^2+a2^2)^(-7/2)-7*a2^3*(r2^2+a2^2)^(-9/2));
I1*r1^2*(-5*a1*(r1^2+a1^2)^(-9/2)+30*a1^3*(r1^2+a1^2)^(-11/2)-33*a1^5*(r1^2+a1^2)^(-13/2))+I2*r2^2*(-5*a2*(r2^2+a2^2)^(-9/2)+30*a2^3*(r2^2+a2^2)^(-11/2)-33*a2^5*(r2^2+a2^2)^(-13/2));
I1*r1^2*(35*a1*(r1^2+a1^2)^(-11/2)-385*a1^3*(r1^2+a1^2)^(-13/2)+1001*a1^5*(r1^2+a1^2)^(-15/2)-715*a1^7*(r1^2+a1^2)^(-17/2))+I2*r2^2*(35*a2*(r2^2+a2^2)^(-11/2)-385*a2^3*(r2^2+a2^2)^(-13/2)+1001*a2^5*(r2^2+a2^2)^(-15/2)-715*a2^7*(r2^2+a2^2)^(-17/2))];
end
end

贊一下

回復(fù)此樓

16樓2012-05-10 16:21:18

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

xiegangmai

版主 (職業(yè)作家)

我沒頭銜

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
散金: 2394
紅花: 63
沙發(fā): 4
帖子: 4046
在線: 1330.2小時(shí)
蟲號(hào): 653607
注冊: 2008-11-13
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 仿真模擬

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

引用回帖:

16樓: Originally posted by 521438 at 2012-05-10 16:21:18:
我的程序如下：
function x=fsolveDemon1
x0=;
opt=optimset('Display','off')
x=fsolve(@f,x0,opt);
function F=f(x)
syms I1 I2 r1 r2 a1 a2
r1=x(1);
r2=x(2);
a1=x(3);
a2=x(3);
F=;
end
end

你這樣的程序是根本求解不了的。
fsolve函數(shù)是求解非線性方程組數(shù)值解的，是不能帶符號(hào)的。
B=I*r^2/2*(r^2+z^2)^(-3/2)；這里面的I與I1、I2是什么關(guān)系？
把d3Bdz3，d5Bdz5，d7Bdz7寫出來，分別帶入
I1*d3Bdz3(z=a1)+I2*d3Bdz3(z=a2)=0
I1*d5Bdz5(z=a1)+I2*d5Bdz5(z=a2)=0
I1*d7Bdz7(z=a)+I2*d7Bdz7(z=a2)=0
三個(gè)方程，并寫入f(x)中，I1、I2、r作為x(4)、x(5)、x(6)，即F內(nèi)六個(gè)表達(dá)式。

你先試試吧，不行了把數(shù)學(xué)模型用圖貼出來吧，把已知條件和未知條件描述清楚。

贊一下

回復(fù)此樓

明德厚學(xué)、求是創(chuàng)新

17樓2012-05-10 23:25:01

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

521438

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 82.9
帖子: 19
在線: 13.2小時(shí)
蟲號(hào): 1801287
注冊: 2012-05-07
性別: MM
專業(yè): 電磁場與電路

引用回帖:

17樓: Originally posted by xiegangmai at 2012-05-10 23:25:01:
你這樣的程序是根本求解不了的。
fsolve函數(shù)是求解非線性方程組數(shù)值解的，是不能帶符號(hào)的。
B=I*r^2/2*(r^2+z^2)^(-3/2)；這里面的I與I1、I2是什么關(guān)系？
把d3Bdz3，d5Bdz5，d7Bdz7寫出來，分別帶入
I1*d3Bd ...

這里的I是電流，我的模型是4個(gè)圓線圈對(duì)稱放置組成而成。所帶電流分別為I1,I2。這里的數(shù)學(xué)表達(dá)式B=I*r^2/2*(r^2+z^2)^(-3/2)是單個(gè)圓線圈在軸線上產(chǎn)生的磁場表達(dá)式。具體到我的模型，電流i就為I1,I2.
版主所說的將I1,I2,作為x(4),x(5)這里也只是有3個(gè)方程啊，6個(gè)未知數(shù)，我加了約束條件使得a1=a2，就是5個(gè)未知數(shù)了。我再試試吧，我需要的不是具體的數(shù)值解，只要是能計(jì)算出R1,R2之間的關(guān)系，R1與a1之間的關(guān)系就行！

贊一下

回復(fù)此樓

18樓2012-05-11 11:00:33

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

521438

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 82.9
帖子: 19
在線: 13.2小時(shí)
蟲號(hào): 1801287
注冊: 2012-05-07
性別: MM
專業(yè): 電磁場與電路

http://edu.www.gaoyang168.com/attachment/04/4b/1801287_1336723264_727.jpg
這里我把我的數(shù)學(xué)模型也貼出來，如下：電流方向沒有標(biāo)注，2組線圈分別加反向的電流，即左側(cè)的線圈是順時(shí)針，右側(cè)的電流方向?yàn)槟鏁r(shí)針。

贊一下

回復(fù)此樓

19樓2012-05-11 16:05:14

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

Bruce_lee

銅蟲 (著名寫手)

應(yīng)助: 5 (幼兒園)
金幣: 2632.7
散金: 424
紅花: 4
帖子: 1277
在線: 317小時(shí)
蟲號(hào): 948103
注冊: 2010-01-25
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

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

★ ★
xiegangmai: 金幣+2, 謝謝參與！ 2013-03-25 23:02:38

function F = fsolveDemon1(x,I)
%FUN Summary of this function goes here
% Detailed explanation goes here
F=[I(1).*x(1).^2.*(3*x(3).*(x(1).^2+x(3).^2).^(-7/2)-7*x(3).^3*(x(1).^2+x(3).^2).^(-9/2))+I(2).*x(2).^2.*(3*x(3).*(x(2).^2+x(3).^2).^(-7/2)-7*x(3).^3.*(x(2).^2+x(3).^2).^(-9/2))
I(1).*x(1).^2*(-5*x(3).*(x(1).^2+x(3).^2).^(-9/2)+30*x(3).^3.*(x(1).^2+x(3).^2).^(-11/2)-33*x(3).^5.*(x(1).^2+x(3).^2).^(-13/2))+I(2).*x(2).^2.*(-5*x(3).*(x(2).^2+x(3).^2).^(-9/2)+30*x(3).^3.*(x(2).^2+x(3).^2).^(-11/2)-33*x(3).^5.*(x(2).^2+x(3).^2).^(-13/2))
I(1).*x(1).^2.*(35*x(3).*(x(1).^2+x(3).^2).^(-11/2)-385*x(3).^3.*(x(1).^2+x(3).^2).^(-13/2)+1001*x(3).^5.*(x(1).^2+x(3).^2).^(-15/2)-715*x(3).^7.*(x(1).^2+x(3).^2).^(-17/2))+I(2).*x(2).^2.*(35*x(3).*(x(2).^2+x(3).^2).^(-11/2)-385*x(3).^3.*(x(2).^2+x(3).^2).^(-13/2)+1001*x(3).^5.*(x(2).^2+x(3).^2).^(-15/2)-715*x(3).^7.*(x(2).^2+x(3).^2).^(-17/2))];
end
這是第一個(gè)函數(shù)。將其保存成fsolveDemon1.m文件。
function fsolveDemon2
x0=[0.1,0.1,0.1]; % x0中元素的個(gè)數(shù)與方程中x個(gè)數(shù)對(duì)應(yīng)相等，這里取都取0.1，也可自己修改。
options=optimset('MaxFunEvals',20000,'MaxIter',2000);
I=[1,2]; % 向量I有兩個(gè)值元素。這里權(quán)且取值1和2，自己可改成需要的數(shù)值。
X=fsolve(@(x)fsolveDemon1(x,I),x0,options);
for i=1:50 %迭代次數(shù)50次，可以自己修改。
[X,fval,exitflag]=fsolve(@(x)fsolveDemon1(x,I),X,options);
end
clc;
X=vpa(X,8) %顯示8為數(shù)字結(jié)果，如果想多顯示幾位或少顯示幾位，可以自己修改數(shù)值。
fval
exitflag
end
這是第二個(gè)函數(shù)。將其保存成fsolveDemon2.m文件。
使用時(shí)先將兩個(gè)文件拷貝到matlab能夠調(diào)用的指定文件夾中，然后在命令窗口輸入fsolveDemon2，回車即可。
結(jié)果中fval的值如果趨近于0，并且同時(shí)exitflag=1，這說明X就是所要求的解。

贊一下(2人)

回復(fù)此樓

努力就有收獲！

20樓2013-03-25 01:03:25

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

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

返回列表

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

[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +4	吃吃吃才有意義 2026-03-19	4/200	2026-03-19 13:48 by houyaoxu
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +3	Ncdx123456 2026-03-19	3/150	2026-03-19 13:18 by houyaoxu
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz