Question

matlab符号积分问题？

已知：ym=135.6;rk=511.1;rd=502.7;rm=403.5;re=333.1;rb=271.0;
t0=0.8;te=0.15;a=0.1;
b=rd^2-rm^2+(rk^2-rd^2)*(1-t0+a*rk/(rk-rd))-2*a*(rk^3-rd^3)/(3*(rk-rd));
g1=(rd^2-rm^2+(1-t0+a*rk/(rk-rd))*(r1.^2-rd^2)-2*a*(r1.^3-rd^3)./(3*(rk-rd)))./(b^2-(rd^2-rm^2+(1-t0+a*rk/(rk-rd))*(r1.^2-rd^2)-2*a*(r1.^3-rd^3)./(3*(rk-rd))).^2).^1/2；
syms y1 g1;
y1=int(g1 ,r1).
matlab 没写出表达式，也没提示错误。求帮助！
matlab R2018b ,WIN 7

Answer 1

题主给出符号积分的代码，运行matlab 没写出表达式？其主要的问题是：

1、syms 声明变量的位置，一般都放在执行代码最前面。

2、由于被积函数比较长，可能会得到莫名其妙的答案，但可以用vpa函数得到其解，即

y1=vpa(int(g1 ,r1))

3、修改后运行，可以得到如下结果。

追问

谢谢您的答案，求出了结果。但是求出rd对应的yd值不合理；也无法画出（rd<r1<rk）的曲线。请再帮下忙，谢谢！

追答

你后续给出的问题，缺已知参数（如a、τ0）.

追问

a和t0都和第一次问题中的一样，现在就是求rd对应的yd数值不符合实际值。请帮帮忙啊！

追答

根据你的公式重新编程，用integral数值积分函数得到y值。用for循环语句得到对应r的y值。

追问

请问，能否把你写的编码发给我，学习学习，我是刚入门啊。十分感谢！！

追答

可以的

追问

优质答主，请帮帮忙把程序发给我呗，谢谢！！

追答

昨日己发给你了，在私信上

追问

好的，谢谢！那个悬赏怎么给你啊？

Answer 2

先定义变量，再定义函数。
g1这里是函数，不要定义成变量
比如这里g1是关于自变量r1 的函数，试一试这么写，
syms r1;
g1= ....r1.....
int(g1,r1)

追问

还是不行啊！

追答

你把代码复制过来看看呢？