sqyao_147-CSDN博客

原创 MATLAB----级数

本篇参考于：中国大学慕课，专题六，“7.3级数"1.级数求和2.泰勒级数sum也可以进行级数求和，但是项数必须有限，并且要表示成向量形式，如果级数是无限的且项数非常多，用sum无法解决1.级数求和调用格式：symsum(s,v,n,m)其中，s表示一个级数的通项，是一个符号表达式。v是求和变量，v省略时使用系统的默认变量。n和m是求和变量v的初值和末值。在级数计算中，因为小数都表示有理分数的形式，容易导致分子或分母出现极大整数从而无法计算的情况，出现了具体的级数，..

2021-08-26 14:55:39 5731

原创 MATLAB----符号微积分

本篇参考于：中国大学慕课，专题六，“7.2符号微积分"1.符号函数的极限limit2.符号函数的导数diff3.符号函数的积分int1.符号函数的极限调用格式：limit(f,x,a)f为求极限的函数，x为变量，a为常数。即求函数f关于变量x在a点处的极限。a和x都可以省略，若x省略，则x默认为系统自变量，若a省略，a的默认值是0。如果f的左右极限不等，还可以求单边极限，调用格式为limit(f,x,a,'right')，limit(f,x,a,'left')2...

2021-08-25 15:48:15 1547

原创 MATLAB---非线性方程求解与函数极值运算

本篇参考于：中国大学慕课，专题六，“6.4MATLAB---非线性方程求解与应用”1.单变量非线性方程求解fzero2.非线性方程组求解fsolve3.求函数极值无约束条件最优化求解fminbnd,fminsearch,fiminunc有约束条件最优化求解fmincon1.单变量非线性方程求解注：fzero执行的是一个数值搜索的过程，搜索结果依赖于函数特性和指定的初始值。2.非线性方程组求解fzero使用fssolve求解时，赋初值0.1...

2021-08-24 22:26:06 1827

原创 MATLAB----线性方程组求解

本篇参考于：中国大学慕课，专题六，“6.2线性方程组求解”1.直接法0高斯消去法1.1列主元消去法（计算机标准算法）1.2矩阵三角分解法[L,U]=lu(A)/(L,U,P)=lu(A)2.迭代法2.1雅克比迭代法2.21.直接法1.1 列主元消去法（左除）注意：如果矩阵A是奇异的或接近奇异的，则MATLAB会给出警告信息A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]b=([13,-9,6,0])'x=A\b.

2021-08-23 18:24:06 1218

原创 MATLAB----数值微分与数值积分

本篇参考于：中国大学慕课，专题六，“6.1数值微分与数值积分”1.数值微分的实现diff2.数值积分的实现quad&quadl（quad精度高）integral （上下限可无穷）quadgk （上下限可无穷，适用于震荡函数）trapz （没有具体函数，只有样本数据）3.三重积分数学上没有或很难计算微分的计算解决1.数值微分的实现注意：diff函数计算的是向量元素之间的差分，故差分向量元素的个数比原向量少了一个。同样，对于矩阵来说，差分后的矩阵.

2021-08-23 12:49:07 795

原创 MATLAB----数据拟合

本篇参考于：中国大学慕课，专题五，“5.5数据拟合”曲线拟合常用来预测，或者根据散点拟合函数曲线（与插值不同的是并非所有点都经过该条曲线）多项式你也的函数为polyfit()，其功能为求得最小二乘拟合多项式系数，其调用格式为1）P=polyfit(X,Y,m)2）[P,S]=polyfit(X,Y,m)3）[P,S,mu]=polyfit(X,Y,m)根据样本数据X和Y，产生一个m次多项式P及其在采样点误差数据S，mu是一个二元向量，mu(1)是mean(X)，而mu(2)是std(

2021-08-22 10:42:58 8383

原创 MATLAB----数据插值

本篇参考于：中国大学慕课，专题五，“5.3数据插值”1.一维数据插值interp1()2.二维数据插值interp2()1.一维数据插值调用格式：interp1(X,Y,X1.method)。根据X，Y的值，计算函数在X1处的值。其中X,Y为已知的等长向量，分别表示采样点和采样点的值，X1表示一个标量或向量，表示要插值的点。method插值方法：linear线性插值，默认插值方法nearest最近点插值，插值点为后一个相邻点的函数值pchip分段三次埃尔米特插值。插值较为

2021-08-21 10:53:09 4029

原创 MATLAB----多项式计算

本篇参考于：中国大学慕课科学计算与MATLAB语言，专题五“5.2多项式计算”1.多项式的表示2.多项式的四则运算3.多项式的求导4.多项式的求值5.多项式的求方阵根1.多项式的表示多项式的表示用向量表示，即对于一个n次多项式用向量[an,an-1,an-2....a1,a0]表示。创建多项式向量时，要注意三点：1）多项式系数向量的顺序是从高到低。2）多项式系数向量包含0次项系数，所以其长度为多项式的最高次数加1。3）如果有的项没有，则系数向量相应位置用0补。2.多项式的

2021-08-20 15:31:57 1318

原创 MATLAB----数据统计分析

本篇参考于：中国大学慕课科学计算与MATLAB语言，专题五“5.1数据1.统计分析”1.求最大元素与最小元素max,min2.求平均值与中值mean,median3.求和与求积sum,prod4.累加和与累乘积cumsum,cumprod5.求标准差与相关系数std, corrcoef6.排序sort1.最大元素与最小元素y=max(X)[y,l]=max(X)X是向量，第一个返回向量最大元素；第二个除返回最大元素外，还将最大值元素存入k中。max(A)[

2021-08-20 12:45:05 679

原创 MATLAB----图形修饰处理

参考于：中国大学慕课，专题四“图形修饰处理”视点处理色彩处理裁剪处理1.视点处理1.1 方位角和仰角方位角又称作旋转角，是视点与原点的连线在xy平面上的投影，与y轴负方向形成的夹角，正值表示逆时针，负值表示顺时针仰角又称视角，是视点与原点的连线与xy平面的夹角，正值表示视点在xy平面上方，负值表示视点在xy平面下方1.2 视点view函数view(az,el)其中，az是方位角，el是仰角，默认值方位角-37.5°，仰角30°。[X,Y]=meshgrid

2021-08-19 11:10:55 498

原创 MATLAB----绘制三维曲面

本篇参考于：中国大学慕课科学计算与MATLAB统计，专题四“4.5三维曲面”绘制格网meshgrid格网和曲面surf&mesh其他形式surfc,surfl&meshc,meshz标准三维曲面函数sphere&cylinder三个方程两个变量使用fsurf和fmesh1.meshgrid函数生成网格eg.绘制空间曲线x=2:6y=(3:8)'[X,Y]=meshgrid(x,y)Z=randn(size(X))plot3(X,Y,Z)

2021-08-18 17:47:12 10143

原创 MATLAB----绘制三维曲线

参考于:中国大学慕课科学计算于MATLAB语言专题四“4.4三维曲线”1.plot3函数plot3(x,y,z,选项)plot3用来绘制三维曲线，与plot用法类似。当x，y，z为同型矩阵时，以x，y，z对应列绘制曲线，曲线条数为列的数量；当x，y，z有向量也有矩阵时，向量的长度应该与矩阵相符。%同型矩阵t=0:0.01:2*pi t=t'x=[t,t,t,t...

2021-08-18 15:28:23 12626 1

原创 R语言----绘制点与点连线图

引用于：公众号“R语言和统计”，原标题为“这篇《新英格兰医学杂志》论文中的美丽的图到底是怎么画的？附R代码‘> mydata<-data.frame(+ ID=c(rep(1:10,2),rep(11:20,2)),+ treatment=rep(c("v","p"),each=20),+ time=rep(c("before","after","before1","after1"),each=10),+values=c(rnorm(10,mean=100,sd=30),rnorm(

2021-08-17 18:48:50 6960 2

原创 MATLAB----其他形式的二维曲线

1.对数坐标图（用于更清楚的描述数据）（1）semilogx(x1,y1,选项1,x2,y2,选项2...)x轴为常用对数刻度，y轴为线性刻度（2）semilogy(x1,y1,选项1,x2,y2,选项2...) y轴为常用对数刻度，x轴为线性刻度（3）loglog(x1,y1,选项1,x2,y2,选项2...) x轴和y轴都为常用对数刻度eg. 绘制1/x的直角线性坐标图和三种对数坐标图subplot(2,2,1)x=0:0.1:10y=1./xplot(x,y)ti...

2021-08-17 11:10:54 252

原创 MATLAB-----绘制二维曲线

不务正业就爱做图，还是菜鸡，请多指教参考了中国大学慕课”科学计算与MATLAB语言“第4.2节1.不同频率的正弦图x=linspace(0,2*pi,100) %产生0-2pi 100个向量plot(x,[sin(x);sin(2*x);sin(3*x)]) %画title({'\bf 不同频率的正弦曲线'}) %加粗加标题legend('sin(x)','sin(2x)','sin(3x)') %加图例，与plot的顺序一致xlabel('Variable X') %x标签yla..

2021-08-16 20:52:47 2448

原创全是搬运的知识点，谢谢各位作者

0.ggplot常用函数及参数介绍引用于：ggplot常用参数汇总_R语言中文社区-CSDN博客1.ggplot2的theme()函数引用于：R语言学习实战——解决边际分布图_深竹清风的博客-CSDN博客2.RcolorBrewer包的使用引用于：R语言绘图配色-RcolorBrewer包学习及进阶 - 简书 (jianshu.com)谢谢各位大佬分享~...

2021-08-12 21:35:29 60

原创 R语言----使用ggExtra包进行边际绘图

气死我了！这个灰色怎么删掉！！！#导入基本包，RColorBrewer为颜色包，ggExtra为绘制边际图library(ggplot2)library(RColorBrewer)library(ggExtra)#设置颜色条中Set2中的五个颜色，对应5个MonthmyColors<-c(brewer.pal(5,"Set2"))ggplot(airquality,aes(Wind,Temp,col=factor(Month)))+geom_point(size=2,al.

2021-08-12 21:23:24 1890

原创 2021-08-12

制作二维散点核密度图（适合数据连续变化）theme_set(theme_bw())gg<-ggplot(airquality,aes(Wind,Temp))+geom_point(aes(col=Month),size=2,alpha=0.5)+geom_smooth(method="loess",se=F)+labs(title="Wind and Temp",subtitle="scatterplot",x="Wind",y="Temp",caption="Source:air

2021-08-12 19:02:52 61

原创 R语言----制作散点图并添加拟合线

本此练习使用的是R语言自带数据集cars~library(ggplot2)ggplot(cars,aes(speed,dist))+#使用geom_point()函数绘制散点图，其中alpha为设置透明度geom_point(size=4,color="steelblue2",alpha=0.5)+#使用geom_smooth()函数绘制拟合线，其中lm为线性拟合，se设置为FALSE取消置信区间geom_smooth(color="steelblue2",method="lm"

2021-08-10 16:20:56 15058

原创 R语言----制作数据分布图（直方图+概率密度曲线）

R语言--制作group1和group2的分布图并输出tiff格式300dpi保存制作分布图我选择的是制作直方图+概率密度曲线~ggplot(data,aes(value))+#使用facet_grid()函数以group为类左右排列facet_grid(group~.)+#使用geom_histogram()制作密度分布图为aes(x=value,y=..density..)数量分布图为y=..count..geom_histogram(bins=30,aes(x=value,

2021-08-07 22:37:40 12076 1

sqyao_147的博客