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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 用matlab求解方程式
汕头大学海洋科学接受调剂
111/2返回列表12下一页
查看: 3610  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

钟声入翠微

铜虫 (小有名气)

[求助] 用matlab求解方程式 已有1人参与

假设a和b已知，a,b,x,y>0,求解x和y
a=[1+(0.0172*x)^2]^(-y/2)
b=y*arctan(0.0172*x)
希望可以提供程序代码，或者解！非常感谢！
回复此楼

» 猜你喜欢
1楼 2021-07-05 16:17:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
钟声入翠微: 金币+40, ★★★★★最佳答案 2021-07-08 21:19:15
引用回帖:
9楼: Originally posted by 钟声入翠微 at 2021-07-08 14:23:57
谢谢！（30000,0.05）这个初始值是自己根据实际情况设置，还是随意设置的呢？这次我运行出来了，但是感觉结果不是我想象中的那样。
...

这个需要尝试，可以用 fsolve 函数随机赋值求解。

结果不理想就要考虑方程是否正确，或者a,b取值是否合理。

function example_0708
clear;
a=0.709133065;
b=0.08555;
x0 = [1,1];
x = fsolve(@root2d,x0);
fprintf('\tx=%.8f\n',x(1));
fprintf('\ty=%.8f\n',x(2));

function F = root2d(x)
F(1) = (1+(0.0172*x(1))^2)^(-x(2)/2)-a;
F(2) = x(2)*atan(0.0172*x(1))-b;
end

end
赞一下(1人) 回复此楼

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

钟声入翠微
数值计算
10楼2021-07-08 16:30:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zf015222
。

发自小木虫Android客户端
2楼2021-07-05 16:42:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
这个尽量不要做伸手党，求解符号方程组，看看solve函数用法

发自小木虫Android客户端
赞一下 回复此楼

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

钟声入翠微
数值计算
3楼2021-07-05 19:49:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
送红花一朵
引用回帖:
3楼: Originally posted by 独孤神宇 at 2021-07-05 19:49:57
这个尽量不要做伸手党，求解符号方程组，看看solve函数用法

我是去百度了solve用法，但是运行的时候一直显示错误，
我最先给a，b赋值进行计算，但是一直显示运算符错误
syms x y;
>>[x ,y]=solve(' 0.709133=[1+(0.0172*x)^2]^(-y/2)','0.08555=y*arctan(0.0172*x)');
然后使用了未赋值
clc,clear
syms a b y x
sola=solve(a==(1+(0.0172*x)^2)^(-y/2),x)   %待求解的变量是x
sol=solve(b==y*atan(0.0172*x),y)  %待求解的变量是y
还是运算不出想要的东西
因为只会一点点matlab，就一直解决不了问题。
赞一下 回复此楼
4楼2021-07-06 09:02:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 钟声入翠微 at 2021-07-06 09:02:45
我是去百度了solve用法，但是运行的时候一直显示错误，
我最先给a，b赋值进行计算，但是一直显示运算符错误
syms x y;
>>=solve(' 0.709133=^(-y/2)','0.08555=y*arctan(0.0172*x)');
然后使用 ...

clear;
syms a b
syms x y positive
eqns=[(1+(0.0172*x)^2)^(-y/2)-a==0, y*atan(0.0172*x)-b==0];
[X,Y]=solve(eqns,[x y])

这个可能是方程有点问题吧。 a  b  的具体数值是多少
赞一下 回复此楼

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

钟声入翠微
数值计算
5楼2021-07-06 12:51:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by 独孤神宇 at 2021-07-06 12:51:49
clear;
syms a b
syms x y positive
eqns=;
=solve(eqns,)

这个可能是方程有点问题吧。 a  b  的具体数值是多少...

0.709133065        0.08555
0.245947581        0.23527
0.679383641        0.63827
0.32593606        0.51032
0.208709677        0.53469

a和b的值如上，比较小，我也不清楚能不能解出来。
赞一下 回复此楼
6楼2021-07-06 14:42:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
送红花一朵
引用回帖:
5楼: Originally posted by 独孤神宇 at 2021-07-06 12:51:49
clear;
syms a b
syms x y positive
eqns=;
=solve(eqns,)

这个可能是方程有点问题吧。 a  b  的具体数值是多少...

我自己又代入了一下，得到这样的结果，是什么意思呢？

Warning: The solutions are parameterized by the symbols: z, z1. To include parameters and conditions in the solution, specify the
'ReturnConditions' option.
> In solve>warnIfParams (line 510)
  In solve (line 367)
Warning: The solutions are valid under the following conditions: 1/((1849*z^2)/6250000 + 1)^(z1/2) -
6387302814580195/9007199254740992 == 0 & z1*atan((43*z)/2500) - 1711/20000 == 0 & 0 < z & 0 < z1. To include parameters and
conditions in the solution, specify the 'ReturnConditions' option.
> In solve>warnIfParams (line 517)
  In solve (line 367)

X =

z


Y =

z1
赞一下 回复此楼
7楼2021-07-06 14:53:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by 钟声入翠微 at 2021-07-06 14:53:48
我自己又代入了一下，得到这样的结果，是什么意思呢？

Warning: The solutions are parameterized by the symbols: z, z1. To include parameters and conditions in the solution, specify the
'ReturnCondit ...

初始值变化很大，用 vpasolve 函数赋初始值

clear;
a=0.709133065;b=0.08555;
syms x y positive
eqns=[(1+(0.0172*x)^2)^(-y/2)-a==0, y*atan(0.0172*x)-b==0];
[X,Y]=vpasolve(eqns,[x y],[30000;0.05])
赞一下 回复此楼

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

钟声入翠微
数值计算
8楼2021-07-06 17:42:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
送红花一朵
引用回帖:
8楼: Originally posted by 独孤神宇 at 2021-07-06 17:42:15
初始值变化很大，用 vpasolve 函数赋初始值

clear;
a=0.709133065;b=0.08555;
syms x y positive
eqns=;
=vpasolve(eqns,,)...

谢谢！（30000,0.05）这个初始值是自己根据实际情况设置，还是随意设置的呢？这次我运行出来了，但是感觉结果不是我想象中的那样。

发自小木虫Android客户端
赞一下 回复此楼
9楼2021-07-08 14:23:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 钟声入翠微 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +6 关忆北. 2026-04-14 6/300 2026-04-15 13:32 by Llllss12
[考研] 297，工科调剂? +7 河南农业大学-能 2026-04-14 7/350 2026-04-15 13:26 by 西北望—风沙
[考研] 求调剂推荐 +5 小聂爱学习 2026-04-14 5/250 2026-04-15 09:01 by blunt110
[考研] 279学硕食品专业求调剂院校 20+6 孤独的狼爱吃羊 2026-04-12 28/1400 2026-04-14 15:44 by zs92450
[考研] 22专硕求调剂 +8 haoyun上岸 2026-04-11 10/500 2026-04-13 22:15 by zhq0425
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 085404 293求调剂 +9 勇远库爱314 2026-04-08 9/450 2026-04-12 02:24 by 秋豆菜芽
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 08:57 by zhq0425
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 一志愿矿大，材料工程专硕314分，0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
信息提示
关闭
请填处理意见
关闭