| 查看: 577 | 回复: 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); |
回复此楼» 猜你喜欢
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有8人回复
-
一志愿华南师大 070300（化学）304分求调剂
已经有3人回复
-
一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。
已经有3人回复
-
一志愿重庆大学085700资源与环境专硕，总分308求调剂
已经有3人回复
-
354求调剂
已经有8人回复
-
304求调剂
已经有6人回复
-
一志愿西南交大，求调剂
已经有5人回复
-
299求调剂
已经有5人回复
-
材料专业求调剂
已经有6人回复
-
311求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- c程序设计伴侣 : 帮你更好地理解谭浩强老师的那本书以及更多！
已经有168人回复
- 【2011最新美国犯罪纪录片】【凝望深渊/坠入深渊】【高清1280版BD-RMVB/中字】
已经有31人回复
- 【分享】CAD十个实用问答
已经有5人回复
- 【转帖】电脑用久了就会变慢。教你补救的方法
已经有15人回复
2楼2015-09-12 09:27:06
3楼2015-09-14 18:19:16
4楼2015-09-15 08:56:18
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.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
