24小时热门版块排行榜

>论坛更新日志 (3856)
>文献求助 (403)
>虫友互识 (328)
>导师招生 (313)
>博后之家 (155)
>硕博家园 (154)
>休闲灌水 (138)
>考博 (124)
>论文投稿 (110)
>基金申请 (90)
>招聘信息布告栏 (54)
>考研 (54)
>教师之家 (48)
>公派出国 (43)
>绿色求助(高悬赏) (37)
>找工作 (37)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

若为君故-

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 23.1
帖子: 13
在线: 3.2小时
虫号: 4000018
注册: 2015-07-31
专业: 金融学

[求助] 这个程序为什么一直出错已有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);

回复此楼

» 猜你喜欢

到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有7人回复
申请2026年博士已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有6人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

c程序设计伴侣 : 帮你更好地理解谭浩强老师的那本书以及更多！已经有168人回复
【2011最新美国犯罪纪录片】【凝望深渊/坠入深渊】【高清1280版BD-RMVB/中字】已经有31人回复
【分享】CAD十个实用问答已经有5人回复
【转帖】电脑用久了就会变慢。教你补救的方法已经有15人回复

1楼2015-09-11 21:51:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.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中乘除幂和点乘点除点幂的区别。

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)))')

最后，考虑到楼主方程组的复杂性，应该是找不到解析解的。