利用复数乘法几何意义绘制五角星

王福昌 2018-3-22 22:59

利用复数的乘法可完成平面上的旋转变换。每次逆时针转144°，则连转5次即可回到原点，即一笔画出五角星。用复数乘法的意义，可以很方便地计算出每个顶点的坐标，从而绘出五角星。下面给出了MATLAB源代码。 theta = 144/180*pi; n = 1:6; z = exp(i*theta*n); x = real(z); y = imag(z); fill( ...

5631 次阅读|没有评论

F分布密度曲线和上alpha分位数

王福昌 2018-2-26 23:21

clear,clc,close all; figure('color', );hold on x = 0:0.01:5; df = ; for n = 1:length(df) y = fpdf(x,df(n,1),df(n,2));plot(x,y,'linewidth',1); end text(1.22,0.6,'$$\leftarrow n_1=10,n2=40$$','interpreter','latex'); text(2,0.19,&# ...

8997 次阅读|没有评论

t 分布的密度曲线

王福昌 2018-2-26 23:11

clear,clc,close all; figure('color', );hold on x = -5:0.1:5; df = ; for n = 1:length(df) y = tpdf(x,df(n));plot(x,y,'linewidth',1); end text(-1,0.15,'$$\leftarrow n=1$$','interpreter','latex'); text(-0.8,0.25,'$$\lefta ...

3230 次阅读|没有评论

卡方分布的密度曲线

王福昌 2018-2-26 22:57

密度曲线 上alpha分位点 clear,clc,close all; ...

10319 次阅读|没有评论

修正箱线图的绘制

王福昌 2018-2-24 08:00

绘制修正箱线图的具体方法见盛骤等，《概率论与数理统计》第四版，高等教育出版社P134-135. 下面给出图形和程序：红*表示异常值。 %%% 修正箱线图 x = ;% ; xs = sort(x); n = length(x); p = ; xnp = zeros(size(p)); np = floor(n*p); r = n*p-np; id = (r ==0); xnp(id) = (xs(np(id))+xs(np(id)+1 ...

7370 次阅读|没有评论

箱线图的绘制

王福昌 2018-2-23 17:51

参考教材《概率论与数理统计》第四版， 高等教育出版社 盛骤等编 P133 程序为： clear,clc,close all x = ; xs = sort(x); %%% 绘制箱线图 n = length(x); p = ; xnp = zeros(size(p)); np = floor(n*p); r = n*p-np; id = (r ==0); xnp(id) = (xs(np(id))+xs(np(id)+1))/2; ...

3013 次阅读|没有评论

教材和MATLAB中关于求一元数据 p (0<p<1)分位数的求法对比

王福昌 2018-2-22 23:07

参考教材《概率论与数理统计》第四版， 高等教育出版社 盛骤等编 P132 根据书中定义求得的结果与根据MATLAB的quantile和prctile得到的结果有一些差别，原因是计算方法不同。我编写如下MATLAB程序进行对比，找找其中的差别吧？ 维基百科： https://en.wikipedia.org/wiki/Quantile 提供了9个公式 ...

2829 次阅读|没有评论

复变函数映射 w = z^2 的一种几何解释

王福昌 2018-2-19 00:05

图形根据教材解释绘制，贴出MATLAB代码来，供大家参考。欢迎提供宝贵意见！ 设 z = x + i y, w = u + i v, 则 w = z^2 = (x^2-y^2) + i*2*x*y 反映了两对变量 u,v 和x,y 之间的对应关系，因为涉及到 4 个变量，无法用一个平面内的几何图形表示出来， 必须把它看成两个复平面上的点集之间的 ...

12074 次阅读|没有评论

复数 1+i 的 4 次方根在复平面分布图

热度 1 王福昌 2018-2-17 23:14

准备开学上复变函数，图1.9. 图后有 MATLAB源代码，绘制箭头需下载 arrow3.m clear,clc,close all z = 1 + i; r = 2^(1/8); t = 0:0.05:6.3; plot(r*cos(t),r*sin(t),'k');hold on% 画圆 text(1,1, ,'interpreter','latex','fontsize',14); for k = 0:3 w(k+ ...

8685 次阅读|3 个评论 热度 1

本页有 1 篇博文因作者的隐私设置或未通过审核而隐藏