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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 用matlab求解方程式
111/2返回列表12下一页
查看: 3721  |  回复: 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下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭