abingchem

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 2 (幼兒園)
金幣: 4871
帖子: 1672
在線: 117.6小時(shí)
蟲(chóng)號(hào): 206219

[交流] 【討論】怎么樣能高效些？

有個(gè)分段函數(shù)：
         | x + y, x>=0  and y >=0
         | x^2 + y, x>=0 and y<0
f(x,y) = | x + y^2,  x <0 and y >=0
         | x^2 + y^2, x<0 and y <0
如果參數(shù)x，y只是單個(gè)值，這樣用if語(yǔ)句就可以處理。

但如果x、y都是數(shù)組的話，怎么樣處理呢？x=-10:10;y=-10:10;

如果擴(kuò)展到x，y都是矩陣，有沒(méi)有什么好辦法？[x,y]=meshgrid([-10:0.1:10],[-10:0.1:10])

我能想到的最原始的辦法是
for  ....
   for ....
   ...
   end
end

回復(fù)此樓

» 猜你喜歡

26考研求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以已經(jīng)有5人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有7人回復(fù)
303求調(diào)劑已經(jīng)有3人回復(fù)
有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～已經(jīng)有7人回復(fù)
梁成偉老師課題組歡迎你的加入已經(jīng)有10人回復(fù)
一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) 已經(jīng)有8人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有3人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有4人回復(fù)

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

【課件】李開(kāi)復(fù)-Become A Good Speaker-教你如何掌握高效能演講技巧已經(jīng)有66人回復(fù)
【求助】石墨烯復(fù)合材料電容器性能評(píng)價(jià)【有效期至2011年06月17日】已經(jīng)有3人回復(fù)
【分享】如何有效的開(kāi)發(fā)自己的潛能【已搜索無(wú)重復(fù)】已經(jīng)有97人回復(fù)
【求助】如何計(jì)算能帶有效質(zhì)量？已經(jīng)有8人回復(fù)
【求助】如何選擇較小配體實(shí)現(xiàn)高效Suzuki偶聯(lián)反應(yīng)，并且盡量降低自偶聯(lián)可能性？已經(jīng)有9人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-01-04 16:16:46

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

xiegangmai

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

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
帖子: 4046
在線: 1330.2小時(shí)
蟲(chóng)號(hào): 653607

★ ★ ★
abingchem(金幣+2):經(jīng)典做法~ 2011-01-05 11:41:39
robert2020(金幣+3):辛苦了！ 2011-01-11 11:04:45

引用回帖:

Originally posted by abingchem at 2011-01-04 16:16:46:
有個(gè)分段函數(shù)：
         | x + y, x>=0  and y >=0
         | x^2 + y, x>=0 and y<0
f(x,y) = | x + y^2,  x <0 and y >=0
         | x^2 + y^2, x<0 and y <0
如果參 ...

不知道這樣對(duì)不對(duì)，樓主驗(yàn)證一下，和你用循環(huán)計(jì)算的結(jié)果比較比較。

CODE:

clear

clc

fxy = @( x, y )x .* ( x >= 0 & y >= 0 ) + y .* ( x >= 0 & y >= 0 ) +...

    x .^ 2 .* ( x >= 0 & y < 0 ) + y .* ( x >= 0 & y < 0 ) + ...

    x .* ( x < 0 & y >- 0 ) + y .^ 2 .* ( x < 0 & y >- 0 ) + ...

    x .^ 2 .* ( x < 0 & y < 0 ) + y .^ 2 .* ( x < 0 & y < 0 );

%% 如果x、y都是數(shù)組的話

x = -10 : 0.1 : 10;

y = -10 : 0.1 : 10;

z1 = fxy( x, y );

figure

plot3( x, y, z1 )

%% 擴(kuò)展到x，y都是矩陣

[ x, y ] = meshgrid(  -10 : 0.1 : 10, -10 : 0.1 : 10 );

z2 = fxy( x, y );

figure

surf( x, y, z2 )

贊一下(1人)

回復(fù)此樓

2樓2011-01-04 18:03:53

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

abingchem

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 2 (幼兒園)
金幣: 4871
帖子: 1672
在線: 117.6小時(shí)
蟲(chóng)號(hào): 206219

引用回帖:

Originally posted by xiegangmai at 2011-01-04 18:03:53:

fxy = @( x, y )x .* ( x >= 0 & y >= 0 ) + y .* ( x >= 0 & y >= 0 ) +...
x .^ 2 . ...

這個(gè)計(jì)算的表達(dá)式好像要進(jìn)行很多次的運(yùn)算，而且里面有大量是重復(fù)計(jì)算了，對(duì)吧？這樣效率會(huì)高嗎？如果[x,y]很大的話，時(shí)間比較長(zhǎng)。。。

順便問(wèn)一下，MATLAB7以后用哪個(gè)函數(shù)統(tǒng)計(jì)浮點(diǎn)運(yùn)算次數(shù)呢？

贊一下

回復(fù)此樓

3樓2011-01-05 11:44:09

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

lds588

鐵桿木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 15 (小學(xué)生)
金幣: 10949.6
帖子: 1444
在線: 206.3小時(shí)
蟲(chóng)號(hào): 890720

abingchem(金幣+1):多謝參與~ 2011-01-05 13:36:26

學(xué)習(xí)了~~~~~

回復(fù)此樓

4樓2011-01-05 12:50:41

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

xiegangmai

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

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
帖子: 4046
在線: 1330.2小時(shí)
蟲(chóng)號(hào): 653607

★ ★ ★
abingchem(金幣+3):多謝xiegangmai 2011-01-05 13:37:08
robert2020(金幣+3):辛苦了！ 2011-01-11 11:04:32

引用回帖:

Originally posted by abingchem at 2011-01-05 11:44:09:

這個(gè)計(jì)算的表達(dá)式好像要進(jìn)行很多次的運(yùn)算，而且里面有大量是重復(fù)計(jì)算了，對(duì)吧？這樣效率會(huì)高嗎？如果[x,y]很大的話，時(shí)間比較長(zhǎng)。。。

順便問(wèn)一下，MATLAB7以后用哪個(gè)函數(shù)統(tǒng)計(jì)浮點(diǎn)運(yùn)算次數(shù)呢？

你用tic、toc命令測(cè)試一下，相對(duì)循環(huán)方式，當(dāng)x，y越大，計(jì)算時(shí)間越省。
當(dāng)x、y是單個(gè)數(shù)或一維數(shù)組時(shí)，用循環(huán)方式可能速度還快一些。

MATLAB6以后的版本中就沒(méi)flops命令了，早期版本可以用flops統(tǒng)計(jì)浮點(diǎn)運(yùn)算次數(shù)，而更新版本中已經(jīng)不再提供該功能。

這是MATLAB幫助中的說(shuō)明，我的MATLAB 2009b。
Remove all existing instances of flops. With the incorporation of LAPACK in MATLAB version 6, counting floating-point operations is no longer practical.

贊一下(1人)

回復(fù)此樓

5樓2011-01-05 12:53:05

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

abingchem

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 2 (幼兒園)
金幣: 4871
帖子: 1672
在線: 117.6小時(shí)
蟲(chóng)號(hào): 206219

引用回帖:

Originally posted by xiegangmai at 2011-01-05 12:53:05:

你用tic、toc命令測(cè)試一下，相對(duì)循環(huán)方式，當(dāng)x，y越大，計(jì)算時(shí)間越省。
當(dāng)x、y是單個(gè)數(shù)或一維數(shù)組時(shí)，用循環(huán)方式可能速度還快一些。

MATLAB6以后的版本中就沒(méi)flops命令了，早期版本可以用flops統(tǒng)計(jì)浮點(diǎn)運(yùn) ...

恩，運(yùn)行結(jié)果就是這樣子的。數(shù)組大一些的時(shí)候效果更好！

贊一下

回復(fù)此樓

6樓2011-01-05 13:43:05

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

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

[考研] 材料與化工求調(diào)劑 +4	為學(xué)666 2026-03-16	4/200	2026-03-17 14:28 by houyaoxu
[碩博家園] 深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4	xujiaoszu 2026-03-11	9/450	2026-03-17 10:29 by xujiaoszu
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5	薛云鵬 2026-03-13	5/250	2026-03-17 10:15 by Sammy2
[考研] 304求調(diào)劑 +7	小熊joy 2026-03-14	7/350	2026-03-17 08:53 by 霧散后相遇lc
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-16	3/150	2026-03-16 10:07 by 求調(diào)劑zz
[考研] 326求調(diào)劑 +4	上岸的小葡 2026-03-15	5/250	2026-03-16 08:39 by Linda Hu
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[基金申請(qǐng)] 有必要更換申報(bào)口嗎 20+3	fannyamoy 2026-03-11	3/150	2026-03-14 00:52 by zhanghaozhu
[考研] 341求調(diào)劑 +4	番茄頭--- 2026-03-10	4/200	2026-03-13 23:12 by JourneyLucky
[考研] 285 求調(diào)劑資源與環(huán)境一志愿北京化工大學(xué) +3	未名考生 2026-03-10	3/150	2026-03-13 23:04 by JourneyLucky
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評(píng)估A類院�；瘜W(xué)工程與工藝，誠(chéng)求調(diào)劑 +7	小劉快快上岸 2026-03-11	7/350	2026-03-13 16:06 by ruiyingmiao
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 296求調(diào)劑 +3	大口吃飯身體健 2026-03-13	3/150	2026-03-13 10:31 by 學(xué)員8dgXkO