[mathematical modeling] ship three degree of freedom MMG model based on MATLAB [including Matlab source code 1925]

One 、 How to get the code

How to get the code 1：
The complete code has been uploaded to my resources ：【 mathematical modeling 】 be based on matlab Ship three degrees of freedom MMG Model 【 contain Matlab Source code 1925 period 】

How to get the code 2：
By subscribing to Ziji Shenguang blog Paid column , With proof of payment , Private Blogger , This code is available .

remarks ：
Subscribe to Ziji Shenguang blog Paid column , Free access to 1 Copy code （ The period of validity From the Subscription Date , Valid for three days ）;

Two 、 Partial source code

function main
M = 20000000;       % twenty million  X  twenty million   matrix 
f1 = @()sinfun1(M);
f2 = @()sinfun2(M);
f3 = @()sinfun3(M);
f4 = @()sinfun4(M);
 
sprintf('f1- No pre allocation for Cycle time :%f',timeit(f1))
sprintf('f2- With preallocation for Cycle time :%f',timeit(f2))
sprintf('f1/f2 Multiple ：%f',timeit(f1)/timeit(f2))
sprintf('f3- No pre allocated vector time :%f',timeit(f3))
sprintf('f4- Time with preallocated vector :%f',timeit(f4))
sprintf('f2/f4 Multiple ：%f',timeit(f2)/timeit(f4))
sprintf('f1/f4 Multiple ：%f',timeit(f1)/timeit(f4))
 
function y = sinfun1(M)
x = 0:M - 1;
for k = 1:numel(x)
    y(k) = sin(x(k) / (100*pi));
end
 
function y = sinfun2(M)
x = 0:M-1;
y = zeros(1, numel(x));
for k = 1:numel(x)
    y(k) = sin(x(k) / (100*pi));
end

3、 ... and 、 Running results

Four 、matlab Edition and references

1 matlab edition
2014a

2 reference
[1] Li Xin .MATLAB mathematical modeling [M]. tsinghua university press .2017
[2] Wang Jian , Zhao Guosheng .MATLAB Mathematical modeling and simulation [M]. tsinghua university press .2016
[3] Yu Shengwei .MATLAB Mathematical modeling classic case practice [M]. tsinghua university press .2015