MaYuHuaw的博客

数值分析高斯消元法（python）解线性方程组雅可比迭代法（python）解线性方程组高斯-赛德尔迭代法（python）解线性方程组逐次超松弛迭代法SOR（python）解线性方程组

1.什么是核反冲法和核反应法？核反冲法： 当快中子在核上散射时，便引起核反冲，再通过探测器记录反冲核。根据弹性散射时能量守恒和动量守恒定律核反应法：中子与核发生反应后放出带电粒子或γ射线，而且反应能足够大，在电离室和闪烁晶体中能产生 足够大的脉冲而被记录。2.简述3He正比计数管的原理及应用原理：利用核反应将3He充入电离室或正比计数管，根据电离脉冲大小可以测定反应产 总能量。应用：3.简述金刚石中子探测器的基本原理和优缺点缺点：CVD生长技术产生的多晶薄膜有相对较高的晶界密度和.

1. 什么是宏观截面？给出平均自由程、平均迁移自由程和自由程的定义与区别。

1. 什么是中子星？

中子物理思考题-中子与核相互作用

文章目录

文章目录梯度下降法（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说明：《数值分析》中说最速下降法是个不稳定的算法，误差会在传递中放大。所以先写的共轭梯度，这个就晚一天

使用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

高斯消元法：1.消元：通过初等行变换一步步把系数矩阵化为一个上三角矩阵，2.回代：一个个解出xn···x2, x1一般性解法：设Ax=b，A \in R^{n\times n},若A的所有顺序主子式均不为零，则基本高斯消元无需换行进行到底，得到唯一解，其消元和回代的计算公式为：...

雅可比迭代法，有两种迭代方式：数值计算迭代、矩阵方程迭代代码中的数值迭代实为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

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

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

1.此处的路径说明我们现在所在的文件夹位置。接下来所编写的代码，数据等文件都会存放在这个路径下，在当前文件夹中可以看到现有的文件。2. 编辑器：在编辑器中可以写代码3. 命令行：可以写命令、代码，也用于显示编写代码的运行结果4. 工作区：可以看到代码中的变量、值...

%%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.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行的值

%%字符串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

%%分支语句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.8.26 %function [返回值] = 函数名（输入变量）% 函数体语句function [c,d] = calculate(a,b) c = cos(a) d = sqrt(b)%当代码写完后保存，自动生成一个和函数名相同名字的函数文件 %调用：在命令行输入函数名就可以调用。示例：% [value1,value2] = calculate(10,20) %m脚本文件中调用：方法同上，但是需要函数文件和脚本文

函数文件：%%函数重载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 ;

如何提高运算速度？使用数组进行运算！避免for循环。示例：计算一个分段函数f(x)= x² , x>1； 1 , -1<x<=1；3+2x , x<=-1对比计算时间：①for循环计算、②矩阵进行计算、脚本代码：%%使用数组(一维、多维)计算，提高计算速度，避免for循环一个个计算元素。%计算分段函数（piecewise function）f(x)= x² , x>1% 1 ,

% 函数句柄：% 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...

% 曲线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 .*

% % 本征值与本征向量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]),求出系数矩阵、增广矩阵的秩，如果秩相等有解，

线性函数求根、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.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都是向量% 返回拟合

matlab的2016a或2016b的字母ab代表上下半年的意思。打开matlab，可以在上面看到现在所处的路径，一般是软件安装时的路径。当我们开始编程时可以更改路径，把新文件夹的路径粘贴进来就可以了，新项目的程序文件就会保存在新文件夹中。然后新建一个脚本，就可以在脚本中编写程序啦！编辑器：可以在这儿编写程序。命令行窗口：可以在这儿做一些测试，像python中的console一样。命令行窗口也用于显示程序输出结果。...

%行注释，代码后面加;可以不输出结果%变量命名规则和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'

%符号;可加可不加，都不影响代码运行，只是阻止输出而已% 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

%二维平面绘图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(

%帮助与文档 (有时可能需要联网)%doc 函数 %可以查看函数的说明文档doc ifdoc switch%help 函数 %可以在命令行查看函数说明help if

高斯制→国际制（cgs→mks）国际制→高斯制(mks→cgs)

典型问题例1:典型问题-第二类边界条件例2:

下面这个是已知一个解，求另一个解

应用：1）傅里叶正弦展开

欢迎使用Markdown编辑器你好！ 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章，了解一下Markdown的基本语法知识。新的改变我们对Markdown编辑器进行了一些功能拓展与语法支持，除了标准的Markdown编辑器功能，我们增加了如下几点新功能，帮助你用它写博客：全新的界面设计 ，将会带来全新的写...