| 查看: 1053 | 回复: 5 | |||
[求助]
matlab解超越方程组求助！
|
|
dy = [ y(1)*(2*7.6e-23*y(2)*0.006-2*2.3e-21*y(3)*0.00058-2*1e-20*(1.2e22-y(3))*0.00058-2.64e-20*y(1)-(log(1/0.85)+0.1))/1.56e-9; -4*7.6e-23*y(1)*y(2)/(1.88e-19); 4*y(1)*(-2.3e-21*y(3)+1e-20*(1.2e22-y(3)))/(1.88e-19)]; dy是对t求导， t的范围在0-1e-7或者0-1e-8; 限制：y(1) = max(y(1),1); 初始条件：y(1) = 1; y(2) = 0; y(3) = 0; 我整理方程得到了这种方程组，不会解 麻烦高手们把Matlab程序等解题过程写的详细一些！谢谢！ 为了防止上述方程描述错误 我在下面附件用word把方程写出来了！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : matlab解超越方程组求助.docx
2015-05-26 10:29:20, 12.03 K
» 猜你喜欢
德国亥姆霍兹Hereon中心招收两位医用镁合金腐蚀与LPSO相变方向2026公派博士生
已经有0人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，推荐试试
已经有11人回复
-
物理学I论文润色/翻译怎么收费?
已经有204人回复
-
2026-CJ开始申报了
已经有1人回复
-
西安电子科学大学杭州研究院刘丽香教授招收智能多模态传感器和微型储能器件方向博士
已经有13人回复
-
重庆交通大学光子学微结构与器件课题组招收2026年硕士研究生信息
已经有0人回复
-
一志愿郑大材料学硕298分，求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 新人求助解MATLAB超越方程，方程中含有复数
已经有9人回复
- matlab对含贝塞尔函数的超越方程求根
已经有5人回复
- 求助：超越方程的解法
已经有13人回复
- 求解一个方程组，急！！！
已经有19人回复
- matlab如何求解多元多次方程，包括绘图
已经有7人回复
- matlab解方程组求助攻
已经有4人回复
- 多元一次方程的matlab解法，要求解为正整数
已经有10人回复
- 求助用matlab或者（1st0pt)编程解方程组
已经有3人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- matlab解偏微分方程求助
已经有12人回复
- matlab求解高阶方程组
已经有3人回复
- 有偿求助：matlab怎么解联立方程组？
已经有11人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- matlab求解色散方程，求助
已经有12人回复
- 求助Matlab解方程问题
已经有6人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
- 【求助】复超越方程的解法
已经有19人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-05-26 19:00:43
3楼2015-05-26 21:53:10
|
我把原来的式子、参数、和应该得到的图形上传上来了，并附我写的程序，麻烦请看一下哪里出了问题，为什么得不到想要的结果，请大神写清具体步骤！求赐教！并附全部金币，只有8个了。  参数1.png  参数2.jpg  参数3.jpg  待解式子.png  结果图（应该为I的分布图）.png |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 参数和代码.docx
2015-05-27 19:01:30, 13.03 K
4楼2015-05-27 19:04:19
|
再贴一下文件中代码，方面观看。 lc clear close all y0 = [1e9;4e20;0]; tspan=[0 1e-7]; tic [t,y] = ode23('rate_eq',tspan,y0); toc y(:,1) = max(y(:,1),1); figure subplot(3,1,1); plot(t,y(:,1)); xlabel(' 时间(s)'); ylabel('光子数密度(m^{-3})'); subplot(3,1,2); plot(t,y(:,2)); xlabel('时间(s)'); ylabel('反转粒子束密度(m^{-3})'); subplot(3,1,3); plot(t,y(:,3)); xlabel('时间(s)'); ylabel('基态粒子数密度{-3})'); function dy = rate_eq(t,y) sigma = 7.6e-23; sigma_o = 1e-20; sigma_s = 2.3e-21; n1 = 1.972; n2 = 3.48; L = 0.1; l = 0.006; d = 0.00058; R = 0.85; c = 2.997963e8; lc = n1*l+n2*d+(0.23-l-d); tr = 2*lc/c; n0 = 1.2e22; beta = 2.6e-10;h = 6.623e-34;lamda = 1.06e-6;niu = c/lamda; B = 6*beta*h*niu*c*(0.3)^2;%此项可以为0； y(1) = max(y(1),1); dy = [ y(1)*(2*sigma*y(2)*l-2*sigma_s*y(3)*d-2*sigma_o*(n0-y(3))*d-B*y(1)-(log(1/R)+L))/tr; -4*sigma*y(1)*y(2)/(h*niu); 4*y(1)*(sigma_o*(n0-y(3))-sigma_s*y(3))/(h*niu)]; 请虫友帮忙！ |
5楼2015-05-27 19:08:59
6楼2015-05-28 10:31:24
