版塊導(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)劑公告

返回列表

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ù)此樓

» 猜你喜歡

中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有215人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(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é)理學(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)有1人回復(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的回帖

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的回帖

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的回帖

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的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 hustboshao 的主題更新

返回列表

普通表情龍兔虎貓

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

[考研] 275求調(diào)劑 +15	Micky11223 2026-03-25	20/1000	2026-03-29 20:44 by 唐沐兒
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[考研] 材料學(xué)碩333求調(diào)劑 +11	北道巷 2026-03-24	11/550	2026-03-29 08:31 by fmesaito
[考研] 0856求調(diào)劑 +13	zhn03 2026-03-25	14/700	2026-03-29 08:13 by fmesaito
[考研] 322求調(diào)劑 +7	宋明欣 2026-03-27	7/350	2026-03-28 21:27 by sanrepian
[考研] 復(fù)試調(diào)劑 +3	raojunqi0129 2026-03-28	3/150	2026-03-28 15:27 by 落�？伤�
[考研] 292求調(diào)劑 +14	鵝鵝鵝額額額額?/a> 2026-03-25	15/750	2026-03-28 08:45 by WYUMater
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +5	dongfang59 2026-03-27	5/250	2026-03-28 04:53 by wxiongid
[有機(jī)交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 一志愿西北大學(xué) 總分282 英語(yǔ)一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 求調(diào)劑 +3	劉柯@ 2026-03-24	4/200	2026-03-27 11:28 by shangxh
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓
[考研] 343求調(diào)劑 +4	贈(zèng)我一本書 2026-03-23	4/200	2026-03-27 00:40 by wxiongid
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5	星沉uu 2026-03-26	6/300	2026-03-26 19:02 by macy2011
[考研] 081700 調(diào)劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 一志愿河工大 081700 276求調(diào)劑 +4	地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 086003食品工程求調(diào)劑 +6	淼淼111 2026-03-24	6/300	2026-03-25 10:29 by 3Strings
[考研] 化工專碩求調(diào)劑 +3	question挽風(fēng) 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4	17501029541 2026-03-23	6/300	2026-03-24 10:21 by 戴圍脖的小蚊子