24小时热门版块排行榜    

查看: 572  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 若为君故*-* 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +11 墨墨漠 2026-03-18 12/600 2026-03-20 19:43 by 丁丁*
[考研] 265求调剂 +10 梁梁校校 2026-03-19 10/500 2026-03-20 19:42 by 丁丁*
[考研] 0703化学调剂 ,六级已过,有科研经历 +13 曦熙兮 2026-03-15 13/650 2026-03-20 19:35 by Dream007008
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +3 solanXXX 2026-03-20 3/150 2026-03-20 16:41 by fxue1114
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 材料学硕297已过四六级求调剂推荐 +6 adaie 2026-03-19 6/300 2026-03-20 14:12 by 星空星月
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +3 葵梓卫队 2026-03-18 5/250 2026-03-19 19:35 by 给你你注意休息
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 一志愿中海洋材料工程专硕330分求调剂 +7 小材化本科 2026-03-18 7/350 2026-03-19 10:46 by Linda Hu
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
信息提示
请填处理意见