| 查看: 1077 | 回复: 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
» 猜你喜欢
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有256人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有11人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 新人求助解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
- 金币: 68129.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
