24小时热门版块排行榜    

查看: 753  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 若为君故*-* 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 4/200 2026-11-02 03:50 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-02 03:18 by xzcv0Tehd8lQ
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 xvJS1BYDxvBS 2026-10-31 12/600 2026-11-02 03:11 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-02 03:11 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:58 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-02 02:45 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:27 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 21:52 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
信息提示
请填处理意见