weixin_47616793的博客

import numpy as npimport open3dimport mayavi.mlabmesh=[[[0]*3]*10000]*14FDI=[17,16,15,14,13,12,11,21,22,23,24,25,26,27];#for i in range(14):#mesh[i]=open3d.io.read_triangle_mesh(‘C:/Users/lixiaohu/Desktop/auto切牙数据/1234567/150451-0/tooth_0_’+str(FDI[i

import numpy as npimport open3dfrom vispy import ioimport numpy as npimport osimport mayavi.mlabmesh=[[[0]*3]*10000]*14FDI=[17,16,15,14,13,12,11,21,22,23,24,25,26,27];for i in range(14):mesh[i]=open3d.io.read_triangle_mesh(‘C:/Users/lixiaohu/Deskt

这C:\Users\lixiaohu>cd DesktopC:\Users\lixiaohu\Desktop>dir *.zip驱动器 C 中的卷没有标签。卷的序列号是 1A68-AC22C:\Users\lixiaohu\Desktop 的目录2021/10/25 09:25 8,630,059 357176_test.zip2021/05/17 17:07 32,989,448 Netch-完整版.zip2021/07/16 10:32

clear alla=[1,0,1,1,0,1,0,1,1,1,0];k=4n=length(a(????); %求输入数组长度h=zeros(1,n-k); %解码和译码的n含义不同check=zeros(1,k); %四个组的验证数组che=0;if k4if a(1)~=rem(a(3)+a(5)+a(7)+a(9)+a(11),2); %四个组，如果有任意一个组验证不通过，则相应数组的变量变为1，标志位che变成1，方便后面修改

这个图片是一个快速确定归到哪一个类别的汉明码的归类图。汉明码编码时要再源码上插入校验位，在1，2，4，8.，16…等2的幂位置上插入即可。现在问题来了。插了这么多空，怎么填值进去呢？谁对这个值产生贡献呢？归类图就解决这样的问题，例如对9号位置，只需要查表，得到它归为P1和P8类，所以它会为该两位的取值做贡献。怎么填写贡献呢？利用奇偶校验，把归类为该位上的所有数字加起来，根据奇偶性给0或1.就是填写方法。...

cleara=[1,0,1,0,1,1,1];k=4n=length(a(????);z=1;j=1;h=zeros(1,n+k);for i=1:n+k %将输入数组写入全零矩阵待用，且跳过验证if i==2^(z-1) %在2的次幂位置上添加索引，用来填写增加的校验码作为索引的内容z=z+1;continueelseh(i)=a(j) %所以h就是一个膨胀了的数组a，加入了几个校验码j=j+1;endendh =

n=7;k=4;p=cyclpoly(n,k,‘all’);[H,G]=cyclgen(n,p(1,:));a=G(1,:)b=G(2,:)c=G(3,:)d=G(4,:)a+ba+ca+db+cb+dc+da+b+cb+c+da+c+da+b+da+b+c+d运行如下： 1 0 1 1 0 0 0 1 1 1 0 1 0 0 1 1 0 0

target=V0;十万点source=V111;七万点[COEFF,Prealligned_source] = princomp(source);[COEFF,Prealligned_target] = princomp(target);% the direction of the axes is than evaluated and corrected if necesarry.Maxtarget=max(Prealligned_source)-min(Prealligned_source);

syms a b c d e f A B C D E FX=[a,b;c,d;e,f];Y=[A,B;C,D;E,F];% X=[1,2;3,3;4,5]% Y=[1,2;3,4;4,2]muX = mean(X,1)muY = mean(Y,1)X0 = X - repmat(muX, 3, 1)Y0 = Y - repmat(muY, 3, 1)ssqX = sum(X0.^2,1)ssqY = sum(Y0.^2,1)ssqX = sum(ssqX)ssqY = sum(ssq

X=[1,2;3,3]Y=[1,2;3,4][a,b,c]=procrustes(X,Y)Z = c.b * Y * c.T + c.ca =-4.4409e-16b = 1 2 3 3c =T: [2x2 double]b: 0.7906c: [2x2 double]Z =1.0000 2.00003.0000 3.0000c.Tc.bc.cans =旋转矩阵0.3162 0.94870.9487 -0.3162a

V0是原始点云，V1是修剪后的点云，V0经过旋转得到V111，快要与V1才走到中途V111处就不动了，接下来精细化调整向V0。我们从V111上取点P，计算它离V0的哪一个点最近。发现是第79799个点离它最近。接下来。取V0（79799）得点PP...

泊松分布K=[1,2,3,4]random(‘poisson’,K )得到K = 1 2 3 4ans = 2 0 2 2取了四个数，为什么是四个？K=[1,2,3,4]random(‘poisson’,1,K )ans(:,:,1,1) = 2 1ans(:,:,2,1) = 2 1ans(:,:,3,1) = 0 1ans(:,:,1,2) = 1 2ans(:,:,2,2)

def f(p,q) :K = NumberField(x^3 +p*x+q, ‘a’);return K.integral_basis()[print(p,q,f(p,q)) for p in range(1,7,2) for q in range(1,7,2)]1 1 [1, a, a^2]1 3 [1, a, a^2]1 5 [1, a, a^2]3 1 [1, a, a^2]3 3 [1, a, a^2]35[1/3∗a2+1/3∗a+1/3,a,a2]3 5 [1/3*a^2 +

def f®:d=[0]*200e=[0]*200d[0]=ceil®;e[0]=1/(d[0]-r);for i in range(1,8):d[i]=ceil(e[i-1]);e[i]=1/(d[i]-e[i-1]);M=MatrixSpace(ZZ,8,1)d=M(d)return transpose(d)f(N(sqrt(163)))二次域得到类数与上面的连分数有关，周期为n1，。。。，nr时，其和/3-r为二次域的类数，在p为模4余3时成立def f§:K.=Quadr