| 查看: 518 | 回复: 4 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
这个程序为什么一直出错已有2人参与
|
||||
|
编了一个程序但运行出来一直显示??? Error using ==> solve ' L1*（a1./(sqrt(1-a1.^2))）=L2*（a2./(sqrt(1-a2.^2))） ' is not a valid expression or equation. 是为什么呀? clc;clear all; syms x y; L1=1.1822;L2=1.2491;L3=1.2832;L4=1.3180; t1=14.75;t2=14.85;t3=14.9;t4=14.95; w1=[15*(t1-20+x/15)*pi]/180; w2=[15*(t2-20+x/15)*pi]/180; w3=[15*(t3-20+x/15)*pi]/180; w4=[15*(t4-20+x/15)*pi]/180; n=108; c=2*pi.*23.45/360; sgm=c.*sin((2.*pi.*(284+n))/365); a1=sin(y)*sin(sgm)+cos(y)*cos(sgm)*cos(w1); a2=sin(y)*sin(sgm)+cos(y)*cos(sgm)*cos(w2); a3=sin(y)*sin(sgm)+cos(y)*cos(sgm)*cos(w3); a4=sin(y)*sin(sgm)+cos(y)*cos(sgm)*cos(w4); solve('L1*（a1./(sqrt(1-a1.^2))）=L2*（a2./(sqrt(1-a2.^2))）','L3*（a3./(sqrt(1-a3.^2)）)=L4*（a4./(sqrt(1-a4.^2)）)'); x=vpa(x,1); y=vpa(y,1); |
回复此楼» 猜你喜欢
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有7人回复
-
申请2026年博士
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
c程序设计伴侣 : 帮你更好地理解谭浩强老师的那本书以及更多！
已经有168人回复- 【2011最新美国犯罪纪录片】【凝望深渊/坠入深渊】【高清1280版BD-RMVB/中字】
已经有31人回复
- 【分享】CAD十个实用问答
已经有5人回复
- 【转帖】电脑用久了就会变慢。教你补救的方法
已经有15人回复
3楼2015-09-14 18:19:16
2楼2015-09-12 09:27:06
4楼2015-09-15 08:56:18
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68120.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
|
MATLAB的solve求解方程，不需要移项把原方程整理为右边等于零的形式。楼主程序报错，与方程是否整理为右边等于零的形式无关。 楼主原程序存在以下两个错误： 其一，如果对MATLAB括号熟悉的话，一眼扫过就会发现，楼主原程序中： solve('L1*（a1./(sqrt(1-a1.^2))）=L2*（a2./(sqrt(1-a2.^2))）','L3*（a3./(sqrt(1-a3.^2)）)=L4*（a4./(sqrt(1-a4.^2)）)'); 出现了中文状态的括号，MATLAB是不能识别中文状态的字符的。 在MATLAB代码中出现中文字符，代码必然出错。 其二，MATLAB的solve函数是用于求解解析解的函数，方程中乘除幂的书写不需要也不能用点乘点除点幂的形式。可自行了解MATLAB中乘除幂和点乘点除点幂的区别。 最后，考虑到楼主方程组的复杂性，应该是找不到解析解的。 |
5楼2015-09-15 14:18:59