
数学建模与数学语言(Matlab/Lingo/IDL/R)
whucv
C#、WCF、并行编程
三维建模与图形图像处理(OpenGL/OSG)、VC/MFC
superliuwhu@126.com
展开
-
SAS进行主成分分析(Principal Component Analysis,PCA)
进行主成分分析主要步骤如下: 1. 指标数据标准化(SPSS软件自动执行),目的是消除不同量纲的影响; 2. 指标之间的相关性判定; 3. 确定主成分个数m; 4. 主成分Fi表达式; 5. 主成分Fi命名;影响中国蔬菜产量的若干因素,收集了1990-2002年13年数据,中国蔬菜生产、消费与贸易研究汪晓银 使原创 2012-09-17 22:49:52 · 27827 阅读 · 1 评论 -
SAS进行多元非线性回归
多元非线性回归方程重要方法是转化为线性回归方程.转化时应首先选择适合的非线性回归形式,并将其线性化。对于实际问题,首先应对原始数据进行作图或通过观察,选择适当函数进行拟合。已知1978~2006年全国GDP(y),第一产业x1、第二产业x2、工业生产总值x3、第三产业生产总值x4,请建立y对x1~x4的回归模型。 data ex;原创 2012-09-12 21:13:49 · 11326 阅读 · 1 评论 -
SAS进行多元线性回归
Pr=Probability(概率)多元线性回归yi=β0+β1*xi1+β2*xi2+……+ei;原假设Ho为βj=0,即线性回归的系数为0.通常使用Pr>|t|小于α时拒绝原假设Ho,认为系数不为0;否则接受原假设Ho,认为系数为0,系数没有通过检验。某水稻糙米含铬量的观测值x11.3711.349.670.7617.6715.9原创 2012-09-12 19:56:42 · 20644 阅读 · 1 评论 -
统计描述与SAS过程
1、统计报表TABULATE过程制表不是列出观测值,而是计算观测的分类统计量,绘制统计量的表格。对于数据汇总比较有用.PROC TABULATE DATA=数据集名; CLASS 分类变量; VAR 分析变量; TABLE 页维说明,行维说明,列维说明/统计量选项;RUN;统计名包括:N,NMISS,MEAN,STD,MIN,MAX,RANGE,SUM,US原创 2012-09-12 12:00:18 · 3286 阅读 · 0 评论 -
创建和运行 SAS 程序
SAS 程序包含两类程序步:DATA 步和 PROC 步。这两类程序步构成了所有 SAS 程序,它们既可以单独使用也可以结合使用。DATA 步通常用于创建或修改 SAS 数据集,但也可用来生成定制报表。例如,可以使用 DATA 步执行以下任务:将数据放入 SAS 数据集 计算新变量的值 检查并更正数据中的错误 通过对现有数据集取子集、合并和更新来生成新的原创 2012-09-11 12:20:20 · 8208 阅读 · 0 评论 -
matlab画曲面图和显示latex标记
一般来说学习数学建模,常用的软件有五种,分别是Matlab、 Lingo、 Mathematica、 SPSS和SAS。对于数学规划的模型,建议大家使用Lingo软件求解比较方便,对于其它问题,如时间序列模型,你使用什么软件求解都可以,关键看个人的喜好和对某种软件的熟悉程度,例如你可以使用SPSS, SAS或Eviews, R软件等求解时间序列模型。画出x^2+(y-5)^2=16绕x轴旋转一原创 2012-09-05 21:09:34 · 7191 阅读 · 0 评论 -
matlab中map工具箱(Mapping toolbox)使用
Mapping Toolbox是Matlab提供的,一整套包含许多函数跟图形用户界面的工具箱;它可以帮助用户分析几何空间方面的数据;并以地图的形式显示出来!worldmap命令help worldmap WORLDMAP 创建世界上给定区域的坐标 WORLDMAP 区域 或 WORLDMAP(区域名) 新建一个空的地图坐标,这个投影模型和限制与此区域相适应 . REGION 可能是一个字符原创 2012-09-01 22:48:28 · 49949 阅读 · 2 评论 -
LINGO调用VC编写的函数动态库实验
Visual C++编写函数,然后用动态库的形式输出该函数供LINGO调用.某公司生产A,B和C三种产品,售价分别是12元、7元和6元.生产每件A产品需要1小时技术服务、10小时直接劳动、3千克材料;生产每件B产品需要2小时技术服务、4小时直接劳动、2千克材料;生产每件C产品需要1小时技术服务、5小时直接劳动、1千克材料.现在最多能提供100小时技术服务、700小时直接劳动、 400千原创 2012-08-31 22:50:09 · 3533 阅读 · 2 评论 -
matlab中将公式以latex显示到图表
>> syms M>> syms A>> syms B>> f=dsolve('Dy*10000=M-B*y-A','y(0)=100','t') f = exp(-1/10000*B*t)*(-M+A+100*B)/B-(-M+A)/B >> latex(f)ans ={e^{-{\frac {1}{10000}}\,Bt}} \left( -M+A+100\,B原创 2012-08-14 21:29:02 · 1529 阅读 · 0 评论 -
投入产出模型分析 分配平衡方程组和消耗平衡方程组
zeros 生成0矩阵 eig 特征值、特征向量ones 生成1矩阵 diag 对角矩阵eye 生成单位矩阵 trace 方阵的迹linspace 生成等距行向量 rank 矩阵的秩rand原创 2012-08-04 17:43:15 · 7604 阅读 · 0 评论 -
simulink信号分支和多信号输入到一个示波器
分支(利用鼠标右键画线)首先在sine波形和积分器间建立一条连接线,要想分支,将鼠标移动到当前连接线间任意一点,按下鼠标右键,移动到示波器信号输入端,松开鼠标右键。示波器默认为一个输入端,怎样增加输入信号接口数量。双击示波器图标,打开Scope窗口,在工具栏第二个图标scope parameters单击,更改Number of axes的数目为n,这里是需要接受的输入信号数目。原创 2012-07-12 13:19:44 · 11097 阅读 · 0 评论 -
matlab求解线性方程组
直接解法Ax=l,求x。利用逆矩阵、矩阵除法、x=inv(A'A)*A'llind.m文件如下,%求解线性方程组% x1 + 2x2 + 3x3 + 9x4 = 5%2x1 + 2x2 + 5x3 + 4x4 = 2%3x1 + 5x2 + x3 + 5x4 = 3%7x1 + 4x2 + 2x3 -10x4 = 8A=[1 2 3 9;2 2 5 4;3 5 1 5;原创 2012-07-11 12:13:38 · 4105 阅读 · 0 评论 -
simulink计算积分
>> simulink弹出如下界面,选择File|New|Model,对满足数学规律过程进行仿真,结果用示波器表示,仿真时间t=10uint。需要积分模块、正弦模块做数据处理,示波器模块用作显示。添加正弦波、积分和scope原创 2012-07-12 12:20:17 · 13934 阅读 · 0 评论 -
matlab显示图形窗口解决方案
%matlab 2010a>>figure%后没出现图形显示窗口直接plot也没有显示%输入如下语句显示图形窗口,关掉matlab后重新打开,恢复正常可以直接plotclose allx = linspace(0,2*pi,100);y = sin(x);figure,plot(x,y);原创 2012-07-03 17:22:23 · 4709 阅读 · 0 评论 -
matlab数据拟合与线性回归
函数拟合知道y=a+b*exp(-x)观测数据x 0.0 0.3 0.8 1.1 1.6 2.3y 0.82 0.72 0.63 0.60 0.55 0.5求a、b,线性回归或拟合问题。建立矛盾方程组AX=y其中X=[a,b]'这里将a、b看做待求量,X=A\y或由最小二乘准则X=inv(A'A)*A'yt=[ 0.0 0.3 0.8 1.1 1.6原创 2012-07-03 22:22:15 · 2181 阅读 · 0 评论 -
给定n个三维坐标点,matlab求平面度
File ImportData,打开点云文件,默认变量名为文件名(可以重命名变量)//导入数据导出的38502个点,这里plane.txt//planex=plane(:,1);y=plane(:,2);z=plane(:,3);A=[x,y,ones(size(x))];B=z;abc=A\B;a=abc(1);b=abc(2);c=abc(3);F=(a*x+b*y+c-原创 2012-03-05 22:46:24 · 5258 阅读 · 2 评论