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

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

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

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

hustboshao

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1308.9
紅花: 1
帖子: 137
在線: 78小時(shí)
蟲號(hào): 1567118
注冊(cè): 2012-01-06
性別: GG
專業(yè): 固體力學(xué)

[求助] 很簡(jiǎn)單的一個(gè)matlab程序運(yùn)行時(shí)間過(guò)長(zhǎng)，得不到結(jié)果已有1人參與

主程序：
clc
clear all
global kesai omega0 m rou D a1 a3 K C %%主系統(tǒng)基本參數(shù)
global epsilong_T s_33 C_pzt C_pzt1 K_pzt K_pzt1 k0 d_33 s_0 l_0 theta L R %壓電材料參數(shù)
kesai=0.0013;m=0.44;omega0=62.83;rou=1.2;D=0.015;
a1=2.3;a3=-18;
K=1736.947916;
C=0.07187752;
s_0=0.005*0.005*pi; %cross section area of the PZT rod
l_0=0.01; %length of the PZT rod
epsilong_T=3400*8.854e-12; % permittivity under constant stress 常應(yīng)力下的介電常數(shù)
s_33=20.7e-12; %compliance of the piezoceramic under constant electric 恒定電場(chǎng)下的彈性柔順系數(shù)
d_33=593e-12; %壓電常數(shù)
C_pzt=epsilong_T*s_0/l_0; %capacitance between the electrodes of the PZT rod with no external force無(wú)外力下的電容
K_pzt=(1/s_33)*s_0/l_0; %stiffness of the short-circuited PZT rod 短路剛度
k0=d_33*(1/(sqrt(epsilong_T*s_33))); % electromechanical coupling cofficient 機(jī)電耦合系數(shù)
%C_pzt1=C_pzt*(1-k0^2);%capacitance of the PZT rod under constant strain 常應(yīng)變下的電容
C_pzt1=1.2e-7;
%K_pzt1=K_pzt*(1/(1-k0^2)); % stiffness of the PZT rod with open electrodes 回路剛度
%K_pzt1=C_pzt1/(epsilong_T*s_33*(1-k0^2)^2);
K_pzt1=K*0.25;
%theta=k0/(1-k0^2)*(1/(sqrt(epsilong_T*s_33)));% electromechanical coupling factor
theta=1.55e-3;
L=0.000001; % 電感
R=10; % 電阻
step_size=0.01; %時(shí)間步長(zhǎng)
precision=1e-9; %精度量級(jí)1
options=odeset('reltol',precision);
N=20;
NN=N/0.01;
tic
U=5
%[T1,Y1]=ode45(@uncontrol,[0:step_size:N],[0.001 0],options,U);
[T2,Y2]=ode45(@control,[0:step_size:N],[0.001 0 0.001 0],options,U);
toc
%w_1=max(Y1(NN/2:end,1))/0.015
w_2=max(Y2(NN/2:end,1))
w_2_Q=max(Y2(NN/2:end,3))
%figure(1)
%plot(T1,Y1(:,1)/0.015,'r')
%hold on
%figure(2)
%plot(T2,Y2(:,1)/0.015,'b')
%hold on
%figure(3)
%plot(T2,Y2(:,3)/0.015,'G')
%hold on
figure(4)
plot(T2,Y2(:,1)/0.015,'-r',T2,Y2(:,3)/0.015,'--b')
xlabel('Time','FontName','Times New Roman','FontSize',20)
ylabel('w_1,w_2','FontName','Times New Roman','FontSize',20)
legend('uncontrol','control')
hold on

子程序：
function dy=control(t,y,U)
global kesai omega0 m rou D a1 a3 K C %%主系統(tǒng)基本參數(shù)
global epsilong_T s_33 C_pzt C_pzt1 K_pzt K_pzt1 k0 d_33 s_0 l_0 theta L R %壓電材料參數(shù)
dy=[y(2);
(rou*D*U*a1*y(2)*0.5+rou*D*a3*y(2)^3./U*0.5-C*y(2)-K*y(1))./m-K_pzt1*y(1)./m+theta*y(3)./m;
y(4);
-R*y(4)./L-y(3)./(L*C_pzt1)+theta*y(1)./L];
end

回復(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ù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有89人回復(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)有26人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對(duì)論效應(yīng)的影響已經(jīng)有17人回復(fù)
大家看看，我該如何投稿到哪個(gè)期刊？《從金庸的武學(xué)視角談?dòng)?jì)算機(jī)編程思維的培養(yǎng)》已經(jīng)有10人回復(fù)

I am what i am

1樓 2015-11-09 15:56:56

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

hustboshao

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1308.9
紅花: 1
帖子: 137
在線: 78小時(shí)
蟲號(hào): 1567118
注冊(cè): 2012-01-06
性別: GG
專業(yè): 固體力學(xué)

引用回帖:

2樓: Originally posted by getengqing at 2015-11-09 16:49:24
這么一堆，誰(shuí)會(huì)仔細(xì)看啊，建議樓主問(wèn)詳細(xì)一點(diǎn)的問(wèn)題···

不好意思，參數(shù)寫的有點(diǎn)多，主程序其實(shí)就是一個(gè)ODE45函數(shù)，很簡(jiǎn)單的，但是不知道為什么算不出來(lái)

贊一下

回復(fù)此樓

I am what i am

3樓2015-11-09 16:59:43

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

查看全部 5 個(gè)回答

getengqing

木蟲 (正式寫手)

應(yīng)助: 104 (高中生)
金幣: 3759.8
散金: 30
紅花: 5
帖子: 592
在線: 110.5小時(shí)
蟲號(hào): 2986702
注冊(cè): 2014-02-23
性別: GG
專業(yè): 自然語(yǔ)言理解與機(jī)器翻譯

這么一堆，誰(shuí)會(huì)仔細(xì)看啊，建議樓主問(wèn)詳細(xì)一點(diǎn)的問(wèn)題···

贊一下

回復(fù)此樓

一起交流學(xué)習(xí)/分享優(yōu)秀資源

2樓2015-11-09 16:49:24

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

lsjc77

銅蟲 (初入文壇)

應(yīng)助: 8 (幼兒園)
金幣: 263.3
紅花: 6
帖子: 34
在線: 124.6小時(shí)
蟲號(hào): 1002599
注冊(cè): 2010-04-21
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
hustboshao: 金幣+10, ★★★很有幫助 2015-11-12 10:15:42

換成ode23s或者15s，沒(méi)仔細(xì)看你的程序，但是發(fā)現(xiàn)很多很小的數(shù)值，感覺(jué)你的問(wèn)題應(yīng)該是剛性的，剛性方程不能用ode45解。
還有一個(gè)建議就是無(wú)量綱化，讓你要解的參數(shù)在 0 1 之間，或者至少在一個(gè)數(shù)量級(jí)上。
不過(guò)你要解的方程很少，如果只是用ode45時(shí)間長(zhǎng)，那換成ode23s/15s應(yīng)該會(huì)很快給你一個(gè)結(jié)果。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

hustboshao

4樓2015-11-10 14:54:40

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

hustboshao

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1308.9
紅花: 1
帖子: 137
在線: 78小時(shí)
蟲號(hào): 1567118
注冊(cè): 2012-01-06
性別: GG
專業(yè): 固體力學(xué)

送紅花一朵

引用回帖:

4樓: Originally posted by lsjc77 at 2015-11-10 14:54:40
換成ode23s或者15s，沒(méi)仔細(xì)看你的程序，但是發(fā)現(xiàn)很多很小的數(shù)值，感覺(jué)你的問(wèn)題應(yīng)該是剛性的，剛性方程不能用ode45解。
還有一個(gè)建議就是無(wú)量綱化，讓你要解的參數(shù)在 0 1 之間，或者至少在一個(gè)數(shù)量級(jí)上。
不過(guò)你要 ...

無(wú)量綱化是個(gè)好辦法，謝啦

贊一下

回復(fù)此樓

I am what i am

5樓2015-11-12 10:16:14

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

查看全部 5 個(gè)回答

普通表情龍兔虎貓

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

[考研] 求調(diào)劑推薦材料 304 +21	荷包蛋hyj 2026-03-26	21/1050	2026-04-01 21:09 by lijunpoly
[考研] 085601材料工程找調(diào)劑 +20	oatmealR 2026-03-29	21/1050	2026-04-01 21:00 by lijunpoly
[考研] 一志愿西安交大材料學(xué)碩（英一數(shù)二）347，求調(diào)劑到高分子/材料相關(guān)專業(yè) +7	zju51 2026-03-31	9/450	2026-04-01 19:35 by CFQZAFU
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +10	腳滑的守法公民 2026-03-27	10/500	2026-04-01 17:23 by pies112
[考研] 調(diào)劑推薦 +9	清酒714 2026-03-26	10/500	2026-04-01 16:25 by yanflower7133
[考研] 0710生物學(xué)考研調(diào)劑 +3	李多米lee. 2026-03-27	4/200	2026-04-01 16:21 by zzchen2000
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +11	傾____蕭 2026-03-27	12/600	2026-04-01 15:27 by oooqiao
[考研] 環(huán)境工程調(diào)劑 +9	hyzzzzzzz. 2026-04-01	9/450	2026-04-01 14:20 by salamander`
[考研] 一志愿同濟(jì)大學(xué)323分（080500）求調(diào)劑 +4	yikeniu 2026-04-01	4/200	2026-04-01 14:06 by asdfzly
[考研] 材料工程085601數(shù)二英一335求調(diào)劑 +5	雙馬尾痞老板2 2026-03-31	5/250	2026-03-31 19:07 by Wang200018
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +9	真名有冰 2026-03-29	10/500	2026-03-31 18:01 by xhai2011
[考研] 求調(diào)劑 +9	akdhjs 2026-03-31	10/500	2026-03-31 17:14 by zhyzzh
[考研] 083000學(xué)碩274求調(diào)劑 +12	Li李魚 2026-03-26	12/600	2026-03-31 10:01 by cal0306
[考研] 292求調(diào)劑 +13	是妍子也是研子 2026-03-30	13/650	2026-03-30 18:01 by 小徐0109
[考研] 材料化工340求調(diào)劑 +3	jhx777 2026-03-30	3/150	2026-03-30 17:54 by JourneyLucky
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7	Shern—- 2026-03-28	7/350	2026-03-30 16:31 by nothing投稿中
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +6	niko- 2026-03-27	6/300	2026-03-30 10:25 by herarysara
[考研] 298求調(diào)劑 +4	種圣賜 2026-03-28	4/200	2026-03-29 08:42 by q1092522407
[考研] 求調(diào)劑 +7	爭(zhēng)取九點(diǎn)睡 2026-03-28	8/400	2026-03-28 21:07 by 爭(zhēng)取九點(diǎn)睡
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002