蒲丰投针matlab代码

用蒲丰投针法在计算机上计算π值,取a=4、l=3。

用蒲丰投针法在计算机上计算π值，取a=4、l=3。

C语言编写程序如下：

#include<stdio.h>

#include<math.h>

#include<stdlib.h>

#include<time.h>

int main()

{float x,theta,a,b,s,pi;

int n=0,i;

srand((int)time(NULL));

for(i=0;i<10000;i++)

{x=rand();

theta=rand();

printf("x=%f\ntheta=%f\n",x,theta);

a=4*x/(RAND_MAX+1.0);

b=3*sin(3.142*theta/(RAND_MAX+1.0)); if(a<=b)

{n++;}

}

s=n/10000.0;

pi=2.0*3/4/s;

printf("pi=%f",pi);

getchar();

}

运行结果如下图所示：

蒲丰氏投针问题的模拟过程，随机数发生器也是自编的，以供大家参考和提出建议。谢谢。（seed1和seed2最好选择3和5，为了使投针次数达到1000000，CVF进行如下设置Project->settings->link->

output，将stack allocations reserve:设为1000000000） ??

program getpi

implicit none

real,parameter::a=5,L=4,pi=3.14159 integer::n1,i,counter=0

real,allocatable::R1(:),R2(:) real::theta,x,pi1

write(*,*) 'input the size of the array:' read(*,*) n1

allocate(R1(n1)) allocate(R2(n1)) call random(n1,R1,R2) do i=1,n1

x=a*(2*R1(i)-1)

010 020 021 022 023

北京/密云/平谷/延庆 广州/增城/从化 上海 天津/静海/蓟县/宁河 重庆(渝中 高新 南岸)万州/壁山/开县 荣昌/铜梁、綦江 024 沈阳 025 南京/溧水县/高淳县 027 武汉 028 成都/双流/都江堰/崇州/新津/金堂/大邑/邛崃/ 郫县/彭州/蒲江/邛崃 029 西安 052 遵义 053 安顺 088 丽江 310 邯郸/成安县/永年县 311 石家庄/鹿泉/辛集/晋州/新乐/正定/无极/藁城 灵寿/元氏/赵县 312 保定/蠡县/高阳/雄县/高碑县/容城/安新/徐水 /定州/满城/定兴/涿州/安国/清苑 313 张家口 314 承德 315 唐山/玉田县/唐海/迁西 316 廊坊/霸州/香河县/文安/大城/固安/三河/大厂 317 沧州/黄骅/青县/沧县/南皮/泊头/任丘 318 衡水/安平县/枣强 319 邢台/宁晋/清河/南宫/新河 335 秦皇岛/昌黎/抚宁 351 太原 354 晋中 355 长治 370 商丘 371 郑州/新郑/荥阳/新密 372 安阳/安阳县 373 新乡/卫辉/辉县 374 许昌/长葛/禹州 375 平顶山 376 信阳 377 南阳 378 开封/尉

[matlab图像处理] 阈值分割

%迭代式阈值分割 otsu阈值分割 二值化 close all;%关闭所有窗口 clear;%清除变量的状态数据 clc;%清除命令行 I=imread('rice.png'); subplot(2,2,1); imshow(I);

title('1 rice的原图'); %迭代式阈值分割

zmax=max(max(I));%取出最大灰度值 zmin=min(min(I));%取出最小灰度值 tk=(zmax+zmin)/2; bcal=1; [m,n]=size(I); while(bcal)

%定义前景和背景数 iforeground=0; ibackground=0;

%定义前景和背景灰度总和 foregroundsum=0; backgroundsum=0; for i=1:m

for j=1:n

tmp=I(i,j); if(tmp>=tk)

%前景灰度值

iforeground=iforeground+1;

多元线性回归和BP神经网络及决策向量机之间的比较，个人理解：

多元线性回归：就是多个属性的线性组合，在组合时，通过不断调节每个属性的权重来使多元线性回归函数更多的适用于多个样本。

BP神经网络：通过使用最快速下降法，通过反向传播来不断调整网络中的权值和阈值，使网络的误差平方和最小。 决策向量机：它仍是对每个样本操作，使得所有样本距离最终生成的拟合曲线的间隔最小化。 算法比较：

pmm1BP目标函数： J?(?j)22j?1j

?yd权值调整：

k?1?wij????Jk?1?wij决策向量机目标函数：min1/2w^2

支持向量机(Support vector machines，SVM)与神经网络类似，都是学习型的机制，但与神经网络不同的是SVM使用的是数学方法和优化技术。

学习效率的比较：

导入数据： File->import data

参数优化常用方法：

[train_pca,test_pca] = pcaForSVM(train_data,test_data,97);//主元分析

[bestCVmse,bestc,bestg,ga_option]=gaSVMcgForRegress(train_label,train_pca);

多元线性回归和BP神经网络及决策向量机之间的比较，个人理解：

多元线性回归：就是多个属性的线性组合，在组合时，通过不断调节每个属性的权重来使多元线性回归函数更多的适用于多个样本。

BP神经网络：通过使用最快速下降法，通过反向传播来不断调整网络中的权值和阈值，使网络的误差平方和最小。 决策向量机：它仍是对每个样本操作，使得所有样本距离最终生成的拟合曲线的间隔最小化。 算法比较：

pmm1BP目标函数： J?(?j)22j?1j

?yd权值调整：

k?1?wij????Jk?1?wij决策向量机目标函数：min1/2w^2

支持向量机(Support vector machines，SVM)与神经网络类似，都是学习型的机制，但与神经网络不同的是SVM使用的是数学方法和优化技术。

学习效率的比较：

导入数据： File->import data

参数优化常用方法：

[train_pca,test_pca] = pcaForSVM(train_data,test_data,97);//主元分析

[bestCVmse,bestc,bestg,ga_option]=gaSVMcgForRegress(train_label,train_pca);

[matlab图像处理] 阈值分割

%迭代式阈值分割 otsu阈值分割 二值化 close all;%关闭所有窗口 clear;%清除变量的状态数据 clc;%清除命令行 I=imread('rice.png'); subplot(2,2,1); imshow(I);

title('1 rice的原图'); %迭代式阈值分割

zmax=max(max(I));%取出最大灰度值 zmin=min(min(I));%取出最小灰度值 tk=(zmax+zmin)/2; bcal=1; [m,n]=size(I); while(bcal)

%定义前景和背景数 iforeground=0; ibackground=0;

%定义前景和背景灰度总和 foregroundsum=0; backgroundsum=0; for i=1:m

for j=1:n

tmp=I(i,j); if(tmp>=tk)

%前景灰度值

iforeground=iforeground+1;

交织码的Matlab代码

从广义上说，交织编码属于分组编码，但由于采用存储器实现交织，它同时又具有卷积编码的记忆特征。在发端，它将分组长度L的信息序列通过交织器重排，在收端通过解交织器恢复原序列，以求使信道中突发错误变为无记忆随机独立差错。交织编码是一种信道改造技术，它通过信号设计将一个原来属于突发差错的有记忆信道改造为基本上是独立差错的随机无记忆信道。交织编码作为克服衰落信道中突发性干扰的有效方法，通常和其它用于纠正无记忆独立差错的信道编码相结合构成级联码（见下文），广泛应用于当代移动通信。

交织码通常表示为（M，N），分组长度L=MN，交织方式用M行N列的交织矩阵表示。一般，交织方式分为分组交织和卷积交织。分组交织的交织矩阵按列写入，按行读出；去交织矩阵按行写入按列读出。卷积交织的交织矩阵和去交织矩阵的写入与读出均按行进行。

本次仿真采用（7，4）卷积交织编码。 性能参数如下：

?123171152721??82241812628???交织矩阵：??

15932519137????2216104262014??分组卷积交织前序号：

[ 1 2 3 4 5 6 7 8 9 10 12 13 14 15 16 17 18 19 20 21 2

卡尔曼滤波matlab 代码

kalman滤波matlab代码

%kalman filter卡尔曼滤波 clear clc

A = [1,1;0,1]; B = [1/2,1]'; C = [1,0];

x1(1)= 100; %初始化 x2(1)= 10;

x = [x1(1),x2(1)]'; z=C*x;

P = [1,0;0,1]; Q=[2,0;0,1] R = 10; g=0.98; u=-g;

I=eye(2);

for k=2:20

xk=A*x+B*u; %KF xg1(k)=xk(1); xg2(k)=xk(2);

z(k)=C*xk+wgn(1,1,10);

P=A*P*A'+Q; %KF Kk=P*C'/(C*P*C'+R); %KF x=xk+Kk*(z(k)-C*xk); %KF x1(k)=x(1); x2(k)=x(2);

e1(k)=x1(k)-xg1(k); e2(

clear all;

clc;

tic;

X=zeros(1000,4);Y=zeros(1000,4);U=zeros(1000,1);

M=3;

a=[28 18 74 74 70 72 60 36 12 18 14 90 78 24 54 62 98 36 38 32]; b=[42 50 34 6 18 98 50 40 4 20 78 36 20 52 6 60 14 58 88 54]; Aeq=zeros(20,80);

for j=1:20

for k=1:20

if j==k

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

for k=21:40

if k==j+20

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

for k=41:60

if k==j+40

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

for k=61:80

if k==j+60

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

end

Aeq;

A=zeros(4,80);

for i=1:4

for k=1:80

if i==1

for k=1:20

A(i,k)=1;

end

end

if i==2

for k=21:40

A(i,k)