brqhl_ing

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線: 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專業(yè): 物理電子學(xué)

[交流] 【求助】即FOR循環(huán)之后的問(wèn)題已有3人參與

function dx = Rossler(t,x，k1,k2)
dx1=-(1+0.015)*x(2)-x(3)+k1*(x(4)-x(1))+k2*(-x(4)-x(1))*(x(4)-x(1));
dx2=(1+0.015)*x(1)+0.15*x(2);
dx3=0.5+x(3)*(x(1)-10);
dx4=-(1+0.015)*x(5)-x(6)-k1*(x(4)-x(1))-k2*(-x(4)-x(1))*(x(4)-x(1));
dx5=(1+0.015)*x(4)+0.15*x(5);
dx6=0.5+x(6)*(x(4)-10);
dx=[dx1;dx2;dx3;dx4;dx5;dx6];
這是上個(gè)FOR循環(huán)里面的方程只不過(guò)這次控制項(xiàng)不同有k1 k2
同樣定義誤差函數(shù)e=((x(:,4)-x(:,1)).^2+(x(:,5)-x(:,2)).^2+(x(:,6)-x(:,3)).^2).^0.5;
如何找出k1 k2 取值使e=0 這個(gè)e可以是耦合系統(tǒng)經(jīng)過(guò)一定時(shí)間達(dá)到誤差函數(shù)e=0畫(huà)出k1 k2使e=0的區(qū)域（k1 k2相圖）

回復(fù)此樓

» 猜你喜歡

食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
311（085601）求調(diào)劑已經(jīng)有4人回復(fù)
305求調(diào)劑已經(jīng)有4人回復(fù)
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑已經(jīng)有7人回復(fù)
一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑已經(jīng)有12人回復(fù)
289求調(diào)劑已經(jīng)有11人回復(fù)
生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 已經(jīng)有6人回復(fù)
321求調(diào)劑已經(jīng)有7人回復(fù)
343求調(diào)劑已經(jīng)有4人回復(fù)

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

FORTRAN新手求助主程序循環(huán)問(wèn)題已經(jīng)有10人回復(fù)
【交流】●物理學(xué)打假行動(dòng)，請(qǐng)起訴愛(ài)因斯坦● 已經(jīng)有58人回復(fù)
【整理】今天整理硬盤(pán)上的資料，偶然發(fā)現(xiàn)的好玩的小東西就直接發(fā)上來(lái)吧已經(jīng)有105人回復(fù)

1樓 2010-07-13 20:59:06

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

xiegangmai

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

我沒(méi)頭銜

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

★
anyuezhiji(金幣+1):謝謝提供意見(jiàn) 2010-07-14 07:38:28

這個(gè)問(wèn)題好像不能用優(yōu)化來(lái)做，不能寫(xiě)出k1、k2的目標(biāo)函數(shù)。

有個(gè)笨一些的方法，假設(shè)知道k1，k2的范圍，先取k1，k2較大的步長(zhǎng)，雙重循環(huán)計(jì)算e。找出e接近0的時(shí)候k1，k2的較小范圍，再將k1，k2的步長(zhǎng)縮小，找到使e=0的k1，k2。

先試試吧。
看看高手有沒(méi)有別的辦法。

贊一下(1人)

回復(fù)此樓

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

2樓2010-07-13 23:51:05

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

brqhl_ing

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

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線: 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專業(yè): 物理電子學(xué)

引用回帖:

Originally posted by xiegangmai at 2010-07-13 23:51:05:
這個(gè)問(wèn)題好像不能用優(yōu)化來(lái)做，不能寫(xiě)出k1、k2的目標(biāo)函數(shù)。

有個(gè)笨一些的方法，假設(shè)知道k1，k2的范圍，先取k1，k2較大的步長(zhǎng)，雙重循環(huán)計(jì)算e。找出e接近0的時(shí)候k1，k2的較小范圍，再將k1，k2的步長(zhǎng)縮小，找到使 ...

這個(gè)k1 k2的取值范圍是知道的一半都是在0到1之間。就是去k1 k2一個(gè)計(jì)算步長(zhǎng) 相當(dāng)于把空間分成很多點(diǎn)組成的計(jì)算這些點(diǎn)的e的情況

贊一下

回復(fù)此樓

3樓2010-07-14 08:32:33

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

brqhl_ing

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

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線: 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專業(yè): 物理電子學(xué)

robert2020:匆匆又到暑期畢業(yè)季，大家有點(diǎn)忙...煩請(qǐng)耐心等待... 2010-07-16 16:21:08

大家快來(lái)討論討論如何解決這個(gè)問(wèn)題啊

贊一下(1人)

回復(fù)此樓

4樓2010-07-14 21:02:09

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

brqhl_ing

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

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線: 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專業(yè): 物理電子學(xué)

大家快來(lái)討論討論如何解決這個(gè)問(wèn)題啊!大俠們快來(lái)�。�！

贊一下

回復(fù)此樓

5樓2010-07-16 10:49:58

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

brqhl_ing

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

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線: 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專業(yè): 物理電子學(xué)

不能沉了�。。「呤诌€沒(méi)出現(xiàn)呢

贊一下

回復(fù)此樓

6樓2010-07-17 08:56:54

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

free--fire

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

應(yīng)助: 3 (幼兒園)
金幣: 543.5
散金: 570
帖子: 199
在線: 51.7小時(shí)
蟲(chóng)號(hào): 1051983
注冊(cè): 2010-07-05
性別: GG
專業(yè): 力學(xué)中的基本問(wèn)題和方法

★ ★
brqhl_ing(金幣+10): 2010-07-17 11:40:27
anyuezhiji(金幣+2):感謝應(yīng)助，歡迎常來(lái)！ 2010-07-17 17:49:56
nono2009:代碼建議用[code][/code]，就不會(huì)有顯示問(wèn)題。 2010-07-18 11:11:38

引用回帖:

Originally posted by brqhl_ing at 2010-07-13 20:59:06:
function dx = Rossler(t,x，k1,k2)
dx1=-(1+0.015)*x(2)-x(3)+k1*(x(4)-x(1))+k2*(-x(4)-x(1))*(x(4)-x(1));
dx2=(1+0.015)*x(1)+0.15*x(2);
dx3=0.5+x(3)*(x(1)-10);
dx4=-(1+0.015)*x(5)-x(6)-k1*(x(4)-x( ...

CODE:

function help100716()

clc;

clear all;

tic;

global k1 k2 ;

n=0:0.01:0.5;

s=length(n);

K=zeros(s^2,3);

m=1;

for k1=0:0.01:0.5

    for k2=0:0.01:0.5

        [T,x]=ode45(@Rossler,[0,600],[0.2;0.3;0.6;0.4;0.7;0.1]);

        e=((x(:,4)-x(:,1)).^2+(x(:,5)-x(:,2)).^2+(x(:,6)-x(:,3)).^2).^0.5;

        if(min(e)<=1e-1)%% set the bar of the 'e' value

            K(m,:)=[k1 k2 min(e)];

            m=m+1;

        end

    end

end

if m>1

figure;

plot(K(1:m-1,1:2));

end

toc;

function dx = Rossler(t,x)

global k1 k2;

dx1=-(1+0.015)*x(2)-x(3)+k1*(x(4)-x(1))+k2*(-x(4)-x(1))*(x(4)-x(1));

dx2=(1+0.015)*x(1)+0.15*x(2);

dx3=0.5+x(3)*(x(1)-10);

dx4=-(1+0.015)*x(5)-x(6)-k1*(x(4)-x(1))-k2*(-x(4)-x(1))*(x(4)-x(1));

dx5=(1+0.015)*x(4)+0.15*x(5);

dx6=0.5+x(6)*(x(4)-10);

dx=[dx1;dx2;dx3;dx4;dx5;dx6];

[ Last edited by nono2009 on 2010-7-18 at 11:10 ]

贊一下(2人)

回復(fù)此樓

我是流，是浮云天上游；我是虎，是大蟲(chóng)臥山谷。

7樓2010-07-17 10:57:57

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

free--fire

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

應(yīng)助: 3 (幼兒園)
金幣: 543.5
散金: 570
帖子: 199
在線: 51.7小時(shí)
蟲(chóng)號(hào): 1051983
注冊(cè): 2010-07-05
性別: GG
專業(yè): 力學(xué)中的基本問(wèn)題和方法

★
adu886886(金幣+1):謝謝提供意見(jiàn) 2010-07-18 22:13:28

應(yīng)該是 K(m,: )=[k1 k2 min(e)];
另外就是計(jì)算量太大了，最好是根據(jù)k1和k2的對(duì)稱性將計(jì)算量縮減一半。
如果可以的話請(qǐng)將結(jié)果告訴大家。
good luck !

贊一下(1人)

回復(fù)此樓

我是流，是浮云天上游；我是虎，是大蟲(chóng)臥山谷。

8樓2010-07-17 11:01:44

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

brqhl_ing

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

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線: 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專業(yè): 物理電子學(xué)

引用回帖:

Originally posted by free--fire at 2010-07-17 10:57:57:

function help100716()

clc;
clear all;
tic;

global k1 k2 ;
n=0:0.01:0.5;
s=length(n);
K=zeros(s^2,3);
m=1;
for k1=0:0.01:0.5
for k2=0:0.01:0.5
[T,x]=ode45(@Rossler,[0,60 ...

非常之感謝我在調(diào)試中！！

贊一下

回復(fù)此樓

9樓2010-07-17 14:57:26

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

brqhl_ing

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

應(yīng)助: 0 (幼兒園)
金幣: 262.9
散金: 11
紅花: 1
帖子: 115
在線: 19.5小時(shí)
蟲(chóng)號(hào): 854693
注冊(cè): 2009-09-23
性別: GG
專業(yè): 物理電子學(xué)

引用回帖:

這個(gè)是我把 k1 k2 n 的值改為0.04 時(shí)出的圖感覺(jué)有點(diǎn)不對(duì)！取值太大計(jì)算時(shí)間太長(zhǎng)。所以我改了下就出現(xiàn)這樣的圖了

贊一下

回復(fù)此樓

10樓2010-07-18 09:15:20

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

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

[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落�？伤�
[考研] 求調(diào)劑 +6	蘆lty 2026-03-25	7/350	2026-03-28 13:13 by 唐沐兒
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3	@taotao 2026-03-28	3/150	2026-03-28 10:26 by JourneyLucky
[考研] 328求調(diào)劑 +7	嗯滴的基本都 2026-03-27	7/350	2026-03-28 04:19 by fmesaito
[考研] 307求調(diào)劑 +8	超級(jí)伊昂大王 2026-03-24	9/450	2026-03-27 15:34 by 超級(jí)伊昂大王
[考研] 0856調(diào)劑 +5	求求讓我有書(shū)讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 312求調(diào)劑 +9	上岸吧ZJY 2026-03-22	13/650	2026-03-27 11:24 by sanrepian
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚(yú)的貓
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4	@taotao 2026-03-26	5/250	2026-03-27 08:10 by hypershenger
[考研] 317求調(diào)劑 +7	蛋黃咸肉粽 2026-03-26	7/350	2026-03-27 02:29 by fmesaito
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4	為學(xué)666 2026-03-24	4/200	2026-03-25 11:27 by BruceLiu320
[考研] 318求調(diào)劑 +3	plum李子 2026-03-23	3/150	2026-03-25 09:42 by 霧散后相遇lc
[考研] 材料專碩找調(diào)劑 +5	哈哈哈吼吼吼哈 2026-03-23	5/250	2026-03-24 19:07 by 了了了了。。
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	6/300	2026-03-23 21:45 by 子系博
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4	allen-yin 2026-03-23	6/300	2026-03-23 15:04 by 汪��？！
[考研] 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc