数学建模社区-数学中国 › 标签 › 程序

标签: 程序

相关帖子	版块	作者	回复/查看	最后发表

求各类经典算法的MATLAB程序	MATLAB论坛	醒过沉沦 2011-8-12	22 9960	2027507950 2018-1-25 20:03
乘公交，看奥运程序	全国研究生数学建模竞赛(GMCM)	beau1988 2011-8-22	29 10876	一束阳光 2015-9-5 23:17
请求一个整数规划的matlab程序	最优化算法	huiwen 2011-8-25	11 11259	一束阳光 2015-9-5 23:19
Prims求最小生成树程序	待定区	嘿jlpz 2011-8-26	4 5262	wanglei8328 2018-2-3 19:25
MATLAB GUI界面设计	全国大学生数学建模竞赛(CUMCM)	qinghaiyin 2011-9-2	7 615	alair006 2012-2-7 16:28
数学建模必用matlab程序	Matlab 资料库视频教程讲义代码	wenxinzi 2011-9-6	55 22532	630785319 2018-2-2 15:51
一个灰色的程序	全国大学生数学建模竞赛(CUMCM)	yyt0228 2011-9-8	13 6505	qq_1502348740 2018-1-26 23:09
一个灰色的程序	全国大学生数学建模竞赛(CUMCM)	yyt0228 2011-9-8	17 7240	2027507950 2018-1-25 20:06
求助：用matlab求函数的最小值的程序	MATLAB论坛	棒棒 2011-9-11	7 8888	一束阳光 2015-9-5 23:18
board数的打分程序	Matlab 资料库视频教程讲义代码	jt202010 2011-9-17	0 1462	jt202010 2011-9-17 20:51
matlab制作动画的原程序	Matlab 图像处理视频流处理	我能 2011-9-22	169 27956	MarvinJ 2017-12-4 13:10
高斯主列消元的可执行程序（用于32位系统）	Matlab GUIDE	hahaysw 2011-10-4	0 3641	hahaysw 2011-10-4 23:32
这儿的一切，你懂的（ww）	比赛感想	Dianer 2011-10-11	12 9095	jiquan 2016-11-20 08:57
Mathematica 24点程序！	Mathematica论坛	cjilyy 2011-10-14	0 3916	cjilyy 2011-10-14 09:45
求sas小程序！急！	SAS论坛	祸劫 2011-10-18	5 3337	071029 2012-9-14 13:48
lingo程序出错，请大虾指点迷经	LINDO\LINGO论坛	¤毛毛虫¤ 2011-10-25	5 2826	时光的向日葵 2013-7-12 11:23
缄默中一个非常重要的——四阶龙格库塔算法极其实现程序	美国大学生数学建模竞赛(MCM/ICM)	zhipenglv 2011-10-28	5 1847	zhipenglv 2011-11-17 22:30
各位matlab高手们，这个程序应该怎么编呢？	MATLAB论坛	我是个白痴 2011-11-4	2 1243	我是个白痴 2011-11-4 23:26
蚁群算法程序	MATLAB论坛	haoren1990 2011-11-6	13 5859	591804904 2017-4-22 21:38
怎样用matlab生成1、X^1、X^2、…、X^n的行矩阵？	数模问题互助	bb556 2011-11-8	2 5741	851354452 2011-11-29 22:57

更多...

相关日志

分享层次分析法程序: 毕洋洋 2016-4-13 17:17; 程序： MATLAB program for AHP clc,clear fid=fopen( 'txt.txt' , 'r' ); n1=6;n2=4; a= ; % Read criterion layer judgment matrix end for i=1:n1 str1=char( ;' ]); str2=char( ;' ]); eval(str1); for j=1:n2 tmp=str2num(fgetl(fid)); eval(str2); % Judgment matrix of reading scheme layer end end ri= ; % Consistency index =eig(a); lamda=max(diag(y)); num=find(diag(y)==lamda); w0=x(:,num)/sum(x(:,num)); cr0=(lamda-n1)/(n1-1)/ri(n1) for i=1:n1 =eig(eval(char( ))); lamda=max(diag(y)); num=find(diag(y)==lamda); w1(:,i)=x(:,num)/sum(x(:,num)); cr1(i)=(lamda-n2)/(n2-1)/ri(n2); end cr1, ts=w1*w0, cr=cr1*w0,w0,w1; 213 次阅读|0 个评论

分享那位大神看到，来帮忙编一个非线性方程组解的程序: 失梦人 2014-4-6 09:10; 我要求一个非线性方程组，sqrt((x+a)^2+y^2+z^2)-sqrt(x^2+y^2+z^2)=340*t1 sqrt((x+a)^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t2 sqrt(x^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t3 要求输入参数为a，t1，t2，t3，变量。程序返回方程组的解x，y，z。 functionF=myfun(x) F= ; x0= ; options=optimset('Display','iter'); =fsolve(@myfun,x0,options) x=x(1);y=x(2);z=x(3); 假如t1,2,3和a都是已知数了，那么程序可以直接输出x，y，z。但我如何让他们以输入参数的形式出现呢。就比如说，他们是实际测量得到的数据，数据量很大的时候，总不能每次都重新输一遍程序，假如t1,2,3和a是可以matla独处的，那么在哪里赋值，怎样为之赋值才能让他们以方程组自变量的形式出现呢; 157 次阅读|0 个评论

分享各位前辈，帮小弟看个程序，关于MH抽样的，不知道什么问题，算不下去: kuai5237 2014-4-5 20:19; a= ; w=a(1:500,1); n=500; pdf=@(x)power(2*pi*8.49*8.49,-n/2)*exp(-0.5*sum(power((w-x)/8.49,2)))*((x-50)/50); proppdf=@(x,y)unifpdf(y-x,-1000,1000);%proposalpdfq(y|x) proprnd=@(x)unifrnd(50,100);%therangeis2timesoftheonesetbytheminimunandmaximummeasuredvalues nsamples=15000; =mhsample(78,nsamples,'pdf',pdf,'proprnd',proprnd,'proppdf',proppdf,'thin',21,'burnin',2000);; 个人分类: MH抽样|111 次阅读|1 个评论

分享作业一···: zuster 2012-9-25 17:25; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%以下程序均由一介遗石编写，并在Matlab R2010b成功运行%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %此程序通过使用diag函数打印矩阵 %清空环境变量 clc; clear all; %预制内存空间 B=zeros(5); B=zeros(5)+diag( ); B(3:5,1:3)=B(3:5,1:3)+diag( ); B(2:5,1:4)=B(2:5,1:4)+diag( ); B(1:4,2:5)=B(1:4,2:5)+diag( ); B(1:3,3:5)=B(1:3,3:5)+diag( ) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %此程序用rand创建一个5*6的数组,寻找该组中大于0的元素，并重新赋值为1 %清空环境变量 clc; clear all; x=rand(5,6); for i=1:5 for j=1:6 if x(i,j)0 x(i,j)=1; end end end x %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %此函数用于求解线性方程组 function x=Solving_LinearEquations(a,b) %By YanHong Song %2012年9月22日 %a为待求线性方程组的系数矩阵 %b为待求非齐次线性方程组的常数矩阵，为列向量 %调用格式：x=Solving_LinearEquations(a,b) c= ; %c为待求方程组的增广矩阵 n=size(a,1); m1=rank(a); %分别求矩阵a和矩阵b的秩，确定解的类型 m2=rank(c); if m1==m2 if m1==n %Method1：矩阵除法 disp('原方程的特解是：') x1=a\b %Method2：Gauss消去分解法 =lu(a); x2=u\(l\b) %Method3:Cholesky分解法(此法要求待求方程组的系数矩阵为对称正定矩阵) =chol(a); x3=r2\(r1\b) %Method4:QR分解 =qr(a); x4=R\(Q\b) else disp('原方程的一个特解为：') x0=a\b disp('原方程的基础解系为：') xx=null(a) end else disp('原方程的无解！') end %============================================ %题三调用格式： %============================================ clc;clear all; a= ; b= '; x=Solving_LinearEquations(a,b) %============================================ 结果得到：原方程的特解是： x1 = 2.5353 1.2033 6.8299 -0.1452 x2 = 2.5353 1.2033 6.8299 -0.1452 x3 = 2.8284 2.1213 2.4749 1.7678 x4 = 2.5353 1.2033 6.8299 -0.1452 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %此程序分别用for循环和while循环求解题四方程 %method 1：for循环 sum=0; for i=1:10 sum=sum+i^i; end sum %method 2:while循环 sum=0; i=1; while i11 sum=sum+i^i; i=i+1; end sum 两种方法得到的结果为：sum =1.0405e+010; 810 次阅读|0 个评论

分享仿真matlab程序: 紫藤恋晴天 2012-8-15 22:23; 这俩天的程序我有点抓狂，刷的题我也很无耐，谁有matlab软件的看到后，帮忙运行一下，最后的一个 “3、解调信号及频谱”中的频谱没弄出来，希望能帮忙的友人帮忙看看哪有错误吧问题二： 2 、单边带调制与解调：载波频率： 250KHZ ；要求： 1、原始信号及频谱 2、调制信号及频谱 3、解调信号及频谱 4、要求在一幅图形中给出所用结果程序： dt=0.001; fs=1/dt; t=-0.2:dt:0.2; fc=250; m=sinc(200*t); m=m.*m; subplot(421); plot(t,m); axis( ); title( ' 原始信号 ' ); c=cos(2*pi*fc.*t); d=sin(2*pi*fc.*t); u=(m.*c-imag(hilbert(m)).*d)/2; v=(m.*c+imag(hilbert(m)).*d)/2; subplot(422); plot(t,u, 'r' ); axis( ); title( ' 上边带调制信号 ' ); subplot(423); plot(t,v); title( ' 下边带调制信号 ' ); axis( ); F_m=dt*fft(m); F_u=dt*fft(u); F_v=dt*fft(v); f=linspace(-fs/2,fs/2,length(m)); subplot(424); plot(f,abs(fftshift(F_m))); title( ' 原始信号频谱 ' ); subplot(425); plot(f,abs(fftshift(F_u))); title( ' 上边带调制信号频谱 ' ); subplot(426); plot(f,abs(fftshift(F_v))); title( ' 下边带调制信号频谱 ' ); y=v.*c; F_y=dt*fft(y); subplot(427); plot(f,abs(F_y)); title( ' 解调信号频谱 ' ); H=zeros(size(f)); f_cutoff=100; l=fs/length(m); n=floor(f_cutoff/l); H(1:n)=4*ones(1,n); l1=length(m)-n; H(length(m)-n:length(m)-l)=4*ones(l,n); Z=F_y.*H; z=fs*real(ifft(Z)); subplot(428); plot(t,z); axis( ); title( ' 解调信号 ' ); 图形： dt=0.001; fs=1/dt; t=-0.2:dt:0.2; fc=250; m=sinc(200*t); m=m.*m; subplot(421); plot(t,m); axis( ); title( ' 原始信号 ' ); c=cos(2*pi*fc.*t); d=sin(2*pi*fc.*t); u=(m.*c-imag(hilbert(m)).*d)/2; v=(m.*c+imag(hilbert(m)).*d)/2; subplot(422); plot(t,u, 'r' ); axis( ); title( ' 上边带调制信号 ' ); subplot(423); plot(t,v); title( ' 下边带调制信号 ' ); axis( ); F_m=dt*fft(m); F_u=dt*fft(u); F_v=dt*fft(v); f=linspace(-fs/2,fs/2,length(m)); subplot(424); plot(f,abs(fftshift(F_m))); title( ' 原始信号频谱 ' ); subplot(425); plot(f,abs(fftshift(F_u))); title( ' 上边带调制信号频谱 ' ); subplot(426); plot(f,abs(fftshift(F_v))); title( ' 下边带调制信号频谱 ' ); y=v.*c; F_y=dt*fft(y); subplot(427); plot(f,abs(F_y)); title( ' 解调信号频谱 ' ); H=zeros(size(f)); f_cutoff=100; l=fs/length(m); n=floor(f_cutoff/l); H(1:n)=4*ones(1,n); l1=length(m)-n; H(length(m)-n:length(m)-l)=4*ones(l,n); Z=F_y.*H; z=fs*real(ifft(Z)); subplot(428); plot(t,z); axis( ); title( ' 解调信号 ' ); （在最下面及右下角少一个频谱图形）分析：调制信号图形基本相同，但是幅度发生了变化，将幅度进行调整，及可得基本相同的调调制信号，根据由傅里叶变换得到的频谱进行分析，结果基本符合调制原理，但是由于各种噪音以及误差的影响，部分地方有失真; 486 次阅读|0 个评论

分享 B题交巡警服务平台的设置与调度两点之间距离计算matlab程序: ぐ追ょ寻っ 2011-9-10 13:26; clear clc A=xlsread('C:\Documents and Settings\Administrator\桌面\B\zuobiao'); C=xlsread('C:\Documents and Settings\Administrator\桌面\B\qizhibiaohao') ; %n=input('输入站点个数n=') n=100； DD=zeros(n,n)+inf; for i=1:n P=find(A(:,1)==i); leng=length(P); for j=1:leng k=P(j); m=C(k,2); DD(i,m)= sqrt( (A(i,2)-A(m,2))^2+ (A(i,3)-A(m,3))^2 ); DD(m,i)=DD(i,m); end end; 968 次阅读|0 个评论

分享数学建模MATLAB必备程序（一）最短路: guanqinren 2011-9-8 16:27; function =dijkstra(D,s) %Dijkstra最短路算法Matlab程序用于求从起始点s到其它各点的最短路 %D为赋权邻接矩阵 %d为s到其它各点最短路径的长度 %DD记载了最短路径生成树 =size(D); d=inf.*ones(1,m); d(1,s)=0; dd=zeros(1,m); dd(1,s)=1; y=s; DD=zeros(m,m); DD(y,y)=1; counter=1; while length(find(dd==1))m for i=1:m if dd(i)==0 d(i)=min(d(i),d(y)+D(y,i)); end end ddd=inf; for i=1:m if dd(i)==0d(i)ddd ddd=d(i); end end yy=find(d==ddd); counter=counter+1; DD(y,yy(1,1))=counter; DD(yy(1,1),y)=counter; y=yy(1,1 dd(1,y)=1; end; 个人分类: 算法|0 个评论

更多...

帐号		密码		只需要一步，快速开始		注册地址	找回密码

标签: 程序

相关帖子

相关日志

QQ

电话咨询

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务|