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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 用matlab求解方程式
111/2返回列表12下一页
查看: 3886  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +20 Lanmanbaby 2026-08-09 33/1650 2026-08-12 11:26 by lantianda
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
关闭
请填处理意见
关闭