黑松山资源网 Design By www.paidiu.com
使用scipy.optimize模块的root和fsolve函数进行数值求解线性及非线性方程,下面直接贴上代码,代码很简单
from scipy.integrate import odeint import numpy as np import matplotlib.pyplot as plt from scipy.optimize import root,fsolve #plt.rc('text', usetex=True) #使用latex ## 使用scipy.optimize模块的root和fsolve函数进行数值求解方程 ## 1、求解f(x)=2*sin(x)-x+1 rangex1 = np.linspace(-2,8) rangey1_1,rangey1_2 = 2*np.sin(rangex1),rangex1-1 plt.figure(1) plt.plot(rangex1,rangey1_1,'r',rangex1,rangey1_2,'b--') plt.title('$2sin(x)$ and $x-1$') def f1(x): return np.sin(x)*2-x+1 sol1_root = root(f1,[2]) sol1_fsolve = fsolve(f1,[2]) plt.scatter(sol1_fsolve,2*np.sin(sol1_fsolve),linewidths=9) plt.show() ## 2、求解线性方程组{3X1+2X2=3;X1-2X2=5} def f2(x): return np.array([3*x[0]+2*x[1]-3,x[0]-2*x[1]-5]) sol2_root = root(f2,[0,0]) sol2_fsolve = fsolve(f2,[0,0]) print(sol2_fsolve) # [2. -1.5] a = np.array([[3,2],[1,-2]]) b = np.array([3,5]) x = np.linalg.solve(a,b) print(x) # [2. -1.5] ## 3、求解非线性方程组 def f3(x): return np.array([2*x[0]**2+3*x[1]-3*x[2]**3-7, x[0]+4*x[1]**2+8*x[2]-10, x[0]-2*x[1]**3-2*x[2]**2+1]) sol3_root = root(f3,[0,0,0]) sol3_fsolve = fsolve(f3,[0,0,0]) print(sol3_fsolve) ## 4、非线性方程 def f4(x): return np.array(np.sin(2*x-np.pi)*np.exp(-x/5)-np.sin(x)) init_guess =np.array([[0],[3],[6],[9]]) sol4_root = root(f4,init_guess) sol4_fsolve = fsolve(f4,init_guess) print(sol4_fsolve) t = np.linspace(-2,12,2000) y1 = np.sin(2*t-np.pi)*np.exp(-t/5) y2 = np.sin(t) plt.figure(2) a , = plt.plot(t,y1,label='$sin(2x-\pi)e^{-x/5}$') b , = plt.plot(t,y2,label='$sin(x)$') plt.scatter(sol4_fsolve,np.sin(sol4_fsolve),linewidths=8) plt.title('$sin(2x-\pi)e^{-x/5}$ and $sin(x)$') plt.legend()
以上这篇python scipy求解非线性方程的方法(fsolve/root)就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
黑松山资源网 Design By www.paidiu.com
广告合作:本站广告合作请联系QQ:858582 申请时备注:广告合作(否则不回)
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
黑松山资源网 Design By www.paidiu.com
暂无评论...
稳了!魔兽国服回归的3条重磅消息!官宣时间再确认!
昨天有一位朋友在大神群里分享,自己亚服账号被封号之后居然弹出了国服的封号信息对话框。
这里面让他访问的是一个国服的战网网址,com.cn和后面的zh都非常明白地表明这就是国服战网。
而他在复制这个网址并且进行登录之后,确实是网易的网址,也就是我们熟悉的停服之后国服发布的暴雪游戏产品运营到期开放退款的说明。这是一件比较奇怪的事情,因为以前都没有出现这样的情况,现在突然提示跳转到国服战网的网址,是不是说明了简体中文客户端已经开始进行更新了呢?
更新日志
2024年10月08日
2024年10月08日
- 《丁当 20首重量级歌曲 Fu Good 下一站天后 2CD》[WAV/分轨][650MB]
- 歌莉雅.2012-My.Voice.Story【环星】【FLAC分轨】
- 群星.1993-一曲成名·青春无悔【飞碟】【WAV+CUE】
- 刘德华.2000-爱无知(金装版)【NEWMELODY】【WAV+CUE】
- 许巍《时光漫步》1CD[FLAC/分轨][804.3MB]
- 群星《2024第一季度百度歌曲排行抖音榜top100》1CD[FLAC/分轨][2.3GB]
- 张柏芝《全新经验》[FLAC/分轨][450MB]
- 唐朝.2008-浪漫骑士【音乐家】【WAV+CUE】
- 张芸京.2016-失败的高歌【泡耳音乐】【WAV+CUE】
- 群星.1991-音乐工厂1·皇后大道东【音乐工厂】【WAV+CUE】
- 群星.1992-音乐工厂2·首都【音乐工厂】【WAV+CUE】
- 群星.1994-音乐工厂3·儿童乐园【音乐工厂】【WAV+CUE】
- 胡鸿钧.2021-ex:CHANGE【星梦娱乐】【WAV+CUE】
- 班得瑞原装进口《第一张新世纪专辑:仙境》1CD[APE/CUE分轨][292.3MB]
- 伍佰《摇滚教父 伍佰黄金精选 2CD》[WAV+CUE][990MB]