三维柱形效果图源码程序 - matlab其它源码 - 谷速源码
下载频道> 资源分类> matlab源码> 其它源码> 三维柱形效果图源码程序

标题:三维柱形效果图源码程序
分享到:

所属分类: 其它源码 资源类型: 文件大小: 1.37 KB 上传时间: 2016-01-24 22:52:24 下载次数: 6 资源积分:1分 提 供 者: matlab源码 三维柱形效果图源码程序
内容:
三维柱形效果图源码程序,程序员在编程的过程中可以参考学习使用,希望对IT程序员有用,此源码程序简单易懂、方便阅读,有很好的学习价值!
部分代码如下:
%二维的,以第二种方法抽样
clc;
clear;
close all;
warning off;
 
 
mu=[0,0];
sigema=[1,0.3;0.3,1];
T=5000;
prosigema=1;
thetamin=[-3,-3];
thetamax=[3,3];
seed=1;rand('state',seed);randn('state',seed);
state=zeros(2,T);
theta1=unifrnd(thetamin(1),thetamax(1));
theta2=unifrnd(thetamin(2),thetamax(2));
t=1;
state(1,t)=theta1;
state(2,t)=theta2;
%===run==============
while t<T 
    t=t+1;
    new_theta1=normrnd(theta1,prosigema);
    pratio=mvnpdf([new_theta1,theta2],mu,sigema)/...
        mvnpdf([theta1,theta2],mu,sigema);
    alpha=min([1,pratio]);
    u=rand;
    if u<alpha
        theta1=new_theta1;
    end
    
    new_theta2=normrnd(theta2,prosigema);
    pratio=mvnpdf([theta1,new_theta2],mu,sigema)/...
        mvnpdf([theta1,theta2],mu,sigema);
    alpha=min([1,pratio]);
    u=rand;
    if u<alpha
        theta2=new_theta2;
    end
    state(1,t)=theta1;
    state(2,t)=theta2;
end
 
%=====figure=======
figure(1);clf;
subplot(1,2,1);
nbins=12;
thetabins1=linspace(thetamin(1),thetamax(1),nbins);
thetabins2=linspace(thetamin(1),thetamax(2),nbins);
counts=hist3(state','Edges',{thetabins1,thetabins2});
hist3(state','Edges',{thetabins1,thetabins2});
%bar3(counts,[thetabins1,thetabins2]);
xlabel('\theta1');ylabel('\theta2');zlabel('\counts');
az=61;el=30;view(az,el);
 
%=============目标分布=========
subplot(1,2,2)
nbins=50;
thetabins1=linspace(thetamin(1),thetamax(1),nbins);
thetabins2=linspace(thetamin(2),thetamax(2),nbins);
[theta1grid,theta2grid]=meshgrid(thetabins1,thetabins2);
zgrid=mvnpdf([theta1grid(:),theta2grid(:)],mu,sigema);
zgrid=reshape(zgrid,nbins,nbins);
surf(theta1grid,theta2grid,zgrid);
xlabel('\theta_1');ylabel('\theta_2');zlabel('pdf(\theta_1,theta_2)');
view(az,el);
xlim([thetamin(1),thetamax(1)]);ylim([thetamin(2),thetamax(2)]);
figure(2)
plot(state(1,:));
figure(3)
plot(state(2,:));
    

文件列表(点击上边下载按钮,如果是垃圾文件请在下面评价差评或者投诉):

三维柱形效果图源码程序/
三维柱形效果图源码程序/NewWorks/
三维柱形效果图源码程序/NewWorks/f_n_2.m
三维柱形效果图源码程序/NewWorks/example_n_2.m

关键词: 三维 效果图 源码

Top_arrow
回到顶部
联系方式| 版权声明| 招聘信息| 广告服务| 银行汇款| 法律顾问| 兼职技术| 付款方式| 关于我们|
网站客服网站客服 程序员兼职招聘 程序员兼职招聘
沪ICP备19040327号-3
公安备案号:沪公网安备 31011802003874号
库纳格流体控制系统(上海)有限公司 版权所有
Copyright © 1999-2014, GUSUCODE.COM, All Rights Reserved