![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Matlab
文章平均质量分 77
asdzxcXCasd1
这个作者很懒,什么都没留下…
展开
-
Chebyshev多项式, Legendre多项式,与Chebyshev多项式零点插值
Chebyshev多项式使用通项公式function v = che_1(a,b,n)syms x;t = ((b-a)*x + a + b)/2;v = expand(cos(n*acos(t)));end使用递推关系function v = che_2(a,b,n)syms x;t = ((b-a)*x + a + b)/2;t_0 =原创 2013-10-15 19:41:53 · 8522 阅读 · 1 评论 -
数值分析第三章上机实习题
第一题没什么好说的.用polyfit拟合就可以了.clear;f=@(x) 1./(1+25*x.*x);x = -1 + 0.2*(0:10);y = f(x);P = polyfit(x,y,3);u = -1:0.01:1;v_1 = f(u);v_2 = polyval(P,u);plot(x,y,'ro',u,v_1,'k-',u,v_2,'b--');原创 2013-10-19 11:47:07 · 3965 阅读 · 1 评论 -
数值分析第四章上机实习题
function [I,h] = compoundT(fx,a,b,n)%复化梯形公式%输入fx-sym形式的函数,a/b-积分上下限,n-将[a,b]分为n等分%输出I-积分,h-步长u = linspace(a,b,n+1);h = (b-a)/n;%这个就是简单一除.......v = subs(fx,u);findnan = find(isnan(v)==1);%考虑到会有原创 2013-10-31 21:02:52 · 4789 阅读 · 0 评论 -
数值分析第七章非线性方程MATLAB程序
二分法,最简单的,貌似没要求....function [x, count] = bisection(fx,a,b,error,count)if(nargin == 4)%调用时不需要输入count,计数用,用户不需要知道 count = 1;endx = (a+b)/2;fa = subs(fx,a);f2 = subs(fx,x);err = b - a;if(er原创 2013-11-16 13:47:15 · 2054 阅读 · 1 评论 -
Google面试题:扔玻璃珠Matlab实现
题目:某幢大楼有100层。你手里有两颗一模一样的玻璃珠。当你拿着玻璃珠在某一层往下扔的时候,一定会有两个结果,玻璃珠碎了或者没碎。这幢大楼有个临界楼层。低于它的楼层,往下扔玻璃珠,玻璃珠不会碎,等于或高于它的楼层,扔下玻璃珠,玻璃珠一定会碎。玻璃珠碎了就不能再扔。现在让你设计一种方式,使得在该方式下,最坏的情况扔的次数比其他任何方式最坏的次数都少。也就是设计一种最有效的方式。假设第原创 2013-11-27 09:55:15 · 3004 阅读 · 0 评论