找回密码
 注册
关于网站域名变更的通知
查看: 584|回复: 2
打印 上一主题 下一主题

请问程序应该如何写?感谢!

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-6-19 11:30 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
如图所示,要在matlab中求这种非线性较强的积分,做出(x,t)的图像,程序应该如何写?我调用的quadv数值积分的函数,也没有成功,麻烦大神指点一下,非常感谢!
$ s1 E4 V  x" U! g0 K
4 G& P5 p* b0 r# b

该用户从未签到

2#
发表于 2020-6-19 13:27 | 只看该作者
8 e2 b& B$ v& q
可以用integral函数
9 p& ?+ j) o& _5 s) o
  • syms t real
  • t_2 = 0.0976;
  • L = 0.3;
  • Ee = 2.068e+11;
  • x = 0:1e-3:0.2135;
  • t_1 = x - 0.0226;
  • result = zeros(size(x));
  • result = arrayfun(@(ii) integral(@(t) (t_2-t).*cos(t).*( sin(t_1(ii)) ).^2./(2.*Ee.*L.*( sin(t) + (t_2 - t).*cos(t) )), t_1(ii), t_2 ), [1:1:numel(x)]);
  • plot( x, result, 'ko-')
    9 h9 a$ O  R0 y, \0 s' h0 R
& u. P4 h, R( G& T* Z' e: m! a

点评

大神太牛了,非常感谢!  详情 回复 发表于 2020-6-19 18:10

该用户从未签到

3#
 楼主| 发表于 2020-6-19 18:10 | 只看该作者
leleeda 发表于 2020-6-19 13:27; q, A+ D3 K0 L
可以用integral函数, ~6 }" b  E1 j
  • syms t real
  • $ l# `3 `( Y( f/ x
    大神太牛了,非常感谢!! J4 L$ \$ |. w1 p3 l
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-11-24 13:09 , Processed in 0.156250 second(s), 28 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

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