自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+

小憨的学习博客

数学与计算科学专业的学生一枚,在努力的输出~

  • 博客(65)
  • 资源 (21)
  • 收藏
  • 关注

原创 jupyter-notebook使用指南

记录一些自己使用过程中的问题,以免以后忘记。

2023-06-27 12:06:54 1845

原创 关于clash退出后,华硕电脑连不上网了

今天早上突然浏览器的网页(微软自带、华硕自带,谷歌浏览器)都无法打开,并且手机热点我试了,重启也试了,只是网页打不开,其他软件联网正常。——问题已解决

2023-05-07 09:48:32 15526 16

原创 Python中fealpy库的安装

记录一下老师上课的内容以及操作

2022-05-07 11:58:25 1719 2

原创 Python包pretty_errors

Python包pretty_errors1. 配置打开cmd或者anaconda prompt,输入下面命令。python -m pip install pretty_errors可以通过import模块包来用这个东西,为了避免麻烦,可以直接加到环境里,就不用每次都导入了,用下面这个语句。python -m pretty_errors2.测试结果【VS-Code】x=-100print(-x+50)print(1/0)这个东西还蛮好用吧,简洁明了!...

2022-04-08 16:51:04 1164

原创 用U盘给电脑重装系统记录

电脑重装系统记录1.插入系统U盘,进入电脑的bios界面2.点击USB选项3.进入U盘系统,打开冰封U启动界面,点击win10x64.esd,再次点击win10专业版4.格式化本地磁盘,等待安装,然后自动重启

2022-01-04 19:01:19 473

原创 LaTeX问题汇总

LaTeX问题汇总\linebreak[n]作用\linebreak[n]作用n=0,1,2,3,4

2021-08-22 10:53:49 121

原创 seaborn库——分类图

分类图分类散点图的绘制分布图箱线图小提琴图分类统计图条形图柱状图点线图分类散点图的绘制sns.catplot函数sns.swarmplot函数数据量比较大的时候不适合使用这个函数,因为不美观可以用第一个函数加一个属性 kind="swarm"来生成这个效果swramplot采用了一定的算法,可以让点不会重叠。分布图箱线图sns.boxplot函数箱线图:boxplot或者是catplot(kind=“box”)使用hue参数可以进行分组。小提琴图sns.vio

2021-08-11 11:04:39 364

原创 python报错及解决办法汇总

python报错及解决办法汇总matplotlib'gbk' codec can't decode byte 0x9a in position 37: illegal multibyte sequencematplotlib‘gbk’ codec can’t decode byte 0x9a in position 37: illegal multibyte sequence解决办法:用VS-code直接写这个文件,剪切到配置文件路径,然后重新运行就不报错了猜测是因为txt文件的编码是gbk

2021-08-11 10:19:48 371 1

原创 lingo进阶篇

lingo进阶篇运算符关系运算符逻辑运算符lingo内置函数变量定界函数数学函数集合操作函数运算符x = 2;y = 3*x^10 + 6/(15-x^(1/2));关系运算符B = 10;e = 0.0001;A - e > B;逻辑运算符model:sets:factory / 1..6 / : a;endsetsdata:a = 6,5,4,3,2,1;enddatay = @sum( factory(i) | i#ge#5 : a(i) );

2021-07-26 20:51:01 3565

原创 lingo入门-集合和矩阵

lingo入门-集合和矩阵例题引入矩阵工厂例题引入解法【暴力枚举法】矩阵工厂测试代码sets:factory /1..6/ : a,b;plant /1..3/ : x,y;endsets① factory 和 plant 都是制造矩阵的工厂,但它们是两家不同的工厂。② factory 工厂后面的 /1…6/ 说明它专门生产 1×61\times61×6 的矩阵。factory 工厂最后面出现的 a 和 b,都是 的矩阵。③ plant 工厂后面的 /1…3/ 说明它专门生

2021-07-25 22:18:22 13757

原创 Lingo基本教程

Lingo基本教程lingo基本操作用lingo解方程lingo变量线性规划问题lingo基本操作编译代码和运行lingo强大之处求得运行结果用lingo解方程实例:易错点:不写结尾的分号; 不写乘法符号*lingo变量实例:x>0是lingo中默认的条件,不用写线性规划问题实例:...

2021-07-25 21:56:07 2105

原创 STATA基础教程

STATA基础教程stata基本命令stata基本命令1.使用命令use "文件路径.dta"例如: use “H:\resourses\2021数学建模培训\多元回归案例分析\数据.dta”2.通过导航栏【文件-打开】打开所需数据下图演示如何打开及查看数据表3. 查看数据表信息:describe(简写des)4. 把数据表列出来:list5. 概要统计命令:summarize(简写sum)利用这个命令,可以得到:观测案例数(Obs)、变量的平均值(Mean)、标准差、最小值和最

2021-07-25 13:43:33 9080 1

原创 罚函数——内点罚函数算法

罚函数——内点罚函数算法内点罚函数算法算法算例matlab程序测试文件输出结果内点罚函数算法算法算例matlab程序function [x,val]=neidian(x0,epsilon)m(1)=1;c=10;% c为递增系数for k=1:100 e=m(k); x0=bfgs('ndfz','dndfz',x0); f1=zzzz(x0,e); if(norm(f1)<=epsilon) break else

2021-07-09 17:56:39 3114 2

原创 罚函数法——外点罚函数法

罚函数法——外点罚函数法外点罚函数法算法算例matlab代码测试文件输出结果外点罚函数法算法算例matlab代码function [x,val]=waidian1(x0,epsilon)m(1)=1;c=10;% c为递增系数for k=1:100 e=m(k); x0=bfgs('zg','dzg',x0); f=zg(x0,e); if(f<=epsilon) break else m(k+1)=c*m(k

2021-07-09 17:25:46 3550 3

原创 罚函数法——乘子法

罚函数法——乘子法乘子法乘子法例题代码function [x,mu,lam,output]=multphr(fun,hf,gf,dfun,dhf,dgf,x0)maxk=1000;sigma=2.0;theta=0.8; eta=2.0; k=0; ink=0;epsilon=1e-5;x=x0; he=feval(hf,x); gi=feval(gf,x);n=length(x); l=length(he); m=length(gi);mu=0.1*ones(l,1)

2021-07-09 14:08:44 974 3

原创 最优化方法——外点罚函数法

外点罚函数法还没搞定…function [A,B]=Exterior(f,p,eps)format long;d=1;e1=1e-4;k=1;A(k)=0;B(k)=0;r(k)=1;a=2;while d>eps x1=A(k);x2=B(k); if vpa(p(x1,x2))>0 u1=1; else u1=0; end F=f+r(k)*u1*p^2; Fx1=diff(F,'x1');

2021-06-17 11:56:02 1381

原创 最优化方法——最小二乘法

最小二乘法No.1No.2No.3No.4No.1function[x,res]=nels(A,b)B=A'*A;f=A'*b;L=chol(B,'lower');y=L\f;x=L'\y;res=norm(b-A*x);测试文件clear;clc;A=[2 3 4 5; 4 3 2 1;4 5 6 7;9 5 7 2;4 2 5 3];b=[20 22 35 42 50]';[x,res]=nels(A,b)运行结果No.2function [x,res]=qrls

2021-06-10 17:16:52 935 1

原创 最优化方法——信頼域算法

信頼域算法信頼域方法牛顿型信頼域方法信頼域方法信頼域方法程序function[d,val,lam,i]=trustq(fk,gk,Bk,deltak)n=length(gk);beta=0.6;sigma=0.2;mu0=0.05;lam0=0.05;gamma=0.05;d0=ones(n,1);z0=[mu0,lam0,d0']';zbar=[mu0,zeros(1,n+1)]';i=0;z=z0;mu=mu0;lam=lam0;d=d0;while(i<=15

2021-06-10 16:57:15 215

原创 最优化方法——拟牛顿法

拟牛顿法对称秩1算法BFGS算法DPF算法Broyden族算法对称秩1算法目标函数function f=fun(x)f=2*(x(1)-x(2)^2)^2+(x(2)-2)^2;梯度function gf=gfun(x)gf=[4*(x(1)-x(2)^2);-8*x(2)*(x(1)-x(2)^2)+2*(x(2)-2)];对称秩1算法程序function[k,x,val]=sr1(fun,gfun,x0,epsilon,N)if nargin<5 N=1000;e

2021-06-10 11:51:32 939

原创 最优化方法——共轭梯度法

共轭梯度法线性共轭梯度法FR非线性共轭梯度法线性共轭梯度法线性共轭梯度法程序function[k,x,val]=linecg(A,b,x0,epsilon,N)if nargin<5 N=1000;endif nargin<4 epsilon=1.e-5;endif nargin<3 x0=zeros(length(b),1);endk=0;gk=A*x0-b;dk=-gk;while(k<N) temp=A*dk;

2021-06-10 11:19:30 810

原创 最优化方法——梯度法与牛顿法

梯度法与牛顿法梯度法阻尼牛顿法梯度法目标函数function f=fun(x)f=4*(x(1)^2-x(2))^2+3*(x(1)-1)^2;梯度function gf=gfun(x)gf=[16*x(1)*(x(1)^2-x(2))+6*(x(1)-1); -8*(x(1)^2-x(2))];梯度法的MATLAB实现function[k,x,val]=grad(fun,gfun,x0,epsilon)maxk=5000;beta=0.5;sigma=0.4;k=0;whi

2021-06-03 11:55:19 594 2

原创 最优化方法——线搜索

黄金分割法MATLAB实现MATLAB实现待求解函数function [y]=phi(x)y=3*x^2-2*tan(x);黄金分割法MATLAB实现function [i,s,phis,ds,dphi,G]=golds(phi,a,b,epsilon,delta) % Golden Sectiont=(sqrt(5)-1)/2; h=b-a; %Define interval lengthphia=feval(phi,a); % Write function value

2021-06-03 10:55:03 694

原创 数值计算方法——不动点迭代和牛顿迭代法

不动点迭代和牛顿迭代法不动点迭代feval函数简单迭代法不动点迭代不动点可看成 y=φ(x)与y=xy=φ(x)与y=xy=φ(x)与y=x 的交点feval函数功能是求函数值基本使用格式:y=feval(fhandle, x) %fhandle——函数表达式,x——变量值[y1, y2, ...] = feval(fhandle, x1,..., xn)简单迭代法function [x_reality,n_reality] = Simple_stepit( f_name,x_

2021-06-01 15:24:23 4739

原创 MATLAB——使用问题及解决措施

使用问题及解决措施报错没有为 'function_handle' 输入定义函数报错没有为 ‘function_handle’ 输入定义函数函数需要单独用M文件定义。如下举例:先定义函数fun_ex5.m文件 function f=fun_ex5(x) f=(x^4)+2*x^2-x-3;下面用简单迭代法求解该函数,则用’此处文件’x0=1;tolerance=1e-6;n_limit=100;[x_reality,n_reality] = Simple_stepit( 'fun_e

2021-05-31 11:51:20 615

原创 LaTeX——问题及解决办法

LaTeX——代码记录 $$ \begin{aligned} R_{n}(G)&=E_{z} E_{\sigma}\left[\left.\sup _{u \in K_{1}(p)}\left|\frac{1}{n} \sum_{j=1}^{n} \sigma_{j} \frac{1}{2}\right| \nabla u\right|^{2} \mid z_{i}, z_{n}\right]\\ &=\frac{1}{2} E_{z} E \sigma\left[\left.\

2021-05-31 10:20:43 7969 5

原创 数学建模复习——代码记录

数学建模复习线性规划问题的MATLAB求解线性规划问题的MATLAB求解linprog函数主要用来求线型规划中的最小值问题(最大值的镜像问题,求最大值只需要加个“-”)函数linprog格式:[x,fval,exitflag]=linprog(f,A,b,Aeq,beq,lb,ub,x0,options)f:目标函数A:不等式约束条件矩阵b:对应不等式右侧的矩阵Aeq:等式约束条件矩阵beq:不等式右侧的矩阵Aeq:等式约束条件矩阵beq:对应等式右侧的矩阵lb:x的下界ub:

2021-05-25 16:20:44 1241

原创 Python【小】问题记录

Python【小】问题记录matplotlib相关问题matplotlib相关问题from matplotlib import pyplot as plt x=range(2,26,2)y=[15,13,14.5,2,3,6,11,23,24,5,20,9]#数据的输入是用逗号隔开的!plt.plot(x, y)plt.show()#图形窗口太小了有些问题解答见另一篇文章Python——matplotlib库的学习入门如何设置图形大小?如何保存到本地?如何添加描述信息?plt.

2021-05-25 15:00:00 35

原创 Python——matplotlib库的学习入门

matplotlib库的学习入门安装依赖库NumPy Matplotlib图形中文显示[Matplotlib 教程](https://www.runoob.com/w3cnote/matplotlib-tutorial.html)格式化图形其他绘图函数subplot()bar()numpy.histogram()plt()安装依赖库CMD中输入下面语句pip install matplotlib或者下面这句pip3 install matplotlib -i https://pypi.tuna.

2021-05-24 19:09:25 328 2

原创 数值计算方法——乘幂法和反乘幂法

乘幂法和反乘幂法一级目录二级目录三级目录乘幂法一级目录二级目录三级目录乘幂法function [lamda,x,iter]=chengmi_solve(z,A,epsilon,max_iter)% n=length(A);% z=zeros(n:n);error=1;lamda=0;iter=0;while error>epsilon && iter<max_iter y=A*z; m=max(abs(y))

2021-05-24 11:54:43 1338

原创 视频下载

视频下载使用you-get批量下载视频使用you-get批量下载视频打开CMD pips install you-get安装完后你可以切换目录,比如你想下载到F盘,就F:,CD切一下目录,然后就可以执行下载命令了:you-get url先跳到你想下载的文件夹下于是可以拥有了!如果格式不对,需要一个格式转换使用这个软件 FFmpeg要用到其中的bin文件没有及时记录,未完待续。。。...

2021-05-20 21:33:20 493

原创 C语言习题——统计字符串

统计字符串题目代码题目代码#include <stdio.h>#include<string.h>int main(){ char str[10]; int i,n,num=0,space=0,cha=0,qita=0; printf( "Enter a value :"); gets( str ); n=strlen(str); for(i=0;i<n;i++){ printf( "%d\n",str[i]);

2021-05-19 13:55:23 181

原创 数值计算方法——Jacobi迭代法和G-S迭代法

Jacobi迭代法和G-S迭代法一级目录二级目录矩阵常用操作一级目录二级目录矩阵常用操作A = [1,-1,2,3;5,1,-4,2;3,0,5,2;11,15,0,9];B = [0.43,43,2;-8.9,4,21];%(1)方阵diag(A) %主对角元素triu(A) %上三角tril(A) %下三角inv(A) %矩阵的逆det(A) %行列式的值rank(A) %矩阵的秩norm(A) %矩阵的范数cond(A) %矩阵的条件数trace(A)%矩阵的迹%(

2021-05-17 11:51:33 2630

原创 数值计算方法——Gauss消去实现

Gauss消去实现function A=GaussXQ(A,b)n=length(b);L=zeros(n,n);C=[A b];D=zeros(n,n+1);for j=1:n-1 %执行n-1次,消去成为上三角 if A(j-1,j-1)==0 for k=j:n if A(k,k)~=0 C([j-1,k],;)=C([k,j-1],;); break

2021-05-10 11:52:20 226

原创 anoconda——包管理和环境管理

包管理和环境管理conda命令——安装包conda命令pip命令环境管理包的安装、卸载、更新和查找等等conda命令——安装包使用conda安装时能自动安装相应的依赖包#conda install [包名称列表]conda install scipyconda命令conda config --add channels https://mirrors.tuna.tsinghua.edu.cn/anaconda/pkgs/free/pip命令#安装包pip install numpy

2021-05-03 21:55:35 107

原创 JupyterNotebook的使用

JupyterNotebook的使用JupyterNotebook介绍配置文件夹页面介绍编辑代码![在这里插入图片描述](https://img-blog.csdnimg.cn/20210503212111890.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQ3OTI1ODM2,size_16,color_FFFFFF,t_70)Jup

2021-05-03 21:27:56 1422

原创 计算机网络实验7——网页爬取

Python网页爬取环境实验代码三级目录环境anaconda实验代码import requests #导入库keyword = '数据分析' #设定搜索关键词是python,此处可以修改为任意你想搜索的关键词try: kv = {'q':keyword} r = requests.get("http://www.so.com/s",params=kv)#使用某个浏览器及其网址格式 print(r.request.url) #返回访问的网址except:

2021-04-29 17:55:47 325

原创 数值计算方法——数值积分与数值微分

数值积分与数值微分的MATLAB实现数值积分梯形公式中矩形公式Simpson公式三级目录数值积分梯形公式function f=f(a)%积分函数f=a^2;function s=tx(a,b)%梯形公式fa=f(a);fb=f(b);s=(b-a)/2*(fa+fb);function t=fhtx(a,b,n)xx=linspace(a,b,n+1);s=zeros(1,n);t=0;for i=1:n a=xx(i); b=xx(i+1); s(

2021-04-26 11:16:38 826

原创 数学建模——多元统计分析例题及程序

多元统计分析例题及程序主成分分析简述基本思想计算步骤例题程序因子分析相关性分析回归分析聚类分析主成分分析简述主成分分析(Principal Component Analysis,PCA), 是一种数学降维的统计方法。通过正交变换将一组可能存在相关性的变量转换为一组线性不相关的变量,转换后的这组变量叫主成分。在用统计分析方法研究多变量的课题时,变量个数太多就会增加课题的复杂性。人们自然希望变量个数较少而得到的信息较多。主成分分析作为基础的数学分析方法,其实际应用十分广泛,比如人口统计学、数量地理

2021-04-22 11:38:38 5709

原创 数值计算方法——插值公式的实现

插值公式的实现lagrange插值公式函数Hemite插值公式三级目录lagrange插值公式函数function Ln=lagrange(x,y,w)clcclearx=input('请输入x向量:');y=input('请输入y向量:');w=input('请输入待求值:');m=length(x);n=length(y);if m~=n error('错误,请输入正确维度')endLn=0;for i=1:n z=1;%基函数,每次循环前置为1,这个很重要

2021-04-19 11:52:24 224

原创 LaTeX公式

代码暂存\usepackage{package}\usepackage{amsmath}\usepackage{mathdots}\usepackage{microtype}\usepackage{geometry}\begin{document}\begin{spacing}{2.0}\section{Class F,Consider conv(F)}\\Lemma:\\R(conv(F))=R(F)\\\begin{aligned} & \mathbb

2021-04-16 02:01:29 206

实变函数复习重点.pdf

数学与计算科学——一篇自己整理的适用于实变函数专业期末考试复习的资料

2021-07-03

非线性方程的数值解法.pdf

信息与计算科学专业——一篇适用于数值计算方法期末考试的复习笔记

2021-07-03

函数基本逼近(二)——最佳逼近.pdf

信息与计算科学专业——一篇适用于数值计算方法期末考试的复习笔记

2021-07-03

函数基本逼近(一)——插值逼近.pdf

信息与计算科学专业——一篇适用于数值计算方法期末考试的复习笔记

2021-07-03

矩阵特征值问题的解法.pdf

信息与计算科学专业——一篇适用于数值计算方法期末考试的复习笔记

2021-07-03

数值积分与数值微分.pdf

信息与计算科学专业——一篇适用于数值计算方法期末考试的复习笔记

2021-07-03

线性代数方程组的求解.pdf

信息与计算科学专业——一篇适用于数值计算方法期末考试的复习笔记

2021-07-03

数值计算方法——误差分析.pdf

信息与计算科学专业——一篇适用于数值计算方法期末考试的复习笔记

2021-07-03

概率论——n元随机变量.pdf

信息与计算科学专业——概率论考试复习自己整理的资料

2021-07-03

概率论——抽样分布.pdf

信息与计算科学专业——概率论考试复习自己整理的资料

2021-07-03

概率论——二元分布.pdf

信息与计算科学专业——概率论考试复习自己整理的资料

2021-07-03

概率论——二元随机变量.pdf

信息与计算科学专业——概率论考试复习自己整理的资料

2021-07-03

概率论——概率相关.pdf

信息与计算科学专业——概率论考试复习自己整理的资料

2021-07-03

概率论——估计.pdf

信息与计算科学专业——概率论考试复习自己整理的资料

2021-07-03

概率论——基本概型.pdf

信息与计算科学专业——概率论考试复习自己整理的资料

2021-07-03

概率论——随机变量的数字特征.pdf

信息与计算科学专业——概率论考试复习自己整理的资料

2021-07-03

概率论——一元分布.pdf

信息与计算科学专业——概率论考试复习自己整理的资料

2021-07-03

复变函数论——解析函数的幂级数表示法.pdf

信息与计算科学专业——复变函数考试复习自己整理的资料

2021-07-03

复变函数论——复数与复变函数.pdf

信息与计算科学专业——复变函数考试复习自己整理的资料

2021-07-03

复变函数论——复变函数的积分.pdf

信息与计算科学专业——复变函数考试复习自己整理的资料

2021-07-03

复变函数论——解析函数.pdf

信息与计算科学——复变函数考试复习自己整理的资料

2021-07-03

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除