| 查看: 662 | 回复: 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); |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有12人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 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
100