版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

LESK乐

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1115.9
帖子: 58
在线: 123.1小时
虫号: 2698688
注册: 2013-10-04
专业: 电力电子学

[求助] 用matlab解非线性方程组

如题，方程组在下面，x1 x2 x3 x4的初值非别为0.2， 0.6，0.75，1.45，希望大神帮忙解答
2*(cos(x1)-cos(x2)+cos(x3)-cos(x4))-1=π/4（等号右边是4分之pai）
2*(cos(3*x1)-cos(3*x2)+cos(3*x3)-cos(3*x4))-1=0
2*(cos(9*x1)-cos(9*x2)+cos(9*x3)-cos(9*x4))-1=0
2*(cos(15*x1)-cos(15*x2)+cos(15*x3)-cos(15*x4)))-1=0

[ Last edited by LESK乐 on 2013-10-21 at 14:24 ]

回复此楼

» 猜你喜欢

求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有5人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
需要合成515-64-0，50g，能接单的留言已经有4人回复
自荐读博已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab数值计算时迭代项分母为零的时候怎么办？已经有15人回复
matlab 解方程组并作图 -----望高人指导。。。已经有3人回复
用matlab解非线性方程组编程已经有2人回复
用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？已经有11人回复
matlab如何求解一个非线性微分方程组已经有8人回复
MATLAB求解非线性方程组已经有5人回复
求Matlab解非线性方程代码已经有12人回复
有偿求助：matlab怎么解联立方程组？已经有11人回复
求Matlab高手解决线性方程组的迭代求解问题已经有17人回复
求fortran解非线性方程组已经有4人回复
matlab求解非线性方程组，错误提示怎么解决已经有5人回复
求助 MATLAB解方程组-fslove 已经有7人回复
求助matlab---fsolve解非线性方程组已经有6人回复
matlab中怎么取方程组的一组解并画图已经有5人回复
求助MATLAB解方程组已经有6人回复
matlab求解非线性方程组已经有16人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！已经有18人回复
【求助】非线性方程组求解已经有3人回复
【求助】求教matlab解非线性方程组已经有9人回复
【求助】求一个用牛顿迭代法解非线性方程组的程序已经有10人回复

1楼 2013-10-21 14:22:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhg19

金虫 (小有名气)

应助: 33 (小学生)
金币: 894.6
红花: 3
帖子: 212
在线: 114.5小时
虫号: 1895765
注册: 2012-07-16
性别: GG
专业: 半导体物理

【答案】应助回帖

感谢参与，应助指数 +1

编写xfsolve.m文件：
CODE:
function xfsolve
clear all;clc
x0=[0.2,0.6,0.75,1.45];
x1=fsolve(@abc,x0);
function f=abc(x)
f(1)=2*(cos(x(1))-cos(x(2))+cos(x(3))-cos(x(4)))-1-pi/4;
f(2)=2*(cos(3*x(1))-cos(3*x(2))+cos(3*x(3))-cos(3*x(4)))-1;
f(3)=2*(cos(9*x(1))-cos(9*x(2))+cos(9*x(3))-cos(9*x(4)))-1;
f(4)=2*(cos(15*x(1))-cos(15*x(2))+cos(15*x(3))-cos(15*x(4)))-1;

赞一下

回复此楼

2楼2013-10-21 15:30:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68121.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

CODE:

function solvenonl_Eqs

clear all;clc

x0=[0.2 0.6 0.75 1.45];

x=fsolve(@Eqs,x0)

fval=Eqs(x)

function f=Eqs(x)

x1=x(1);

x2=x(2);

x3=x(3);

x4=x(4);

f(1)=2*(cos(x1)-cos(x2)+cos(x3)-cos(x4))-1-pi/4;

f(2)=2*(cos(3*x1)-cos(3*x2)+cos(3*x3)-cos(3*x4))-1;

f(3)=2*(cos(9*x1)-cos(9*x2)+cos(9*x3)-cos(9*x4))-1;

f(4)=2*(cos(15*x1)-cos(15*x2)+cos(15*x3)-cos(15*x4))-1;

结果：

x =

0.1499 0.5959 0.7473 1.5688

残差fval =

1.0e-012 *

-0.0022 -0.0183 -0.0450 0.3368

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

3楼2013-10-21 15:35:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

LESK乐

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1115.9
帖子: 58
在线: 123.1小时
虫号: 2698688
注册: 2013-10-04
专业: 电力电子学

引用回帖:

3楼: Originally posted by 月只蓝 at 2013-10-21 15:35:54
function solvenonl_Eqs
clear all;clc
x0=;
x=fsolve(@Eqs,x0)
fval=Eqs(x)

function f=Eqs(x)
x1=x(1);
x2=x(2);
x3=x(3);
x4=x(4);
f(1)=2*(cos(x1)-cos(x2)+cos(x3)-cos(x4))-1-pi/4;
f(2)=2*(cos ...

谢谢，不过为什么我用计算器把答案带进去算出来的方程答案不对啊，http://www.ilovematlab.cn/thread-263705-1-1.html这是我之前发的帖子，里面的结果跟你的一样

赞一下

回复此楼

4楼2013-10-21 16:04:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68121.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
LESK乐: 金币+20, ★★★★★最佳答案, 恩，谢谢你 2013-10-21 19:46:38

引用回帖:

4楼: Originally posted by LESK乐 at 2013-10-21 16:04:49
谢谢，不过为什么我用计算器把答案带进去算出来的方程答案不对啊，http://www.ilovematlab.cn/thread-263705-1-1.html这是我之前发的帖子，里面的结果跟你的一样...

我之前回复的残差fval=1.0e-012 *

-0.0022 -0.0183 -0.0450 0.3368
就是把x代回方程组得到的结果，你算出来=-1，肯定是你计算器没按对。

1.png

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

5楼2013-10-21 19:41:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 LESK乐的主题更新

返回列表