lkcool7

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 87.5
散金: 3
帖子: 22
在線: 6.1小時(shí)
蟲號(hào): 3859120
注冊(cè): 2015-05-09
專業(yè): 自動(dòng)化

[求助] 運(yùn)行matlab程序中自定義函數(shù)時(shí)出問題了，還請(qǐng)大神幫忙！謝謝已有1人參與

我自定義了三個(gè)函數(shù)Y1()，Y2()，Y3()，三個(gè)函數(shù)只是內(nèi)部參數(shù)k發(fā)生變
在這里我給出Y1()代碼

function h=Y1(x)

%參數(shù)
k=0.1884;l=25.5;pA=1070;
u=cosh(k*x)+cos(k*x)-((sinh(k*x)+sin(k*x))*(cosh(k*l)-cos(k*l)))/(sinh(k*l)-sin(k*l));

a=sqrt(pA*eval(int(u^2,x,0,l)));
h=u/a;

下面是主程序，目的在于繪出Y1，Y2，Y3曲線圖。
clear
clc
x=0.1:0.1:25.5;

for m=1:1:255
y1(m)=Y1(x(m));
end
for m=1:1:255
y2(m)=Y2(x(m));
end
for m=1:1:255
y3(m)=Y3(x(m));
end
plot(x,y1);
hold on
plot(x,y2);
hold on
plot(x,y3);
hold on

運(yùn)行是程序在調(diào)用Y1()時(shí)出錯(cuò)，提示為

??? Undefined function or method 'int' for input arguments of type 'double'.

Error in ==> Y1 at 4
a=sqrt(pA*eval(int(u^2,x,0,l)));

Error in ==> zhenxing at 6
y1(m)=Y1(x(m));

也就是在Y1（）中調(diào)用int積分函數(shù)時(shí)出錯(cuò)，為什么會(huì)這樣�。�？？？？如果我將a=sqrt(pA*eval(int(u^2,x,0,l)));提前算出，直接將這行替換為a=169.2216;卻可以運(yùn)行，并最后繪制出我要的曲線。但是在我自己定義的函數(shù)中調(diào)用int()居然會(huì)出錯(cuò)？求大神幫助！謝謝

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有94人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

運(yùn)行matlab程序出現(xiàn)警告：TIFF庫警告，怎么解決啊已經(jīng)有0人回復(fù)
matlab程序運(yùn)行報(bào)錯(cuò) 求解已經(jīng)有2人回復(fù)
已有matlab程序如何導(dǎo)入excel數(shù)據(jù)運(yùn)行已經(jīng)有9人回復(fù)
已有matlab程序如何導(dǎo)入excel數(shù)據(jù)運(yùn)行已經(jīng)有2人回復(fù)
求助！哪位大神幫忙運(yùn)行下MATLAB程序，求結(jié)果！已經(jīng)有3人回復(fù)
matlab程序運(yùn)行報(bào)錯(cuò)求助已經(jīng)有7人回復(fù)
求助MATLAB程序已經(jīng)有5人回復(fù)
Matlab中如何運(yùn)行C++程序已經(jīng)有3人回復(fù)
matlab程序已經(jīng)有7人回復(fù)
Matlab程序運(yùn)行結(jié)束提示已經(jīng)有5人回復(fù)
MATLAB程序調(diào)試已經(jīng)有9人回復(fù)
MATLAB程序調(diào)整已經(jīng)有4人回復(fù)
matlab程序運(yùn)行時(shí)快時(shí)慢已經(jīng)有8人回復(fù)
一個(gè)MATLAB程序，運(yùn)行不出來已經(jīng)有1人回復(fù)
一個(gè)matlab程序，運(yùn)行不出來已經(jīng)有2人回復(fù)
matlab程序運(yùn)行出錯(cuò)，求指點(diǎn)修改已經(jīng)有4人回復(fù)
matlab程序只運(yùn)行一次就跳出的問題已經(jīng)有2人回復(fù)
matlab程序運(yùn)行中出現(xiàn)的函數(shù)問題已經(jīng)有1人回復(fù)

1樓 2015-11-12 22:42:06

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

轉(zhuǎn)基因奔奔

金蟲 (初入文壇)

應(yīng)助: 12 (小學(xué)生)
金幣: 827.2
紅花: 2
帖子: 49
在線: 30.8小時(shí)
蟲號(hào): 2956857
注冊(cè): 2014-02-06
專業(yè): 自動(dòng)化

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

感謝參與，應(yīng)助指數(shù) +1

你的u是一個(gè)數(shù)，不是一個(gè)表達(dá)式，對(duì)一個(gè)數(shù)怎么積分嘞？

贊一下(2人)

回復(fù)此樓

2樓2015-11-13 09:31:55

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

lkcool7

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 87.5
散金: 3
帖子: 22
在線: 6.1小時(shí)
蟲號(hào): 3859120
注冊(cè): 2015-05-09
專業(yè): 自動(dòng)化

引用回帖:

2樓: Originally posted by 轉(zhuǎn)基因奔奔 at 2015-11-13 09:31:55
你的u是一個(gè)數(shù)，不是一個(gè)表達(dá)式，對(duì)一個(gè)數(shù)怎么積分嘞？

u=cosh(k*x)+cos(k*x)-((sinh(k*x)+sin(k*x))*(cosh(k*l)-cos(k*l)))/(sinh(k*l)-sin(k*l));
這個(gè)是u的定義，是以x為自變量的函數(shù)啊

贊一下

回復(fù)此樓

3樓2015-11-13 12:15:17

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

轉(zhuǎn)基因奔奔

金蟲 (初入文壇)

應(yīng)助: 12 (小學(xué)生)
金幣: 827.2
紅花: 2
帖子: 49
在線: 30.8小時(shí)
蟲號(hào): 2956857
注冊(cè): 2014-02-06
專業(yè): 自動(dòng)化

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

★ ★ ★ ★ ★
lkcool7: 金幣+5, ★★★很有幫助 2015-11-16 09:57:06

你的x是以輸入變量的形式帶入的，每次循環(huán)都是以一個(gè)數(shù)作為x輸入的。MATLAB會(huì)自動(dòng)解算出u的數(shù)值，而不是把u當(dāng)做一個(gè)表達(dá)式

贊一下

回復(fù)此樓

4樓2015-11-13 16:52:07

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

小木蟲zb

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時(shí)
蟲號(hào): 1227944
注冊(cè): 2011-03-10
性別: GG
專業(yè): 分離過程

樓主最好把公式和你想解決的問題完整的放上來

贊一下

回復(fù)此樓

5樓2015-11-13 20:59:56

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

lkcool7

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 87.5
散金: 3
帖子: 22
在線: 6.1小時(shí)
蟲號(hào): 3859120
注冊(cè): 2015-05-09
專業(yè): 自動(dòng)化

引用回帖:

4樓: Originally posted by 轉(zhuǎn)基因奔奔 at 2015-11-13 16:52:07
你的x是以輸入變量的形式帶入的，每次循環(huán)都是以一個(gè)數(shù)作為x輸入的。MATLAB會(huì)自動(dòng)解算出u的數(shù)值，而不是把u當(dāng)做一個(gè)表達(dá)式

謝謝回答，我已解決了問題，確實(shí)發(fā)覺函數(shù)Y1(x)中的u為數(shù)值而非符號(hào)表達(dá)式。然后我建立了符號(hào)變量 syms x k l 將int(u,x,0,l) 計(jì)算出來
下面是我解決問題后的代碼
function h=Y1(x)
k=0.1884;l=25.5;pA=1070;

u=cosh(k*x)+cos(k*x)-((sinh(k*x)+sin(k*x))*(cosh(k*l)-cos(k*l)))/(sinh(k*l)-sin(k*l));

%p為符號(hào)積分int(u^2,x,0,l)結(jié)果，為符號(hào)表達(dá)式
p=1/4*-4*cosh(k*l)*cos(k*l)^2*exp(3*k*l)-4*sin(k*l)^2*cosh(k*l)*exp(3*k*l)+4*sinh(k*l)*cos(k*l)^3*exp(2*k*l)-sin(k*l)*sinh(k*l)-cosh(k*l)*cos(k*l)+cosh(k*l)^2-4*cos(k*l)^3*exp(k*l)+8*sin(k*l)*sinh(k*l)*k*exp(2*k*l)*l+4*exp(2*k*l)*sinh(k*l)*sin(k*l)^2*cos(k*l)+4*sin(k*l)^2*cos(k*l)*exp(3*k*l)-4*sin(k*l)^2*cos(k*l)*exp(k*l)+4*cos(k*l)^3*exp(3*k*l)-exp(4*k*l)*cosh(k*l)^2+4*sinh(k*l)*cosh(k*l)*exp(3*k*l)*sin(k*l)+4*cosh(k*l)*cos(k*l)*exp(3*k*l)*sin(k*l)-4*sinh(k*l)*cosh(k*l)*sin(k*l)*exp(k*l)-4*exp(2*k*l)*cosh(k*l)^2*k*l-4*sinh(k*l)*cosh(k*l)*cos(k*l)^2*exp(2*k*l)+4*cosh(k*l)*cos(k*l)*sin(k*l)*exp(k*l)+4*exp(2*k*l)*cos(k*l)^2*k*l+exp(4*k*l)*sinh(k*l)*cosh(k*l)+exp(4*k*l)*sinh(k*l)*sin(k*l)+exp(4*k*l)*sin(k*l)*cos(k*l)-exp(4*k*l)*sin(k*l)*cosh(k*l)+4*sin(k*l)^2*cosh(k*l)*exp(k*l)-4*cosh(k*l)^2*sin(k*l)*exp(k*l)+exp(4*k*l)*cosh(k*l)*cos(k*l)-exp(4*k*l)*sinh(k*l)*cos(k*l)-4*cosh(k*l)^2*exp(3*k*l)*sin(k*l)+4*cosh(k*l)*cos(k*l)^2*exp(k*l)+4*sin(k*l)^2*sinh(k*l)*exp(k*l)+4*sin(k*l)^2*sinh(k*l)*exp(3*k*l)+sinh(k*l)*cosh(k*l)-sinh(k*l)*cos(k*l)-sin(k*l)*cosh(k*l)+sin(k*l)*cos(k*l)+2*exp(2*k*l)*sin(k*l)*cos(k*l)-2*exp(2*k*l)*sin(k*l)*cosh(k*l)+2*exp(2*k*l)*sinh(k*l)*cosh(k*l)-2*exp(2*k*l)*sinh(k*l)*cos(k*l))/k/(-cosh(k*l)^2+cos(k*l)^2+2*sin(k*l)*sinh(k*l))*exp(-2*k*l);

a=sqrt(pA*p);
h=u/a;

贊一下

回復(fù)此樓

6樓2015-11-16 10:05:19

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

yanshuiting

新蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1484.8
散金: 50
紅花: 7
帖子: 448
在線: 68.7小時(shí)
蟲號(hào): 2633300
注冊(cè): 2013-09-03
專業(yè): 計(jì)算機(jī)科學(xué)

引用回帖:

3樓: Originally posted by lkcool7 at 2015-11-13 12:15:17
u=cosh(k*x)+cos(k*x)-((sinh(k*x)+sin(k*x))*(cosh(k*l)-cos(k*l)))/(sinh(k*l)-sin(k*l));
這個(gè)是u的定義，是以x為自變量的函數(shù)啊...

樓上說的是對(duì)的，x和u是符號(hào)變量，才能使用 int()函數(shù)，你的代碼里面一是沒有進(jìn)行定義，二是你把它理解成常數(shù)了。所以會(huì)報(bào)錯(cuò)，找不到該函數(shù)（你把int函數(shù)用錯(cuò)了）

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

7樓2015-11-16 10:17: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è): 自然語言理解與機(jī)器翻譯

函數(shù)調(diào)用有問題

回復(fù)此樓

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

8樓2015-11-16 11:41:51

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

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

[考研] 322求調(diào)劑 +3	X123- 2026-03-14	3/150	2026-03-21 18:40 by ColorlessPI
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 299求調(diào)劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊(duì) 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 323求調(diào)劑 +3	洼小桶 2026-03-18	3/150	2026-03-20 22:54 by JourneyLucky
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 生物學(xué)調(diào)劑招人！��！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999