| 24小時(shí)熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 4687 | 回復(fù): 4 | ||
深淵明宇新蟲(chóng) (小有名氣)
|
[求助]
求助matlab程序計(jì)算速度過(guò)慢的原因
|
|
程序代碼如下 function [length]=contactlength(x0) if x0>50||x0<-50 'error:數(shù)據(jù)超出尺寸范圍' else syms x R=300;%非球面頂點(diǎn)曲率半徑 c=1/R; delta=0.1; k=-3.3;%非球面參數(shù) rb=27;%半徑 y=(-1*c*x.^2)./(1+sqrt(1-(1+k)*(c^2)*x.^2)); dy=diff(y); dy2=diff(y,2); dyx0=subs(dy,x0); dy2x0=subs(dy2,x0); r0=(1+dyx0^2)^(3/2)/(abs(dy2x0))%曲率半徑 y0=(-1*c*x0^2)./(1+sqrt(1-(1+k)*(c^2)*x0^2));%非球面接觸點(diǎn)y坐標(biāo) if x0==0 xb=x0; yb=y0+rb-delta; ob=[xb,yb]%圓心二維坐標(biāo) else xb=x0+(rb-delta)/sqrt(1+1/dyx0^2); yb=y0+(-1/dyx0)*(rb-delta)/sqrt(1+1/dyx0^2); ob=[xb,yb]%圓心二維坐標(biāo) end syms x2 y2 ezplot((x2-xb)^2+(y2-yb)^2==rb^2,[-50 50 -10 70]) axis equal; hold on; ezplot(y2==(-1*c*x2^2)/(1+sqrt(1-(1+k)*(c^2)*x2^2)),[-50 50 -10 70]) s=solve(y2==(-1*c*x2^2)./(1+sqrt(1-(1+k)*(c^2)*x2^2)),... (x2-xb)^2+(y2-yb)^2==rb^2,x2,y2); X=double(s.x2); Y=double(s.y2); X=X(~logical(imag(X))); length=max(X)-min(X) end 在實(shí)際運(yùn)算過(guò)程中,要100多秒才能算出結(jié)果,請(qǐng)問(wèn)是哪里有問(wèn)題呢?我該怎么改程序? 本人新手,接觸matlab不久,希望大神們多多指導(dǎo) |
至尊木蟲(chóng) (知名作家)
胡蘿卜青椒

鐵桿木蟲(chóng) (正式寫手)
|
可以考慮使用tic toc,查看具體哪一段代碼耗時(shí)過(guò)長(zhǎng),針對(duì)性優(yōu)化 發(fā)自小木蟲(chóng)Android客戶端 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856求調(diào)劑 +9 | 楒桉 2026-03-28 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +9 | 小熊~~ 2026-03-31 | 9/450 |
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 13/650 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-30 | 6/300 |
|
|
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135,現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|
[考研] 吉大生物學(xué)326分求調(diào)劑 +3 | sunnyupup 2026-03-31 | 3/150 |
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 3/150 |
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
[考研] 292求調(diào)劑 +13 | 是妍子也是研子 2026-03-30 | 13/650 |
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué),化工,生物,生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入! +3 | llss0711 2026-03-28 | 6/300 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲(chóng)085600 2026-03-27 | 8/400 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-25 | 7/350 |
|