function y=piFour(n)
y=0;
for kk=1:n
y=y+(-1)^(kk-1)*(1/(2*n-1));
end
end
%-----以下六行保存为onePer3.m文件
function y=onePer3(n)
y=0;
for kk=1:n
y=y+1/(4^kk);
end
end
============================
命令行测试如下：
>> piFour(10000)
ans =
0.7854
>> piFour(10000)*4
ans =
3.1415
%--------------------
>> onePer3(10000)
ans =
0.3333
>> onePer3(10000)*3
ans =
1追问
括号里的等于四分之π是怎么回事
追答服了你了，学成这样子，若式子中n去无穷大时，式子的和就等于括号里的值。
这是级数展开近似，n越大，精度越高（当然是计算机不溢出的情况）
">> onePer3(10000)
ans =
0.3333
>> onePer3(10000)*3
ans =
1"
这些你是只给了N=10000的吗 ？
(10000)*3 是什么意思 啊
你还是花点心思去学学基础吧，这样学将来吃亏的只会是你自己。
这么浅显的我已经没法再简单地表述或说明了。
命令行里的onePer3（N）中括号的N可以改为任意自然数。onePer3（N）计算的是1/3的近似值，其乘3自然就是1的近似值了。
就如同另外一个piFour（N）计算的是pi/4的近似值，乘上4自然就是计算pi的近似值了，这只是为了方便校验结果而已，
比如第一个问题，
y=1;
for n=1:1:100
y=y+(-1)^n*1/(2*n+1);
end
这样，尽管matlab做为矩阵实验室，用循环函数做有点浪费资源，不过我想不出来别的什么，希望路过大侠指教~
楼主关于这种循环函数方法有不明白的欢迎追问~追问
括号里的等于四分之π是怎么回事
追答极限值啊，n为无穷大的时候
y1=dot((-1).^(n-1),1./(2*n-1));
y2=sum(1./(4.^n));
n为1000，1000时数字换了就行
matlab 编程问题
题1：基本运算 >> A=[22 33 4;1 3 31;54 2 32];B=[1 4 -3;3 2 3;-3 1 8];>> A+6*B ans =[28 57 -14; 19 15 49; 36 8 80]>> A+B*eye(3,3) %eye——单位矩阵 ans =[23 37 1; 4 5 34; 51 3 40]>> A*B...
Matlab编程问题,求大佬指教。
题主使用yalmip写求解线性规划问题的程序，运行始终无结果。为什么没有结果，分析了题主的代码，主要存在着下列几个问题：1、x变量设置与目标函数，约束条件函数不对应，可以这样约定λ11x（1），λ21x（2），λ13x（3），λ24x（4）2、目标函数的变量与约定的变量不一致，可以这样写 f=(250*x(3)...
关于matlab编程的问题
1：对于第一题的微分方程组在matlab中输入[f,g]=dsolve('Df=exp(f*sin(t)+g)','Dg=exp(g*cos(t))+f','t')答案是 Explicit solution could not be found，即没有解析解，因此只能求其近似解，也就是数值解。求解微分方程的数值解需要知道“函数的初值”“自变量的范围”，即f(0)=?,g...
一道Matlab编程题
1. 首先举一个简单的例子：求y=x^2 ，在x为[0,2]上的曲线长度。把下面的复制粘贴进MATLAB syms t x=t;y=t^2;df=@(t)(1+4*t.^2).^0.5; %%MATLAB早期版本不支持@功能 quad(df,0,1)答案ans=1.4789 2. 再回答你的问题：clc clear syms t x=sin (t);y=t^2;z=log(t...
MATLAB 编程问题:矩阵Y=X*C,已知X, Y,怎么求C
X,Y)化为行最简型，检查是否有解，如果有解，写出其唯一解或通解。而化增广矩阵化为行最简型，可以用matlab命令函数rref(X,Y).方法（2），因为系数矩阵X的逆不存在，则可以使用伪逆pinv来求近似值，即C=pinv(X)*Y,因为在实际问题中，我们常常只要有具有相对精度的近似值就可以了。
matlab编程问题
我第一次接触matlab，但老师就要求我们用它来编程了 无从下手 题目：1。设方程的根为X＝[-3,-5,-8,-9],编写MATLAB程序求他们对应的X多项式的系数。2。 设f(x)=1\/((x-2)^2+0.1) +1\/((x-3)^4+0.01),写出一个MATLAB函数程序f31.m,使得调用f1时，x可用距阵代入，得出的f（x）...
matlab编程应用题
在Matlab中分别输入上述三个矩阵，各题的答案如下：1）计算C=b'-A*X'，若C的单元中出现负数，说明相应的食物供给不足。但实际计算得到C=[300,200,150]，所以食物供给充足。2）计算for i=1:4 D(i)=min(C.\/A(:,i))end 最后算得结果分别是200,150,100,100。3）记Y=X(1:3)=[1000,...
Matlab编程问题,求大神帮忙
1、使用sort函数对矩阵元素进行排序 2、使用randint函数生成随机整数，使用find函数找到相应的矩阵元素。程序如下：程序1 x=randn(5,5)x=sort(x,2)程序2 y=randint(5,5,[1 100])y(find(mod(y,3)==0))=1 y(find(mod(y,7)==0))=1 y(find(mod(y,5)==0))=0 结果为：...
matlab编程解方程
x1=-1359.7 y1=2966.7 x2=-640 y2=2126 t=2.6049 k=(y2-y1)\/(x2-x1)b=y2-k*x2 x3=-1359.8 y3=k*x3+b AC=sqrt((y3-y1)^2+(x3-x1)^2)if AC~=t*180\/0.0254 x3=x3-0.1 y3=k*x3+b AC=sqrt((y3-y1)^2+(x3-x1)^2)end x3= y3= 根据题目的要求C只有一个...
MATLAB函数编程题目,一个图形窗口分割成两块,令x属于[0 2pi] ,在第...
x=0:pi\/100:2*pi;subplot(121);plot(x,sin(x),'k-',x,cos(x),'b--');title('自己写')legend('第一土里','第二图例',0);XLabel('横坐标');YLabel('纵坐标');subplot(122)plot(x,0.5*(sin(x)+cos(x)));title('自己写')legend('第一土里','第二图例',0);XLabel('...
