MATLAB编写程序不知道错哪里了

clcclearsymsyy0detamuQEMNy0=0;deta=1.14;mu=1.14;Q=-0.2;E=y0^2*(deta^2-1)+1/(4*y0)*log... clc
clear
syms y y0 deta mu Q E M N
y0=0;
deta =1.14;
mu = 1.14;
Q=-0.2;
E=y0^2*(deta^2 - 1)+1/(4*y0)*log((deta^2 + y0 - 1)/(y0*deta^2))+(deta^2 - 1)/(2*y0)-1/4*log((deta^2 + y0 - 1)/y0)- mu*((y0/2*log(deta))-(deta^2 - 1)/(4*y0^2)-y0*(deta^2 - 1)/2+y0*2/4*log((deta^2 + y0 - 1)/y0)-y0^2/2*log(deta)-y/4*log((deta^2 + y0 - 1)/y0))-Q*y0;
N = y*(deta^2 - 1)*(deta^2 - 1) - 2*y^2*(deta^2 - 1) + 2*y^3*log((deta^2 + y - 1)/y) + mu*((4*log(deta)) + (deta^2 - 1)*(deta^2 - 1) - 2*(deta^2 - 1));
M = 8*y^6*(deta^2 - 1) + 8*y^3*log((deta^2 + y - 1)/(y*deta^2)) + 16*y^3*(deta^2 - 1) - 8*y^4*log((deta^2 + y - 1)/y) - 32*Q*y^5 - 32*y^4*E + mu*(16*y^6*log(deta) - 16*y^5*log(deta) + 16*y^5*(deta^2 - 1) + 8*y^2*(deta^2 - 1) - 8*y^6*log((deta^2 + y - 1)/y) + 8*y^5*log((deta^2 + y - 1)/y));
f=sqrt(N/M);
V=int(f,y,0,1);
eval(V);
运行之后的结果不是ans V,而是f,并且f里面居然含有x,求大神给指点,纠结老久了
展开
 我来答
lhmhz
高粉答主

2017-07-23 · 专注matlab等在各领域中的应用。
lhmhz
采纳数:7264 获赞数:17014

向TA提问 私信TA
展开全部

经运行你编写的MATLAB程序,存在如下问题:

1、由于y0=0,而在E表达式的分母有y0存在,所以E表达式是无意义的。

2、由于你的f表达式比较复杂,建议你用quad()数值积分函数求解。所以你的程序因做如下调整

f=sqrt(N/M);

func=inline(f,'y')

V= quad(@(y)func(y),0,1)

如,y0=0.0001时,运行程序可得

V =   2.9191 + 0.1230i

从运行结果来看,可能你的输入的公式有点问题。

纷扰不过闲事i
2017-07-23 · TA获得超过100个赞
知道答主
回答量:28
采纳率:0%
帮助的人:7.2万
展开全部
你这个程序到底是算什么呢?
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式