该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

上了MATLAB,可是对数学完全不感兴趣。老师弄了个作业,代码附件了。求大佬帮忙理解原理后简单做一个不一样的过程的程序。帮忙注释一下代码。(一开始忘记下载要币了。不好意思)两个.m文件,需要放一个路径好像。

---------主------------------------

%Project1.m

n = 20;

x = 0:0.2:n*pi;

y = sin(x);

A =[x;y];

thta = pi/4;

rot = [cos(thta) -sin(thta);sin(thta) cos(thta)];

%R = (x.^2+y.^2).^(0.5);

B = rot*A;

% x = x1*cos(thta) + y1*sin(thta); y=-x1*sin(thta)+y1*cos(thta);

% Updated equation is follow

% -x1*sin(thta)+y1*cos(thta) = sin(x1*cos(thta) + y1*sin(thta));

% the equation of line is y = tan(thta) * x;

% the solution of root is x = k*pi/(cos(thta)+sin(thta)*tan(thta));

alpha = 1/(cos(thta)+sin(thta)*tan(thta));

p_n = 0:n;

x_1 = alpha*pi*p_n;

y_1 = tan(thta)*x_1;

figure(1);hold on;plot(B(1,:),B(2,:));

plot([0 max(x_1)],[0 max(y_1)]);

plot(x_1,y_1,'*');

hold off;

%

f = @(x)sin(x/alpha);

df = @(x)cos(x/alpha)/alpha;

init_x = x_1+1/4;

n_1 = length(init_x);

solution_x = [];

for i = 1:n_1

[xx] = mnewton(f,df,init_x(i),1e-5);

solution_x = [solution_x xx];

end

solution_y = tan(thta)*solution_x;

figure(2);hold on;plot(B(1,:),B(2,:));

plot([0 max(x_1)],[0 max(y_1)]);

plot(solution_x,solution_y,'o');

hold off;

-------副被调用---------------------

%程序2.4--mnewton.m

function [x,k]=mnewton(f,df,x0,ep,N)

%用途:用牛顿法求解非线性方程f(x)=0

%格式:[x,k]=mnewton(f,df,x0,ep,N) f和df分别为表示f(x)

%及其导数, x0为迭代初值, ep为精度(默认1e-4), N为最大迭代

%次数(默认为500), x,k分别返回近似根和迭代次数

if nargin<5,N=500;end

if nargin<4,ep=1e-4;end

k=0;

while k

x=x0-feval(f,x0)/feval(df,x0);

if abs(x-x0)

break;

end

x0=x; k=k+1;

end

matlab排课方案,求助matlab大神,学校的课程安排太骚了,我们压根就不用学matlab...相关推荐

  1. MATLAB课程表分配问题,求助matlab大神,学校的课程安排太骚了,我们压根就不用学matlab...

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 上了MATLAB,可是对数学完全不感兴趣.老师弄了个作业,代码附件了.求大佬帮忙理解原理后简单做一个不一样的过程的程序.帮忙注释一下代码.(一开始忘记下载 ...

  2. 用matlab画一些骚东西,求助matlab大神,学校的课程安排太骚了,我们压根就不用学matlab...

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 上了MATLAB,可是对数学完全不感兴趣.老师弄了个作业,代码附件了.求大佬帮忙理解原理后简单做一个不一样的过程的程序.帮忙注释一下代码.(一开始忘记下载 ...

  3. ubuntu14.04扩展屏幕后,打开matlab就死机。求大神帮忙解决!!!万分感谢

    ubuntu14.04扩展屏幕后,打开matlab就死机.求大神帮忙解决!!!万分感谢

  4. matlab里exist是什么意思,通达信exist函数用法,exist()函数在matlab中怎么用?求教大神~...

    Q1:exist()函数在matlab中怎么用?求教大神~ exist name 等价于 r=exist(name) ,在程序里面这样更加实用 0 不存在则返回值 1 name 可以是变量名,如果存在 ...

  5. 数据库安全小白紧急求助论坛大神,困扰许久

    数据库安全小白紧急求助论坛大神,困扰许久 mircosoft SQL server 2008 R2怎样才能逃过甚至杜绝外来IP不断尝试SA密码? 求助各位论坛大神,封闭1433端口,设置数据库访问策略 ...

  6. matlab运行支持向量机不出f,求助各位大神关于libsvm,svmpredict总是出不来结果,调试了好久还是不行...

    本帖最后由 maxence 于 2013-8-24 17:55 编辑 请教版各位大神: 小弟使用的是MATLAB2009,使用的是libsvm工具箱,想用SVM方法进行分类. 我的训练数据train_ ...

  7. matlab 排课,Matlab 遗传算法解决智能排课算法 一天四节课,上午两节,下午两

    Matlab 遗传算法解决智能排课算法 一天四节课,上午两节,下午两 Matlab 遗传算法解决智能排课算法 一天四节课,上午两节,下午两节,同一门课不能相邻,特殊课程不能相邻(语文和英语,数学和科学 ...

  8. matlab call lapack,科学网—fortran lapack的dgesvd和matlab SVD结果不一致,请大神指点! - 周锋的博文...

    对于A矩阵,matlab SVD分解结果如下: A = 2.2700   -1.5400    1.1500   -1.9400 0.2800   -1.6700    0.9400   -0.780 ...

  9. matlab 连续两个if,求大神解答一个matlab中的for循环嵌套if选择语句

    公告: 为响应国家净网行动,部分内容已经删除,感谢读者理解. 话题:求大神解答一个matlab中的for循环嵌套if选择语句回答:1.最后缺少一个end;2.这句输出没什么意义:disp A(i)=A ...

  10. 新手求助:大神们帮帮我,关于在ViewPage中添加GridView的问题

    如图我在viewpage中建立了这样三个页面,现在我想在第一个页面中添加一个Gridview布局,但是总是出错,GridView和Viewpage两份代码分开使用没有任何问题,但是一但嵌套在一起就老是 ...

最新文章

  1. OpenMPI 安装
  2. java 简单实现上传 | commons-fileupload
  3. Leetcode1684. 统计一致字符串的数目[C++题解]:字符串O(n^2)简单题
  4. LeetCode_111.二叉树的最小深度
  5. POJ - 2689 Prime Distance(素数区间筛模板)
  6. C#开发 —— 基础知识
  7. 如期而至,GCC 4.9.0正式版发布!
  8. mybatis获取表名——mybatis动态调用表名和字段名#{},${}
  9. 看看80万程序员怎么评论:前端程序员会不会失业?
  10. LCD屏参:手把手教你计算LCD屏参(proch值)
  11. oracle是否启用dataguard,启动和关闭data guard的步骤
  12. 设计模式之美 - 64状态模式
  13. atitit 碎片化学习.docx attilax 总结
  14. 什么是Tableau(BI工具)
  15. 高新技术企业认定专项审计报告包含报告?专审收费标准
  16. 在谷歌浏览器上面安装AxueRP插件。
  17. 《小目标目标检测的解决方法及方式》
  18. MySQL复制+高可用技术
  19. 掌握这几个关键点,玩转IC芯片数据手册和检测攻略
  20. 首个数字银行卡明年发行,广州出台区块链措施支持大湾区

热门文章

  1. OpenCV_(Fit Line with points)用直线拟合一组点
  2. opencv图像分析与处理(16)- 图像压缩中的编码方法:LZW编码
  3. 重磅直播|结构光3D相机光机核心技术及3D成像性能分析
  4. 关于u8, u16, u32, u64
  5. Python标准库映射类型与可散列数据类型的关系
  6. java 执行 jar 包中的 main 方法
  7. python面向对象(其二)
  8. 作业题---设计一个类,求和的方法,求乘积的方法
  9. jqGrid添加详细按钮,单击弹出窗体
  10. 16位伪指令汇编程序查看内存