matlab基础编程实例(200高分急求matlab简单编程)
本文目录
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);
用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;
更多文章:
dnf格林赛罗斯战场在哪 怎么进去?DNF战场是什么的有奖励么奖什么
2024年6月5日 02:58
魔兽争霸rpg单人地图推荐(魔兽RPG里面哪些地图可以单通好玩的)
2024年8月12日 09:51
各位帮帮忙啊,哪位能告诉我现代搜索引擎,比如百度,搜狐,新浪的缺点和优点啊!?请熟知网易的人回答,网易以前的搜索引擎是什么
2024年3月6日 11:40
穿越火线枪战王者官网(穿越火线枪战王者体验服官网地址是什么)
2024年6月18日 05:58
搜狗手游平台(搜狗游戏中心无法注册和登陆,换了几个浏览器了,手机上的uc就可以,自带的浏览器也不行,求助!!)
2024年7月15日 06:27
笔趣阁无弹窗免费阅读网下载(笔趣阁app下载的小说在哪个文件夹里能找到手机文件)
2024年6月28日 20:57
上帝之城监狱帝国下载(上帝之城监狱帝国存档修改攻略 上帝之城监狱帝国怎么修改存档)
2024年7月16日 04:36
淘宝购物时,鼠标放到淘宝主图上会出自动出现现相似宝贝,这个是什么怎么弄的呢?淘宝找同款功能取消了吗
2024年6月26日 03:42
模拟人生免费版下载(怎样下载模拟人生2,请把步骤说清楚谢谢)
2024年6月27日 03:35
eclipse安装教程详细教程(eclipse安装教程 大家一起来看看吧)
2024年8月8日 17:21
环球时报电子版在线阅读(怎样在网上阅读《环球时报》这份报纸)
2024年7月19日 05:50