| 查看: 1148 | 回复: 3 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
求怎么用matlab解这个方程 已有2人参与
|
|||||
|
要解个方程，但我对matlab不熟，请各位大神帮写个程序代码？ 0.5sin(3.6x)cos(x)+3.6cos(3.6x)sin(x)=0 （x>0) 有无穷多解，求100个正的解。 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有229人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用matlab求解 多个y值情况下对应的每个x的解
已经有6人回复
- 如何用matlab解符号方程组
已经有4人回复
- MATLAB编程解方程组并用simulink进行仿真处理
已经有3人回复
- matlab对含贝塞尔函数的超越方程求根
已经有5人回复
- 求如何用有限差分法解matlab中偏微分方程与微分方程联立的方程组
已经有2人回复
- MATLAB求二次规划的最优解
已经有11人回复
- 求助各位MATLAB微分方程互相迭代怎么做？
已经有0人回复
- 用matlab求解复数方程
已经有5人回复
- MATLAB怎么求以下微分方程
已经有0人回复
- matlab解微分方程的问题
已经有2人回复
- matlab 求解 一个方程
已经有6人回复
- matlab求解非线性方程组
已经有3人回复
- 求各位帮我看看这个方程组用matlab怎么解
已经有0人回复
- MATLAB解二阶常微分方程组
已经有0人回复
- 常微分方程解析解
已经有5人回复
- matlab ode15s方程中两个参数怎么调用其他函数带入求解
已经有6人回复
- matlab怎么求二阶常微分方程复杂边界情况的数值解？
已经有0人回复
- 请问用matlab数值方法求解X^5=10的解，怎么求？
已经有1人回复
- matlab计算动力学方程求代码指点
已经有0人回复
- 请教高手：PDE工具箱能解两个偏微分方程相互迭代求解吗
已经有2人回复
- 【求助】用MATLAB 求微分方程
已经有18人回复
- 【求助】用matlab解含第二类边界条件的偏微分方程
已经有1人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68130.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jyypzfswty: 金币+15, ★★★★★最佳答案, 牛人~~灰常感谢~ 2015-12-11 09:09:25
jyypzfswty: 金币+15, ★★★★★最佳答案, 牛人~~灰常感谢~ 2015-12-11 09:09:25
|
用MATLAB的fsolve函数也是可以做的。 代码如下： 找到了147个解，零点和多余的，请自行删除。 x_answer = Columns 1 through 15 0 0.5047 1.3189 2.1562 2.9077 3.2761 3.9627 4.7964 5.6231 6.2228 6.6432 7.4348 8.2731 9.0648 9.4852 Columns 16 through 30 10.0848 10.9116 11.7453 12.4319 12.8003 13.5518 14.3890 15.2032 15.7080 16.2127 17.0269 17.8641 18.6157 18.9841 19.6707 Columns 31 through 45 20.5044 21.3311 21.9307 22.3511 23.1428 23.9811 24.7727 25.1932 25.7928 26.6195 27.4532 28.1398 28.5082 29.2598 30.0970 Columns 46 through 60 30.9112 31.4159 31.9207 32.7348 33.5721 34.3236 34.6920 35.3786 36.2123 37.0391 37.6387 38.0591 38.8507 39.6891 40.4807 Columns 61 through 75 40.9011 41.5008 42.3275 43.1612 43.8478 44.2162 44.9677 45.8050 46.6191 47.1239 47.6286 48.4428 49.2801 50.0316 50.4000 Columns 76 through 90 51.0866 51.9203 52.7470 53.3467 53.7671 54.5587 55.3970 56.1887 56.6091 57.2087 58.0354 58.8691 59.5558 59.9241 60.6757 Columns 91 through 105 61.5129 62.3271 62.8319 63.3366 64.1508 64.9880 65.7396 66.1080 66.7946 67.6283 68.4550 69.0546 69.4750 70.2667 71.1050 Columns 106 through 120 71.8966 72.3171 72.9167 73.7434 74.5771 75.2637 75.6321 76.3837 77.2209 78.0351 78.5398 79.0446 79.8587 80.6960 81.4475 Columns 121 through 135 81.8159 82.5025 83.3362 84.1630 84.7626 85.1830 85.9746 86.8130 87.6046 88.0250 88.6246 89.4514 90.2851 90.9717 91.3401 Columns 136 through 147 92.0916 92.9289 93.7430 94.2478 94.7525 95.5667 96.4039 97.1555 97.5239 98.2105 99.0442 99.8709  附图1.png |
4楼2015-12-10 16:28:32
Op_liu
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 3114.7
- 散金: 51
- 红花: 14
- 帖子: 773
- 在线: 214.8小时
- 虫号: 3482925
- 注册: 2014-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2015-12-10 13:43:13
Op_liu
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 3114.7
- 散金: 51
- 红花: 14
- 帖子: 773
- 在线: 214.8小时
- 虫号: 3482925
- 注册: 2014-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jyypzfswty: 金币+15, ★★★★★最佳答案, 灰常感谢~ 2015-12-11 09:08:33
jyypzfswty: 金币+15, ★★★★★最佳答案, 灰常感谢~ 2015-12-11 09:08:33
|
%%这是非线性方程求解，普通的solve无法求得准确的解，fsolve无法解得100个解； %%我这里用的是数值解法，精确度在为10^(-4)，可改变x来提高精确度，看个人需要； %%思路：求出y=0.5*sin(3.6*x).*cos(x)+3.6*cos(3.6*x).*sin(x)=0的100个解X； %可画出其图像，很容易看出y会在解处变号！！！由此可筛选出解 x=0:0.0001:70;%间隔越小，精确度越高 y=0.5*sin(3.6*x).*cos(x)+3.6*cos(3.6*x).*sin(x);%函数 plot(x,y,'b') hold on n=1; for i=1  length(x)-1)if y(i)*y(i+1)<=0%如果小于0，知道此处y发生变号，那么此处对应的x即解 X(n)=x(i);%将找到的解赋值给另一个矩阵 n=n+1; end end X=X(2:101);%因为问题说前100个正解，不包括0 Y=0.5*sin(3.6*X).*cos(X)+3.6*cos(3.6*X).*sin(X); plot(X,Y,'.r')%画出解的点 |
3楼2015-12-10 13:43:57
