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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 用matlab求解方程式
111/2返回列表12下一页
查看: 3458  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中南大学理学化学 +6 15779376950 2026-03-01 8/400 2026-03-04 16:19 by zhukairuo
[考研] 301求调剂 +3 李LJR 2026-03-04 3/150 2026-03-04 16:15 by sslc1985
[考研] 一志愿314求调剂 +7 202111120625 2026-03-03 7/350 2026-03-04 15:56 by zhukairuo
[考研] 能动297求调剂，本科川大 +4 邵11 2026-03-04 4/200 2026-03-04 15:29 by lylylyheart
[考研] 化学 0703求调剂 总分293 一志愿211 +7 土土小虫 2026-03-03 7/350 2026-03-04 10:15 by qingfeng258
[考研] 材料工程269求调剂 +7 白刺玫 2026-03-02 7/350 2026-03-04 08:47 by zhyzzh
[考研] 0854总分272 +5 打小就是老实人 2026-03-02 6/300 2026-03-04 01:41 by ouhaiyu
[考研] 290求调剂 +9 ErMiao1020 2026-03-02 9/450 2026-03-03 18:03 by linlonghao
[考研] 266材料化工求调剂 +3 哇塞王帅 2026-03-03 3/150 2026-03-03 17:58 by houyaoxu
[考研] 计算机学硕 分数285求调剂 +4 glwshine 2026-03-02 5/250 2026-03-03 14:27 by king呀
[考研] 281电子信息求调剂 +3 jhtfeybgj 2026-03-02 6/300 2026-03-03 12:12 by 热情沙漠
[考研] 290分材料工程085601求调剂 数二英一 +8 llx0610 2026-03-02 9/450 2026-03-02 22:09 by 无际的草原
[考研] 085600材料工程一志愿中科大总分312求调剂 +9 吃宵夜1 2026-02-28 11/550 2026-03-02 20:14 by hypershenger
[考研] 材料284求调剂，一志愿郑州大学英一数二专硕 +15 想上岸的土拨鼠 2026-02-28 15/750 2026-03-02 20:13 by hypershenger
[考研] 303求调剂 +5 今夏不夏 2026-03-01 5/250 2026-03-02 15:01 by 向上的胖东
[考博] 诚招农业博士 +3 心欣向荣 2026-02-28 3/150 2026-03-02 13:33 by 时间不狗
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
信息提示
关闭
请填处理意见
关闭