![](https://img-blog.csdnimg.cn/20190927151132530.png?x-oss-process=image/resize,m_fixed,h_224,w_224)
Python、数值分析、数学物理基础、matlab
文章平均质量分 68
python的各种技巧技能,以及经常出现的问题解决方法途径。
优惠券已抵扣
余额抵扣
还需支付
¥29.90
¥99.00
购买须知?
本专栏为图文内容,最终完结不会低于15篇文章。
订阅专栏,享有专栏所有文章阅读权限。
本专栏为虚拟商品,基于网络商品和虚拟商品的性质和特征,专栏一经购买无正当理由不予退款,不支持升级,敬请谅解。
nutron-ma
欢迎
展开
-
目录catalog
数值分析高斯消元法(python)解线性方程组雅可比迭代法(python)解线性方程组高斯-赛德尔迭代法(python)解线性方程组逐次超松弛迭代法SOR(python)解线性方程组原创 2021-01-16 21:21:22 · 305 阅读 · 1 评论 -
中子物理思考题-中子探测
1.什么是核反冲法和核反应法?核反冲法: 当快中子在核上散射时,便引起核反冲,再通过探测器记录反冲核。根据弹性散射时能量守恒和动量守恒定律核反应法:中子与核发生反应后放出带电粒子或γ射线,而且反应能足够大,在电离室和闪烁晶体中能产生 足够大的脉冲而被记录。2.简述3He正比计数管的原理及应用原理:利用核反应将3He充入电离室或正比计数管,根据电离脉冲大小可以测定反应产 总能量。应用:3.简述金刚石中子探测器的基本原理和优缺点缺点:CVD生长技术产生的多晶薄膜有相对较高的晶界密度和.原创 2021-06-17 13:42:56 · 413 阅读 · 0 评论 -
中子物理思考题-宏观相互作用、输运
1. 什么是宏观截面?给出平均自由程、平均迁移自由程和自由程的定义与区别。原创 2021-06-17 13:40:55 · 322 阅读 · 0 评论 -
中子物理思考题(第一章和第二章 )
1. 什么是中子星?原创 2021-06-16 15:26:37 · 152 阅读 · 0 评论 -
中子物理思考题-中子与核相互作用
中子物理思考题-中子与核相互作用原创 2021-05-30 23:45:53 · 140 阅读 · 0 评论 -
中子物理思考题-中子源
文章目录原创 2021-05-06 23:06:46 · 270 阅读 · 0 评论 -
梯度下降法&共轭梯度法
文章目录梯度下降法(gradient descent)共轭梯度法(conjugate gradient)梯度下降法(gradient descent)实现:系数矩阵使用希尔伯特(hilbert)矩阵H,假设解x=(1,1,1,···,1)T,则常数项为b = Hx。由此H和b建立方程Hx=b。'''梯度下降法(最速下降法)法求解线性方程组Ax=bmayuhua 2021.3.7 s说明:《数值分析》中说最速下降法是个不稳定的算法,误差会在传递中放大。所以先写的共轭梯度,这个就晚一天原创 2021-03-07 18:30:36 · 651 阅读 · 1 评论 -
twilio给自己发短信
使用python的twilio给自己发短信1.申请账号、虚拟电话号1.注册地址:https://www.twilio.comsign up, 人机验证,邮箱验证,验证接收短信的电话号码。获得自己的虚拟号码,获得账号2.pip install twilio# -*- coding: utf-8 -*-"""用于发送短信需要pip install twilio给自己发短信需要申请账号2020.8.29"""from twilio.rest import Client# Your Ac原创 2020-08-29 21:55:24 · 1694 阅读 · 1 评论 -
高斯消元法(python)解线性方程组
高斯消元法:1.消元:通过初等行变换一步步把系数矩阵化为一个上三角矩阵,2.回代:一个个解出xn···x2, x1一般性解法:设Ax=b,A \in R^{n\times n},若A的所有顺序主子式均不为零,则基本高斯消元无需换行进行到底,得到唯一解,其消元和回代的计算公式为:...原创 2020-12-30 12:15:26 · 1334 阅读 · 0 评论 -
雅可比迭代法(python)解线性方程组
雅可比迭代法,有两种迭代方式:数值计算迭代、矩阵方程迭代代码中的数值迭代实为gs迭代,方法二的矩阵迭代有点问题,有空修正。import numpy as npfrom numpy.linalg import *# 维数di = int(input('请输入方程的维数:'))h = np.zeros(shape=(di,di))h_shape = h.shape# 生成系数矩阵hfor i in range(h_shape[0]): for j in range(h_shape[1原创 2020-12-30 12:16:19 · 3738 阅读 · 1 评论 -
高斯-赛德尔迭代法(python)解线性方程组
import numpy as npfrom numpy.linalg import *# 维数di = int(input('请输入方程的维数:'))h = np.zeros(shape=(di,di))h_shape = h.shape# 生成系数矩阵hfor i in range(h_shape[0]): for j in range(h_shape[1]): h[i,j] = 1/(i+1+j+1 -1) # 索引是从0开始的,所以要多加上两个1prin原创 2020-12-30 12:17:26 · 1588 阅读 · 0 评论 -
逐次超松弛迭代法SOR(python)解线性方程组
import numpy as npfrom numpy.linalg import inv# 维数di = int(input('请输入方程的维数:'))h = np.zeros(shape=(di,di))h_shape = h.shape# 生成系数矩阵hfor i in range(h_shape[0]): for j in range(h_shape[1]): h[i,j] = 1/(i+1+j+1 -1) # 索引是从0开始的,所以要多加上两个1pr原创 2020-12-30 12:17:48 · 1373 阅读 · 0 评论 -
1.1 桌面说明
1.此处的路径说明我们现在所在的文件夹位置。接下来所编写的代码,数据等文件都会存放在这个路径下,在当前文件夹中可以看到现有的文件。2. 编辑器:在编辑器中可以写代码3. 命令行:可以写命令、代码,也用于显示编写代码的运行结果4. 工作区:可以看到代码中的变量、值...原创 2020-08-26 22:31:05 · 127 阅读 · 0 评论 -
1.2 数据和变量
%%1.2数据和变量 2020.8.25%分号;逗号,省略号(三个英文句点)...sum = 3+6-8+... %当程序太长时用...换行。记得别跟在数字后面,这里易被认为8点几 5-2%数据显示格式:format short, format long, format rational,用于显示的精度不同format longpiformat shortpi%复数value = 3 + 4*ivalue = (3+4*i)/5*i%预定义变量(ma原创 2020-08-26 22:41:55 · 129 阅读 · 0 评论 -
1.3 数组及其运算
%%数组及其运算 2020.8.26 马玉华a = [1,2,3;4,5,6;7,8,9]%冒号运算:[x1:公差:x2]生成一维数组b = [10:2:50]%linspace(x1,x2,n) 生成x1到x2间的n个数的一维数组linspace(1,8,8)%size() 求矩阵大小size(b)%取出矩阵中的值a(2,3) %第2行第3列的值a(2,[1,2,3]) %第2行第1,2,3列的值a(2,:) %整行:整个第2行的值原创 2020-08-26 22:50:20 · 165 阅读 · 0 评论 -
1.4 字符串操作
%%字符串char, 元组cell, 结构(字典)structure 2020.8.26 马玉华%字符串a = '12'b = double(a) %将字符串中的每一个字符都转换为ASCII码c = char(49,50)d = str2num('12') %字符串转数字。2(two)谐音toclass(d)e = num2str(12) %数字转字符串class(e)%字符串切片f = 'abcd'f(1)f(1:3)%元组和结构暂且不讲%stru原创 2020-08-26 22:55:45 · 121 阅读 · 0 评论 -
2.1.1 分支语句if、switch, 循环语句for、while
%%分支语句if、switch, 循环语句for、while, 其他pause、break、return、error 2020.8.26 马玉华%if分支结构%if 条件% 执行语句%elseif 条件% 执行语句%else% 执行语句%endx = input('请输入x的值:') %会自动识别double和charif x<=10; y = cos(x)elseif x == 10 y = sqrt(x)else y = si原创 2020-08-26 22:58:50 · 300 阅读 · 0 评论 -
2.1.2 M函数文件
%%函数文件基本结构 2020.8.26 %function [返回值] = 函数名(输入变量)% 函数体语句function [c,d] = calculate(a,b) c = cos(a) d = sqrt(b)%当代码写完后保存,自动生成一个和函数名相同名字的函数文件 %调用:在命令行输入函数名就可以调用。示例:% [value1,value2] = calculate(10,20) %m脚本文件中调用:方法同上,但是需要函数文件和脚本文原创 2020-08-26 23:01:19 · 111 阅读 · 0 评论 -
2.1.3 函数重载nargin与nargout
函数文件:%%函数重载nargin与nargout (number argument input 输入变量数目,number argument output)%函数文件如下:function [sum,mean] = sum_mean(a,b)if nargin == 1 % 当输入变量为1个时,执行 sum = a ; mean = a;elseif nargin == 2 % 当输入变量为2个时,执行 sum = a + b ;原创 2020-08-27 22:27:25 · 274 阅读 · 0 评论 -
2.1.4 提高计算速度,使用数组(一维、多维)计算,避免for循环一个个计算元素
如何提高运算速度?使用数组进行运算!避免for循环。示例:计算一个分段函数f(x)= x² , x>1; 1 , -1<x<=1;3+2x , x<=-1对比计算时间:①for循环计算、②矩阵进行计算、脚本代码:%%使用数组(一维、多维)计算,提高计算速度,避免for循环一个个计算元素。%计算分段函数(piecewise function)f(x)= x² , x>1% 1 ,原创 2020-08-27 22:49:10 · 389 阅读 · 0 评论 -
2.1.5句柄函数、匿名函数
% 函数句柄:% fun = @ Mfun 这里的Mfun是函数的M文件表达方式fun0 = @calculatefun0(1,3)% 匿名函数(类似lambda函数):% fun = @(x)表达式fun1 = @(x) x^2 + sin(x+1)fun1(1)calculate函数文件为:function [cos_,sqrt_] = calculate(a,b)cos_ = cos(a)sqrt_ = sqrt(b)end...原创 2020-09-12 11:03:20 · 392 阅读 · 2 评论 -
2.2.1曲线图
% 曲线figure(1) % 建立幕布x_arr = [0:0.1:20*pi];y_arr = exp(-0.1*x_arr).*cos(x_arr);plot(x_arr,y_arr)% 三维曲线figure(2)z_arr = x_arr .* sin(x_arr);plot3(x_arr,y_arr,z_arr)% 三维曲面figure(3)x = [-2:0.1:2];y = x;[x,y] = meshgrid(x,y) %网格化自变量z = x .*原创 2020-09-12 11:06:22 · 145 阅读 · 0 评论 -
3.1本征值和线性方程组
% % 本征值与本征向量a = rand(3) % 产生一个3阶 随机矩阵,(0,1)之间随机分布% eig(a) 方阵a的特征值[eigen_vector,eigenvalue] = eig(a) % 本征向量,本征值。(本征=特征)% % 线性方程组disp(' 求解ax=b。具体为:x+2y=1,3x-2y=4 ') % 唯一解a = [1,2;3,-2];b = [1;4];disp('rank(a),rank([a,b]),求出系数矩阵、增广矩阵的秩,如果秩相等有解,原创 2020-09-12 11:18:58 · 860 阅读 · 0 评论 -
4.1(非)线性函数求根、非线性方程组
线性函数求根、k次多项式拟合% 2020.9.7 %% 线性函数求根、k次多项式拟合% & 函数零点的指令% x = roots(p) 多项式p的根,实根复根都可以求% y=polyval(p,x) 求多项式p(x)的值,x可以写多个值% p = polyfit(x,y,k) 用k次多项式拟合(x,y),返回多项式的降幂系数[]% p = conv = (p1,p2) 返回多项式p1,p2的乘积% [quotient,remainder]=deconv(p1,p2) 求p1除以p2原创 2020-09-12 14:39:23 · 335 阅读 · 0 评论 -
4.3非线性函数最小二乘拟合
%% 非线性函数-最小二乘法拟合 马玉华 2020.9.11% 已知一组自变量数据x和函数值y,以及函数形式y=f(c,x),但是表达式的参数c不知道。对数据xy进行拟合,求解参数c% 基本格式 c=lsqcurvefit(fun,c0,x,y)% 完整格式 [c,Q]=lsqcurvefit(fun,c0,x,y,lb,ub),c0为需要设置的迭代初值,x、y为已知数据(矩阵,x为多元则使用多维数组),lb、ub为参数c的上下限(可不写)。c0、x、y、lb、ub都是向量% 返回拟合原创 2020-09-12 14:50:33 · 2474 阅读 · 4 评论 -
0. 入门:创建项目、编写脚本
matlab的2016a或2016b的字母ab代表上下半年的意思。打开matlab,可以在上面看到现在所处的路径,一般是软件安装时的路径。当我们开始编程时可以更改路径,把新文件夹的路径粘贴进来就可以了,新项目的程序文件就会保存在新文件夹中。然后新建一个脚本,就可以在脚本中编写程序啦!编辑器:可以在这儿编写程序。命令行窗口:可以在这儿做一些测试,像python中的console一样。命令行窗口也用于显示程序输出结果。...原创 2020-08-24 00:20:11 · 286 阅读 · 0 评论 -
1. 数据类型和操作
%行注释,代码后面加;可以不输出结果%变量命名规则和python相同%matlab的索引从1开始,python的索引从0开始%class()查看变量类型class('string_or_char') %字符串charclass(5) %数字double%输入输出n = input('请输入数字or字符串,不可混合:') %会自动识别double或char,不可混合输入class(n)n=5fprintf('this is an output \n'原创 2020-08-26 22:06:39 · 253 阅读 · 0 评论 -
2. for、while循环,if、switch分支
%符号;可加可不加,都不影响代码运行,只是阻止输出而已% for 循环格式:%for 条件% 执行语句%end%例子:求1?+2?+3?+4?+5?sum =0;for n = [1,2,3,4,5]; sum = sum + n^2;endsum% while循环:%while 条件% 执行语句%end%例子:求1+2+3……+10的和n = 1;sum = 0;while n <= 10; sum = sum + n; n原创 2020-08-26 22:11:07 · 133 阅读 · 0 评论 -
3. 绘制二维图像
%二维平面绘图x = [0:0.01:2*pi];y = sin(x);figure %建立一个幕布(有时也可以没有)plot(x,y)title('y=sin(x)')xlabel('x')ylabel('y')xlim([0,2*pi])%一个幕布两条线x = [0:0.01:2*pi];y1 = sin(x);y2 = 2*cos(2*x)+1;figure %建立一个幕布(有时也可以不写)plotyy(x,y1,x,y2)title('y=sin(原创 2020-08-26 22:11:57 · 340 阅读 · 0 评论 -
4. help、doc帮助与文档
%帮助与文档 (有时可能需要联网)%doc 函数 %可以查看函数的说明文档doc ifdoc switch%help 函数 %可以在命令行查看函数说明help if原创 2020-08-26 22:14:39 · 562 阅读 · 0 评论 -
希腊字母表、单位制
高斯制→国际制(cgs→mks)国际制→高斯制(mks→cgs)原创 2020-10-28 22:57:05 · 546 阅读 · 0 评论 -
分离变数法
典型问题例1:典型问题-第二类边界条件例2:原创 2020-05-03 12:51:30 · 1615 阅读 · 0 评论 -
常用微分方程
下面这个是已知一个解,求另一个解原创 2020-05-02 21:02:11 · 1175 阅读 · 0 评论 -
傅里叶级数
应用:1)傅里叶正弦展开原创 2020-05-02 20:52:08 · 173 阅读 · 0 评论 -
电子(自旋、轨道、耦合)磁矩
欢迎使用Markdown编辑器你好! 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章,了解一下Markdown的基本语法知识。新的改变我们对Markdown编辑器进行了一些功能拓展与语法支持,除了标准的Markdown编辑器功能,我们增加了如下几点新功能,帮助你用它写博客:全新的界面设计 ,将会带来全新的写...原创 2019-05-27 21:48:50 · 3737 阅读 · 1 评论