lkcool7

新蟲 (初入文壇)

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

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

我自定義了三個函數(shù)Y1()，Y2()，Y3()，三個函數(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

運行是程序在調(diào)用Y1()時出錯，提示為

??? 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ù)時出錯，為什么會這樣啊？？？？？如果我將a=sqrt(pA*eval(int(u^2,x,0,l)));提前算出，直接將這行替換為a=169.2216;卻可以運行，并最后繪制出我要的曲線。但是在我自己定義的函數(shù)中調(diào)用int()居然會出錯？求大神幫助！謝謝

回復(fù)此樓

» 猜你喜歡

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

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

運行matlab程序出現(xiàn)警告：TIFF庫警告，怎么解決啊已經(jīng)有0人回復(fù)
matlab程序運行報錯求解已經(jīng)有2人回復(fù)
已有matlab程序如何導(dǎo)入excel數(shù)據(jù)運行已經(jīng)有9人回復(fù)
已有matlab程序如何導(dǎo)入excel數(shù)據(jù)運行已經(jīng)有2人回復(fù)
求助！哪位大神幫忙運行下MATLAB程序，求結(jié)果！已經(jīng)有3人回復(fù)
matlab程序運行報錯求助已經(jīng)有7人回復(fù)
求助MATLAB程序已經(jīng)有5人回復(fù)
Matlab中如何運行C++程序已經(jīng)有3人回復(fù)
matlab程序已經(jīng)有7人回復(fù)
Matlab程序運行結(jié)束提示已經(jīng)有5人回復(fù)
MATLAB程序調(diào)試已經(jīng)有9人回復(fù)
MATLAB程序調(diào)整已經(jīng)有4人回復(fù)
matlab程序運行時快時慢已經(jīng)有8人回復(fù)
一個MATLAB程序，運行不出來已經(jīng)有1人回復(fù)
一個matlab程序，運行不出來已經(jīng)有2人回復(fù)
matlab程序運行出錯，求指點修改已經(jīng)有4人回復(fù)
matlab程序只運行一次就跳出的問題已經(jīng)有2人回復(fù)
matlab程序運行中出現(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小時
蟲號: 2956857
注冊: 2014-02-06
專業(yè): 自動化

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

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

你的u是一個數(shù)，不是一個表達式，對一個數(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小時
蟲號: 3859120
注冊: 2015-05-09
專業(yè): 自動化

引用回帖:

2樓: Originally posted by 轉(zhuǎn)基因奔奔 at 2015-11-13 09:31:55
你的u是一個數(shù)，不是一個表達式，對一個數(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));
這個是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小時
蟲號: 2956857
注冊: 2014-02-06
專業(yè): 自動化

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

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

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

贊一下

回復(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小時
蟲號: 1227944
注冊: 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小時
蟲號: 3859120
注冊: 2015-05-09
專業(yè): 自動化

引用回帖:

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

謝謝回答，我已解決了問題，確實發(fā)覺函數(shù)Y1(x)中的u為數(shù)值而非符號表達式。然后我建立了符號變量 syms x k l 將int(u,x,0,l) 計算出來
下面是我解決問題后的代碼
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為符號積分int(u^2,x,0,l)結(jié)果，為符號表達式
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小時
蟲號: 2633300
注冊: 2013-09-03
專業(yè): 計算機科學(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));
這個是u的定義，是以x為自變量的函數(shù)啊...

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

發(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小時
蟲號: 2986702
注冊: 2014-02-23
性別: GG
專業(yè): 自然語言理解與機器翻譯

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

回復(fù)此樓

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

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

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

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

[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +14	曦熙兮 2026-03-15	14/700	2026-03-21 19:12 by ColorlessPI
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 求助 +5	夢里的無言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3	李軼男003 2026-03-20	3/150	2026-03-21 13:28 by 搏擊518
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	5/250	2026-03-20 22:11 by 云游重陽
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染