24小时热门版块排行榜

>论坛更新日志 (4421)
>文献求助 (363)
>虫友互识 (363)
>导师招生 (260)
>论文投稿 (173)
>博后之家 (139)
>硕博家园 (132)
>考博 (100)
>招聘信息布告栏 (99)
>休闲灌水 (96)
>基金申请 (70)
>考研 (67)
>论文道贺祈福 (57)
>教师之家 (55)
>公派出国 (48)
>高分子 (39)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lkcool7

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 87.5
散金: 3
帖子: 22
在线: 6.1小时
虫号: 3859120
注册: 2015-05-09
专业: 自动化

[求助] 运行matlab程序中自定义函数时出问题了，还请大神帮忙！谢谢已有1人参与

我自定义了三个函数Y1()，Y2()，Y3()，三个函数只是内部参数k发生变
在这里我给出Y1()代码

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));

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

运行是程序在调用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（）中调用int积分函数时出错，为什么会这样啊？？？？？如果我将a=sqrt(pA*eval(int(u^2,x,0,l)));提前算出，直接将这行替换为a=169.2216;却可以运行，并最后绘制出我要的曲线。但是在我自己定义的函数中调用int()居然会出错？求大神帮助！谢谢

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有50人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lkcool7

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 87.5
散金: 3
帖子: 22
在线: 6.1小时
虫号: 3859120
注册: 2015-05-09
专业: 自动化

引用回帖:

4楼: Originally posted by 转基因奔奔 at 2015-11-13 16:52:07
你的x是以输入变量的形式带入的，每次循环都是以一个数作为x输入的。MATLAB会自动解算出u的数值，而不是把u当做一个表达式

谢谢回答，我已解决了问题，确实发觉函数Y1(x)中的u为数值而非符号表达式。然后我建立了符号变量 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)结果，为符号表达式
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;