应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 如何用matlab求解方程，请给出程序，问题解决再追加25金币
51/1返回列表
查看: 2625  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

晴天--啦啦

铁杆木虫 (正式写手)

[求助] 如何用matlab求解方程，请给出程序，问题解决再追加25金币 已有3人参与

sin(a)*d-[(b*cos(x-c)+sqrt(d^2-b^2*sin^2(x-c))*sin(x-c)]=0
其中：a=1:0.1*pi:2*pi  , b=0.02 , c=pi/6 , d=20
求未知向量x的值
回复此楼

» 猜你喜欢
等一个晴天
1楼 2022-01-26 23:20:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晴天--啦啦

铁杆木虫 (正式写手)
送红花一朵
引用回帖:
7楼: Originally posted by zyh700600 at 2022-01-31 09:15:00
a取值是0：2pi，得出来的x值范围不也应该是0：2pi么？




因为，您的方程是一个周期性函数，
所以，您的方程有多解。

设定不同的初值，可以得到不同的解。


1.这是趋向 0 的解：

    1.0000     ...

谢谢
回复此楼
等一个晴天
8楼2022-02-02 10:34:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

hollya27

至尊木虫 (文坛精英)

Katsuragi Megumi

【答案】应助回帖

感谢参与，应助指数 +1
a=1:0.1*pi:2*pi
b=0.02
c=pi/6
d=20
syms x

for ii=1:length(a)
    eqn= sin(a(ii))*d-(b*cos(x-c)+sqrt(d^2-b^2*sin(x-c).^2)).*sin(x-c)==0
    S{ii}= solve(eqn,x)
end
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

晴天--啦啦
2楼2022-01-27 01:00:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晴天--啦啦

铁杆木虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by hollya27 at 2022-01-27 01:00:54
a=1:0.1*pi:2*pi
b=0.02
c=pi/6
d=20
syms x

for ii=1:length(a)
    eqn= sin(a(ii))*d-(b*cos(x-c)+sqrt(d^2-b^2*sin(x-c).^2)).*sin(x-c)==0
    S{ii}= solve(eqn,x)
end

谢谢，不过求出来的结果为什么都是[2×1 sym]？我想要具体的数值
赞一下(1人) 回复此楼
等一个晴天
3楼2022-01-27 10:05:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyh700600

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
CODE:
b=0.02;c=pi/6;d=20;

eqn = @(a,x) sin(a)*d-(b*cos(x-c) + sqrt(d^2-b^2*sin(x-c).^2)).*sin(x-c);

options=optimset('Display','off');   % Option to display output
X = @(a) fsolve(@(x) eqn(a,x),0,options);  % Call optimizer

A = 1:0.02*pi:2*pi;  % a取值
XA = arrayfun(@(a) X(a),A);  % A对应的X

[A',XA']

plot(A,XA)   %绘图
赞一下 回复此楼
4楼2022-01-29 11:56:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭