matlab基础编程实例(用MATLAB编写程序,具体如下谢谢啊!!!!)
本文目录
用MATLAB编写程序,具体如下谢谢啊!!!!
matlab自带这三个函数,你可以逐个查看其源代码。G = GCD(A,B)最大公约L=LCM(A,B)最小公倍is=ISPRIME(X)素数判定,素数时is为1,否则为0。 你需要的by函数只要把gcd和lcm拼一下就是了,自己动手吧。--------------------------------------------------------------------------------------------------------------matlab7中三个文件的位置在matlab安装目录下的\toolbox\matlab\specfun\中,这里随手贴出来吧:function c = lcm(a,b)%LCM Least common multiple.% LCM(A,B) is the least common multiple of corresponding elements of% A and B. The arrays A and B must contain positive integers% and must be the same size (or either can be scalar).%% See also GCD.% Copyright 1984-2003 The MathWorks, Inc. % $Revision: 5.10.4.1 $ $Date: 2003/05/01 20:43:49 $if any(round(a(:)) ~= a(:) | round(b(:)) ~= b(:) | a(:) 《 1 | b(:) 《 1) error(’MATLAB:lcm:InputNotPosInt’,... ’Input arguments must contain positive integers.’);endc = a.*(b./gcd(a,b));/////////////////////////////////////////////////////////////////////////////////////////////function = gcd(a,b)%GCD Greatest common divisor.% G = GCD(A,B) is the greatest common divisor of corresponding% elements of A and B. The arrays A and B must contain non-negative% integers and must be the same size (or either can be scalar).% GCD(0,0) is 0 by convention; all other GCDs are positive integers.%% = GCD(A,B) also returns C and D so that G = A.*C + B.*D.% These are useful for solving Diophantine equations and computing% Hermite transformations.%% See also LCM. % Algorithm: See Knuth Volume 2, Section 4.5.2, Algorithm X.% Author: John Gilbert, Xerox PARC% Copyright 1984-2003 The MathWorks, Inc. % $Revision: 5.14.4.1 $ $Date: 2003/05/01 20:43:46 $ % Do scalar expansion if necessaryif length(a) == 1 a = a(ones(size(b)));elseif length(b) == 1 b = b(ones(size(a)));endif ~isequal(size(a),size(b)) error(’MATLAB:gcd:InputSizeMismatch’, ’Inputs must be the same size.’)else siz = size(a); a = a(:); b = b(:);end;if ~isequal(round(a),a) || ~isequal(round(b),b) error(’MATLAB:gcd:NonIntInputs’, ’Requires integer input arguments.’)end for k = 1:length(a) u = ; v = ; while v(3) q = floor( u(3)/v(3) ); t = u - v*q; u = v; v = t; end c(k) = u(1) * sign(a(k)); d(k) = u(2) * sign(b(k)); g(k) = u(3);endc = reshape(c,siz);d = reshape(d,siz);g = reshape(g,siz);/////////////////////////////////////////////////////////////////////////////////////////////function isp = isprime(X)%ISPRIME True for prime numbers.% ISPRIME(X) is 1 for the elements of X that are prime, 0 otherwise.%% See also FACTOR, PRIMES.% Copyright 1984-2003 The MathWorks, Inc. % $Revision: 1.16.4.1 $ $Date: 2003/05/01 20:43:48 $if isempty(X), isp = false(size(X)); return, endif any(X(:) 《 0) || any(floor(X(:)) ~= X(:)) error(’MATLAB:isprime:InputNotPosInt’,... ’All entries of X must be positive integers.’); endisp = false(size(X));n = max(X(:));if n 》 2^32 error(’MATLAB:isprime:InputOutOfRange’,... ’The maximum value of X allowed is 2^32.’);endp = primes(ceil(sqrt(n)));for k = 1:numel(isp) isp(k) = all(rem(X(k), p(p《X(k))));end% p(p《1) would give an empty matrix and all() returns true.% we need to correct isp for this case.isp(X==1 | X==0)=0;
200高分急求matlab简单编程
%X(i+1)=1+Y(i)-1.4X(i)^2%Y(i+1)=0.3X(i)%X0=0,Y0=O%%%%%%%%%%%%%diedai.m%%%%%%%%%%%%function =diedai(x0,y0,n)xo=zeros(1,n);yo=zeros(1,n);for i=1:nif i==1xo(i)=x0;yo(i)=y0;elsexo(i)=1+yo(i-1)+1.4*(xo(i-1))^2;yo(i)=0.3*xo(i-1);endend%%%%%%%%%%%%%%%%%%%%%%%%%%%clearx0=0;y0=0;n=3000;=diedai(x0,y0,n);
请教matlab大神一个简单的编程例子仅有的20分全给你了
matlab中三角函数都是以弧度计量的,比如算sin(30°)应该输入sin(pi/6);
更多文章:
奥维互动地图(手机下载了奥维互动地图,请问怎样从地图上量出两个目标的距离)
2024年6月29日 08:06
红芯浏览器大家怎么看?哪一浏览器的自主可控国产字样统统不见,终于不再蒙人
2024年6月27日 15:03
手机开车模拟驾驶游戏(有什么模拟汽车的手游非常真实的建议推荐一下下载方式,谢谢)
2024年7月19日 07:40
手机版免费个人简历(手机怎么弄个人简历 手机如何弄个人简历)
2024年7月23日 12:10
破解AES 256位加密有多难?AES加密的详细过程是怎么样的
2024年2月9日 19:00
root精灵下载(oppo r9plus怎么才能root,什么root精灵、360一键root、百度一键root,)
2024年7月13日 05:02
怎么使用vray才能渲染室内的东西更快?3dmax找不到vray渲染器
2023年11月30日 18:40
植物大战僵尸2电脑版怎么下载(电脑用腾讯助手下载植物大战僵尸怎么下不了)
2024年4月20日 13:41
电脑软件管理app(电脑下载软件有哪些比较好的应用市场类的软件)
2024年8月2日 20:16
腾讯 使命召唤(腾讯宣布停运《使命召唤OL》,为什么玩家会觉得“被骗了”)
2024年6月5日 21:12