24小时热门版块排行榜    

查看: 570  |  回复: 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);
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxj126

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
移项 转换成equation=0的形式
2楼2015-09-12 09:27:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

若为君故*-*

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by fxj126 at 2015-09-12 09:27:06
移项 转换成equation=0的形式

可以具体说一下怎么改吗?
3楼2015-09-14 18:19:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxj126

木虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 若为君故*-* at 2015-09-14 18:19:16
可以具体说一下怎么改吗?...

就是移项 等式左边=0
4楼2015-09-15 08:56:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

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中乘除幂和点乘点除点幂的区别。
CODE:
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);
[x y]=solve('L1*(a1/(sqrt(1-a1^2)))=L2*9*a2/(sqrt(1-a2^2))','L3*(a3/(sqrt(1-a3^2)))=L4*(a4/(sqrt(1-a4^2)))')

最后,考虑到楼主方程组的复杂性,应该是找不到解析解的。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
5楼2015-09-15 14:18:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 若为君故*-* 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 Mqqqqqq 2026-03-19 4/200 2026-03-20 14:15 by 星空星月
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by 呱呱呱呱叫
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +3 一瓶苯甲酸 2026-03-14 3/150 2026-03-19 15:17 by 尽舜尧1
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 一志愿天大材料与化工(085600)总分338 +5 蔡大美女 2026-03-13 5/250 2026-03-19 10:44 by 是小刘呀~
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
信息提示
请填处理意见