查看: 171|回复: 4|关注: 0

[已解决] 新手求教,integral求积分的问题

[复制链接]

新手

7 麦片

财富积分


050


1

主题

4

帖子

0

最佳答案
不知道我这么写有什么问题?

k=0.04;
b=0.212;
a=2;
R=5.46;
w=4.3175;
A=@(r)2*atan((exp(k/w)-exp(-k/w))/(exp(r/w)-exp(-r/w)));
options=optimset('Display','iter');
r0=fsolve(A,[0;50],options);%此处采用匿名函数法@(x)
%V=pow(m-1,2)+1-2*m*(a*(cos(A)-1)/(pow(r,2)*K)-b*sin(a)/(r*K));
for m = -50:1:50;

F=@(r)sqrt(k*(20-((m-1)^2+1-2*m*(a*(cos(2*atan((exp(k/w)-exp(-k/w))/(exp(r/w)-exp(-r/w))))-1)/(r^2*k)-b*sin(a)/(r*k))))-sqrt(19*k));
Q = integral(F,r0,50000.0);
end

回复主题 已获打赏: 0 积分

举报

新手

7 麦片

财富积分


050


1

主题

4

帖子

0

最佳答案
 楼主| 发表于 2020-11-18 15:13:05 | 显示全部楼层
:hug::hug::hug::hug:

新手

7 麦片

财富积分


050


1

主题

4

帖子

0

最佳答案
 楼主| 发表于 2020-11-18 20:35:00 | 显示全部楼层
:(:(:(:(:(

论坛优秀回答者

5

主题

2623

帖子

810

最佳答案
  • 关注者: 220
发表于 2020-11-19 00:03:57 | 显示全部楼层 |此回复为最佳答案
1.有补充信息,请直接在原帖修改就好;请不要发类似"顶顶顶"之类的无意义回帖,这种灌水回帖除了让他人心生反感之外,对解决你的问题没有任何帮助;

2.MATLAB最好的老师就是帮助文档,建议你逐行运行自己的代码并对照帮助文档检查语法。比如
r0=fsolve(A,[0;50],options)
这句就有问题,你的A是个单一式子的函数,按语法规则x = fsolve(fun,x0,options)给定初值x0只要一个具体数值就够了,你写的是一个向量,这就不合适

3.你这个A函数的图像我大致看了下,类似倒数 1/x的图像,只有自变量在正负无穷处才能令函数极限为零,所以求解r0没什么意义

4.for循环里面求数值积分的函数F,尽量用点乘点除点幂…

5.循环中m数值会变化,F就会变化,每步求出来的Q值可能是不同的,一般要按每步取值进行保存,比如循环外面先写个ii = 1;然后循环中每步令Q( ii ) = 具体计算; ii = ii + 1;
提问请:①准确描述问题②提出你的思考(等着抄作业的一律锁帖)③提供代码文本而非截图④及时反馈
回复此楼 已获打赏: 0 积分

举报

新手

7 麦片

财富积分


050


1

主题

4

帖子

0

最佳答案
 楼主| 发表于 2020-11-19 11:19:42 | 显示全部楼层
TouAkira 发表于 2020-11-19 00:03
1.有补充信息,请直接在原帖修改就好;请不要发类似"顶顶顶"之类的无意义回帖,这种灌水回帖除了让他人心生 ...

谢谢您的指导,对我之前由于自私的以自我为中心的行为深感歉意,在接下来的时间内,我会遵守规定,不发无意义的内容。
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /5 下一条

快速回复 返回顶部 返回列表