 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4688 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
深淵明宇新蟲 (小有名氣)
|
[求助]
求助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) |
回復(fù)此樓» 猜你喜歡
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有60人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有28人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有21人回復(fù)
w136851258
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 7565.1
- 散金: 328
- 紅花: 5
- 帖子: 385
- 在線: 135.7小時(shí)
- 蟲號(hào): 2611295
- 注冊(cè): 2013-08-19
- 專業(yè): 金屬材料
4樓2018-02-16 15:01:21
|
頂一下 發(fā)自小木蟲IOS客戶端 |
2樓2018-01-23 19:17:21
hechunlei
至尊木蟲 (知名作家)
胡蘿卜青椒
- 應(yīng)助: 2 (幼兒園)
- 金幣: 20658.5
- 散金: 114
- 紅花: 2
- 帖子: 5029
- 在線: 277.6小時(shí)
- 蟲號(hào): 2857621
- 注冊(cè): 2013-12-08
- 性別: GG
- 專業(yè): 零件加工制造
3樓2018-02-16 10:25:26
5樓2018-04-16 11:37:50
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 080200專業(yè)277分，求帶走！ +4 | 瓶子PZ 2026-03-31 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 08工科275分求調(diào)劑 +5 | AaAa7420 2026-03-31 | 5/250 |
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +13 | minht 2026-03-29 | 13/650 |
|
|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-30 | 11/550 |
|
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +19 | lllllos 2026-03-30 | 19/950 |
|
|
|
[考研] 263求調(diào)劑 +3 | DDDDuu 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿中海洋材料357 +4 | 麥恩莉. 2026-03-30 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +8 | 11ggg 2026-03-30 | 8/400 |
|
|
|
[考研] 354求調(diào)劑 +3 | lxb598 2026-03-31 | 4/200 |
|
|
|
[考研]
 一志愿大連理工大學(xué)材料求調(diào)劑
+6
|
Gymno 2026-03-30 | 6/300 |
|
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|